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Abstract 


Over the last decades, traffic congestion, car accidents and pollution became daily 
issues. To understand and overcome road traffic problems, scientists from different 
research fields are creating advanced mathematical models. Mathematical models 
help to understand road traffic phenomena, develop optimal road network with 
efficient movement of traffic and minimal traffic congestion. This thesis is devoted 
to macroscopic traffic flow modelling, which describes traffic flow by variables av- 
eraged over multiple vehicles: density, velocity and flow. Macroscopic models 
naturally lead to conservation laws, which are hyperbolic partial differential equa- 
tions. In recent years, this class of equations is more widely considered, but few 
theoretical results are available. This is caused by two main difficulties. The for- 
mer is the non-linear hyperbolic nature of equations, which leads to consider weak 
solutions, instabilities and diffusivity of numerical schemes. The latter is the non- 
uniqueness of weak solutions and the need to introduce exotic functionals to select 
a unique physically reasonable solution. 

In the first chapter, we introduce basic ideas of traffic modelling. First, we 
present the main classification of mathematical models with special attention to 
the level of details. Then we list the differences between the dynamics of traffic 
flow and that of flowing particles. Next, we show the minimal requirements to 
construct a physically reasonable macroscopic traffic flow model. We define three 
macroscopic variables to describe traffic flow, namely (average) density p, (average) 
speed v and (average) flow f. We derive the basic relation between them and 
formulate scalar conservation law. The chapter ends with a short presentation 
of the models under consideration, followed by the results obtained during my 
doctoral studies. 

The second chapter is devoted to a detailed discussion of basic macroscopic 
traffic flow models. The first presented model is the model proposed by Lighthill, 
Witham [68] and Richards [81] (LWR). It describes the dynamics of traffic via a 
scalar conservation law under the hypothesis that v = v(p). We define a rarefaction 
wave, a shock wave and a contact discontinuity for the LWR model, and define the 
Riemann solver RS wr. In the end, we give a list of drawbacks of the LWR model. 


The next considered model is the Aw, Rascle [7] and Zhang [85] model(ARZ). 
The ARZ model consists of two conservation laws, expressing the conservation of 
the number of vehicles and the conservation of the generalized momentum. We 
give the basic properties of the system, such as eigenvalues, eigenvectors and the 
corresponding Lagrangian markers. Next, we construct the Riemann solver RS Anz 
using elementary waves. Finally, we give definitions of weak and entropy solutions 
for the ARZ model corresponding to RSARz. 

In the last part of this chapter, we describe models with phase transition (PT). 
The PT model treats differently traffic with low and high densities, on the basis of 
empirical studies. For this reason, we consider PT model described by the LWR 
model on the set €); corresponding to the low densities and a 2 x 2 system of con- 
servation laws on the set 0), corresponding to the high densities. We present two 
PT models, denoted by PT“ and PT”, and introduced in [15,53]. Then we recall 
from our paper [35] the generalization of these models for cases with a metastable 
phase (Qt. Z 0) and without a metastable phase (4O Qe = 0). Next, we intro- 
duce a notion of admissible solution for the Riemann problems and then Riemann 
solvers RSR and RSs accordingly. The chapter ends with propositions regarding 
consistency and Lj,-continuity for the Riemann solvers RSR and RSs. 

In the third chapter, we describe the LWR model with a local point constraint 
on the flow. More precisely, we consider a situation in which the maximum flow 
of cars is limited at a fixed point along the road. Thanks to such considerations, 
we can model traffic flow through toll gates or construction sites. We define the 
Riemann solver CRS wr and list its main properties. Then we define the entropy 
solution of the Cauchy problem and recall the corresponding existence result. 

The fourth chapter is devoted ARZ model with local point constraint on the 
flow and our results obtained in [42]. In our work we prove the existence of the 
weak solutions, corresponding to a non-conservative Riemann solver, in the class 
of functions with bounded variation. The goal is obtained by showing the con- 
vergence of a sequence of approximate solutions constructed via the Wave Front 
Tracking method. More precisely, we introduce grid, approximate Riemann solver 
CRS ‘pz, by splitting a rarefaction wave and construct approximate Cauchy prob- 
lems. Thanks to the decreasing in time functional Y, we show that the total 
variation of the approximated solution is uniformly bounded. By Helly's theorem 
we obtain convergence of approximated solutions and then we show that the limit 
function is indeed a weak solution to the Cauchy problem for the ARZ model with 
local point constraint on the flow. 

In the fifth chapter, we describe the models PT* and PT? with the local point 
constraint on the flow and present our results obtained in [10,35]. More precisely, 
we introduce Riemann Solvers CR Sg and CRSS, both with a metastable phase and 
without a metastable phase. Then we examine their consistency, Lj, -continuity 
and invariant domains. The remainder of the chapter is devoted to the existence 
result of a weak solution in the class of function with bounded variation for the PT? 
model with a metastable phase. The goal is obtained by showing the convergence of 


a sequence of approximate solutions constructed via Wave Front Tracking method. 
Similarly to the results from the previous chapter, we define grid and approximate 
Riemann solver CRS”. Then we introduce the decreasing in time function T and 
show that the approximate solution has bounded variation, the number of waves 
and interactions is finite in finite time. We apply Helly's theorem and then show 
that the limit function is an entropy solution of the Cauchy problem for the PT* 
model with the metastable phase. 

The sixth chapter is devoted to the results obtained in conference proceedings 
[34,43]. We consider there two macroscopic models on road networks. The former 
is the LWR model with moving constraint on the flow. The concept of moving 
constraint on the flow allows us to model situations in which a truck (or other 
slower vehicle) reduces the flow at its position. From a mathematical point of 
view, the constraint is given by the ordinary differential equation depending on 
the trajectory of the truck. We give a detailed description of the model for a 
unidirectional road, introduce a Riemann solver D, Sjwmg and generalize it for the 
case of road networks. The latter considered model is the PT model introduced in 
the second chapter. We generalize it to the case of road networks by introducing 
an appropriate Riemann solver. 

At last, for the sake of clarity and to ease of comprehension, we defer to the 
appendix technical proofs. 
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Streszczenie 


Na przestrzeni ostatnich dziesięcioleci zatłoczone ulice, wypadki samochodowe oraz 
związane z ruchem samochodowym zanieczyszczenie powietrza stały się codzi- 
ennością. Naukowcy z różnych dziedzin nauki tworzą zaawansowane modele matem- 
atyczne, pomagające zrozumieć zjawiska ruchu drogowego, rozwijać efektywnie 
sieć dróg oraz zmniejszać korki uliczne. Ta rozprawa doktorska poświęcona jest 
makroskopowemu modelowaniu ruchu drogowego, które formułuje zależności między 
uśrednionymi charakterystykami przepływu ruchu takimi jak gęstość, przepływ 
oraz prędkość. Modele makroskopowe w sposób naturalny prowadzą do stosowania 
praw zachowania, które są szczególnymi równaniami różniczkowymi cząstkowymi. 
W ostatnich latach ta klasa równań jest coraz chętniej rozpatrywana, ale wciąż 
dostępnych jest niewiele wyników teoretycznych. Spowodowane jest to dwoma 
poważnymi trudnościami. Pierwszym z nich jest nieliniowa, hiperboliczna natura 
równań, mogąca prowadzić do rozważania słabych rozwiązań, niestabilności lub dy- 
fuzyjności schematów numerycznych. Drugim natomiast jest brak jednoznaczności 
słabych rozwiązań oraz potrzeba rozważania egzotycznych funkcjonałów w celu 
wybrania jednoznacznego, fizycznie uzasadnionego rozwiązania. 

W rozdziale pierwszym prezentujemy wstęp do modelowania ruchu drogowego. 
Na początku podajemy główny podział modeli matematycznych, przy zwróceniu 
uwagi na poziom szczegółowości. Podajemy różnice pomiędzy strukturą płynów a 
strukturą ruchu drogowego, oraz wymieniamy podstawowe założenia potrzebne do 
skonstruowania poprawnego makroskopowego modelu ruchu drogowego. Defini- 
ujemy trzy charakterystyki przepływu ruchu, to jest (średnią) gęstość, (średnią) 
prędkość i (średni) przepływ, wyprowadzamy podstawową zależność pomiędzy nimi 
a następnie formujemy prawo zachowania. Rozdział zakończony jest krótkim przed- 
stawieniem stosowanych modeli oraz wyników uzyskanych podczas trwania studiów 
doktoranckich. 

Drugi rozdział poświęcony został szczegółowemu omówieniu podstawowych mod- 
eli makroskopowych dla ruchu drogowego. Jako pierwszy prezentujemy model za- 
proponowany przez Lighthilla, Withama [68] oraz Richardsa [81] (LWR). Opisuje on 
dynamikę ruchu drogowego poprzez skalarne prawo zachowania wraz z warunkiem 
v = v(p). Następnie wprowadzamy pojęcia fali rozrzedzającej, fali uderzeniowej 
oraz nieciągłości kontaktowej dla modelu LWR oraz definiujemy rozwiązanie za- 
gadnienia Riemanna RSpwrR. Na koniec podajemy listę wad modelu LWR. 

Kolejny opisywany model został zaproponowany przez Aw, Rascla [7] oraz 
niezależnie Zhanga [85| (ARZ). Jest on opisany poprzez układ dwóch praw zachowa- 
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nia. Pierwsze z nich jest prawem zachowania liczby pojazdów, natomiast drugie 
jest uogólnionym prawem zachowania pędu. Na początku rozdziału podajemy pod- 
stawowe własności układu, takie jak wartości i wektory własne oraz odpowiadające 
im znaczniki Lagrange'a. Następnie definiujemy rozwiązanie zagadenienia Rie- 
manna RSARz przy pomocy podstawowych fal. Na koniec podajemy definicje 
słabego rozwiązania oraz rozwiązania entropijnego dla modelu ARZ. 

W ostatniej części tego rozdziału opisujemy modele z przemianą fazowa(PT). 
Modele PT definiują odmiennie ruch drogowy dla dróg niezatłoczonych oraz dróg 
zatłoczonych. Z tego powodu rozważamy modele PT składające się z modelu 
LWR na zbiorze €; odpowiadającym niezatłoczonym drogom oraz układu dwóch 
praw zachowania na zbiorze Q. odpowiadającym zatłoczonym drogom. Przed- 
stawiamy dwie wersje modeli PT, oznaczone przez PT^ i PT? i wprowadzone 
w [15,53]. Następnie podajemy nasze uogólnienia [35] na przypadki odpowied- 
nio bez fazy metastabilnej (Qf N Q. = 0) oraz z fazą metastabilną (Qf N Qe z 0). 
Następnie definiujemy rozwiązania dopuszczalne dla zagadnienia Riemanna oraz 
spełniające tę definicję rozwiązania zagadnienia Riemanna RSR i RSs. Rozdział 
kończymy propozycjami dotyczącymi niezmienniczości oraz ciągłości w przestrzeni 
Lj... funkcji lokalnie catkowalnych dla RSR i RSs. 

W rozdziale trzecim opisujemy model LWR ze stałym ograniczeniem na przepływ. 
Innymi słowy, rozpatrujemy sytuację, w której maksymalny przepływ samochodów 
jest ograniczony w pewnym miejscu na drodze. Dzięki takim rozważaniom możemy 
modelować ruch drogowy na rogatkach, czy w miejscach robót drogowych. Defini- 
ujemy rozwiązanie zagadnienia Riemanna CRS pwg oraz podajemy jego własności. 
Następnie definiujemy entropijne rozwiązanie zagadnienia Cauchy'ego oraz przy- 
pominamy twierdzenie o istnieniu jednoznacznego rozwiązania entropijnego zagad- 
nienia Cauchy ego. 

Rozdział czwarty został poświęcony modelowi ARZ ze stałym ograniczeniem 
na przepłym oraz wynikom własnym uzyskanym w [42]. W naszej pracy udowad- 
niamy istnienie słabego rozwiązania zagadnienia Cauchy'ego w klasie funkcji o 
wahaniu ograniczonym dla modelu ARZ ze stałym ograniczeniem na przepływ. 
W artykule rozpatrujemy rozwiązanie zagadnienia Riemanna dla którego warunek 
Rankine'a-Hugoniota nie jest spełniony dla drugiego równania. Dowód głównego 
twierdzenia opiera się na metodzie Wave Front Tracking. Mówiąc dokładniej, 
wprowadzamy siatkę i definiujemy przybliżone rozwiązanie zagadnienia Riemanna 
CRShpz poprzez podział fali rozrzedzającej. Następnie, dzięki CARS, kon- 
struujemy przybliżone rozwiązanie zagadnienia Cauchy'ego dla tego problemu. 
Przy pomocy malejącego w czasie funkcjonalu Y pokazujemy, że wahanie aproksy- 
mowanego rozwiązania jest ograniczone oraz liczba fal i interakcji między nimi 
jest skończona. Następnie, stosując twierdzenie Helly'ego otrzymujemy istnie- 
nie zbieżnego podciągu aproksymowanych rozwiązań. Ostatecznie wykazujemy, że 
funkcja graniczna jest rzeczywiście słabym rozwiązaniem zagadnienia Cauchy'ego 
dla modelu ARZ ze stałym ograniczeniem na przepływ. 

W rozdziale piątym opisujemy modele P'T* i PT? ze stałym ograniczeniem na 
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przepływ oraz podajemy wyniki własne zawarte w [10,35]. Mówiąc dokładniej, 
wprowadzamy rozwiązania zagadnienia Riemanna CRSR i CRSs. Następnie w 
badamy ich zgodność, ciągłość w przestrzeni Lj,, funkcji lokalnie całkowalnych 
oraz zbiory niezmiennicze. Pozostała część rozdziału została poświęcona wynikowi 

istnienia słabego rozwiązania entropijnego zagadnienia Cauchy'ego w klasie funkcji 

o wahaniu ograniczonym dla modelu PT” z fazą metastabilną. Dowód głównego 

twierdzenia opiera się na metodzie Wave Front Tracking. Podobnie do rozważań w 

poprzednim rozdziale, definiujemy siatkę, definiujemy przybliżone rozwiązanie za- 

gadnienia Riemanna CRS”. Następnie wprowadzamy malejący w czasie funkcjonał 
T i pokazujemy, że wahanie przybliżonego rozwiązania jest ograniczone oraz liczba 

fal i interakcji między nimi jest skończona w skończonym czasie. Stosując twierdze- 

nie Helly'ego otrzymujemy istnienie zbieznego podciagu przybliżonych rozwiązań. 

Na koniec udowadniamy, że funkcja graniczna jest słabym rozwiązaniem entropi- 

jnym zagadnienia Cauchy'ego dla modelu PT" z fazą metastabilną. 

Rozdział szósty został poświęcony wynikom otrzymanym w materiałach kon- 
ferencyjnych [34,43]. W obu pracach rozważaliśmy modele makroskopowe na 
skrzyżowaniach drogowych. Pierwszym z wprowadzonych modeli jest LWR z 
ruchomym ograniczeniem na przepływ. Ruchome ograniczenie na przepływ po- 
maga nam modelować sytuacje, w których ciężarówka bądź inny wolniejszy po- 
jazd redukuje przepływ w swoim otoczeniu. Z matematycznego punktu widzenia, 
ograniczenie jest dane przez równanie różniczkowe zwyczajne zależne od trajek- 
torii tego pojazdu. Podajemy dokładny opis modelu dla drogi jednokierunkowej, 
wprowadzamy rozwiązanie zagadnienia Riemanna BRS;wr oraz uogólniamy je 
na przypadek skrzyżowania drogowego. Drugim rozważanym rodzajem modeli 
jest model PT wprowadzony w rozdziale drugim. Uogólniamy go do przypadku 
skrzyżowania drogowego poprzez wprowadzenie odpowiedniego rozwiązania zagad- 
nienia Riemanna. 

Ostatnią część pracy, dla zachowania przejrzystości tekstu, stanowią dodatki z 
dowodami twierdzeń pomocniczych. 


Resume 


Durant ces 10 dernières années les embouteillages, les accidents de voiture, la pol- 
lution sont devenus des problèmes quotidiens. Afin de mieux comprendre et de 
mieux surmonter les problèmes de trafic routier, les scientifiques de tout domaine 
ont développé des modèles mathématiques avancés. Les modèles mathématiques 
aident à comprendre les phénomènes du trafic routier, à développer des réseaux 
de routes optimaux avec des circulations de véhicules efficaces et des problèmes 
d'embouteillage minimaux. Cette thèse est consacrée à la modélisation macro- 
scopique du trafic routier où celle-ci décrit le trafic avec des variables moyennées sur 
plusieurs voitures : densité, vitesse et flux. Les modèles macroscopiques aboutis- 
sent naturellement aux lois de conservation, qui sont des équations aux dérivées 
partielles hyperboliques. Depuis ces récentes années, cette classe d'équations a 
davantage été considérée, cependant peu de résultats théoriques sont disponibles. 
Cela est dû à deux difficultés principales. La première est la nature hyperbolique 
non linéaire de ces équations, qui conduit à considérer la notion de solutions faibles, 
les instabilités et la diffusion de schémas numériques. La deuxième difficulté est la 
non unicité de la solution faible et la nécessité d'introduire de fonctionnelles exo- 
tiques permettant de sélectionner une unique solution physiquement raisonnable. 

Dans le premier chapitre, nous introduisons les idées de base de la modélisation 
du trafic. Dans un premier temps, nous présentons les classifications principales 
des modéles mathématiques avec une attention particuliére donnée au niveau des 
détails. Nous continuons par lister les différences entre la dynamique de trafic 
routier et celle des particules. Ensuite nous montrons les besoins minimaux pour 
construire un modele de trafic macroscopique physiquement raisonnable. Nous 
définissons trois variables macroscopiques pour décrire le trafic, à savoir la densité 
(moyenne) p, la vitesse (moyenne) v et le flux moyen f. Nous tirons les relations de 
base qui les lient et nous formulons une loi de conservation scalaire. Le chapitre se 
termine par une présentation courte des modèles considérés, suivi de la présentation 
des résultats obtenus durant mes travaux de these. 

Le second chapitre est dédié à une discussion détaillée des modéles macro- 
scopiques de base dans le trafic. Le premier modele présenté est le modele proposé 
par Lightill, Witham [68] et Richards [81] (LWR). II décrit la dynamique du trafic 
par une loi de conservation scalaire sous l'hypothése v — v(p). Nous définissons ce 
que sont une onde de détente, une onde de choc et un contact de discontinuité pour 
le modéle LWR, et nous définition le solver de Riemann. À la fin du chapitre, nous 
donnant une liste de défauts du modele LWR. Le modele considéré ensuite est celui 


d'Aw, Rascle [7] et Zhang [85] (ARZ). Le modele ARZ consiste en deux lois de 
conservation, exprimant la conservation du nombre de vćhicules et la conservation 
du moment généralisé. Nous donnons les propriétés de base du système, telles que 
les valeurs, vecteurs propres et les marqueurs lagrangiens correspondants. Ensuite 
nous construisons le solveur de Riemann RS arz en utilisant les ondes élémentaires. 
Enfin nous définissons les notions de solutions faibles pour le modéle ARZ associé 
à RSARZ. 

Dans la derniere partie de ce chapitre, nous dócrivons le modele de transition 
de phase (TP). le modèle TP traite de manière différente le trafic avec de faibles 
densité et de hautes densités sur la base d’études empiriques. C’est la raison pour 
laquelle, nous considérons le modèle TP décrit par le modèle LWR sur l’ensemble 
Q; correspondant aux faibles densités et à un système 2 x 2 de lois de conserva- 
tion sur l’ensemble (2. correspondant aux hautes densités. Nous présentons deux 
modèles PT, notés PT^ et PT?, et introduits dans [15,53]. Ensuite à partir de 
notre papier [35], nous rappelons la généralisation de ces modèles aux cas avec une 
phase métastable ( Qf N Qe 4 0). Ensuite, nous introduisons la notion de solu- 
tion admissible pour les problèmes de Riemann et les solveurs de Riemann RSR 
and RSs. Le chapitre se termine avec des propositions sur la consistance et la 
continuité Li. pour les solveurs de Riemann RSR and RSs. 

Dans le troisième chapitre, nous décrivons le modele LWR contenant une con- 
trainte ponctuelle locale sur le flux. Plus précisément, nous considérons une situa- 
tion dans laquelle le flux maximal de voitures est limité par un point qui reste fixe 
sur la route. Grace à ces considérations, nous pouvons décrire le trafic à des postes 
de péage et à des chantiers. Nous définisson la solution entropique du probléme de 
Cauchy et nous rappelons le résultat d'existence correspondant. 

Le quatrieme chapitre est dédié au modele ARZ avec une contrainte ponctuelle 
locale sur le flux et à nos résultats obtenus dans [42]. Dans ces travaux, nous 
montrons l'existence de la solution faible correspondant au solver de Riemann non 
conservatif dans la classe des fonctions à variationsbornées. Le but est obtenu 
en montrant la convergence d'une suite de solutions approchées construites à par- 
tir de la méthode de Wave Front Tracking. Plus précisément nous introduisons 
un maillage, un solveur de Riemann approché CR S5; en coupant une onde de 
détente et nous construisons des problemes de Cauchy approchés. Grace de la 
fonctionnelle Y décroissante en temps, nous montrons que la variation totale des 
solutions accrochées est bornée uniformément. Gráce au théoréme de Helly, nous 
obtenons la convergence la solution approchée et nous montrons que la limite est 
bien la solution faible du probléme de Cauchy du modele ARZ avec une contrainte 
ponctuelle locale sur le flux. 

Dans le chapitre 5, nous décrions les modèles PT* et PT? avec contrainte 
ponctuelle locale sur le flux et nous présentons les résultats obtenus dans [10, 35]. 
Plus précisément nous introduisons les solveurs de Riemann CRóg et CRS3, avec 
des phases métastable et non métastable. Ensuite nous étudions leurs consis- 
tances, leurs continuités LL. et leurs domaines invariants. Le reste du chapitre 
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est dédié à l'existence d'une solution faible dans la classe des fonctions à varia- 
tions bornóes pour le modele PT? avec une phase métastable. Le but est obtenu 
en montrant la convergence d'une suite de solutions approchées construites par 
la méthode du Wave Front Tracking. De maniére similaire aux résultats obtenus 
au chapitre précédent, nous définissons un maillage et un solveur de Riemann ap- 
proché CRS". Ensuite, nous introduisons la fonction T décroissante en temps 
et nous montrons que la solution approchée est à variations bornées, le nombre 
d'ondes et d'interactions est fini en temps fini. Nous appliquons le théoréme de 
Helly et nous montrons que la limite est une solution entropique du probléme de 
Cauchy PT” avec phase métastable. 

Le sixiéme chapitre est dédié aux résultats obtenus dans des actes de conférences 
[34,43]. Nous y considérons deux modèles macroscopiques sur réseau. Le premier 
est le modele LWR avec des contraintes mobiles sur le flux. Cette idée de con- 
traintes mobiles sur le flux nous permet de modéliser la situation dans laquelle un 
camion (ou tout autre véhicule lent) réduit le flux à sa position. D'un point de vue 
mathématique, la contrainte est donnée par une équation différentielle ordinaire 
dépendant de la trajectoire du camion. Nous donnons une description détaillée 
du modele pour une route unidimensionnelle, introduisant un solveur de Riemann 
BRS wr et le généralisant à un réseau de routes. Le dernier modèle considéré est 
le modéle PT introduit dans le second chapitre. Nous le généralisons dans le cas 
d'un réseau de routes en introduisant un solveur de Riemann approprié. 

Enfin, dans un soucis de clarté et pour une plus facile compréhension, toutes 
les preuves techniques ont été déplacées dans l'annexe. 
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Vehicular traffic modelling 


1.1 Introduction 


Traffic congestion is one of the most challenging problems for nowadays cities. 
Transportation problems are studied by scientists from different research fields, 
such as Mathematics, Physics, Engineering, Psychology and Sociology. Currently, 
owning a vehicle is not only a convenience, but often is a must. Indeed, the public 
transport is typically inappropriate and forces the daily use of individual cars. As 
a result, we observe a steady escalation in road congestion over the past decades. 
This deeply affects environmental pollution, stress level of the drivers and fuel 
consumption. 

Traffic management systems based on mathematical models may contribute to 
solve some of the mentioned problems. Automatic control of traffic, thanks to 
real-time simulations, can improve the flow of the cars and optimize the routes of 
individual drivers. Therefore, building new roads to minimize traffic congestion, 
which is often not affordable, might be not necessary. However, building new 
roads not always diminishes roads congestion. As Braess’ paradox shows [17], 
expansion of a road network redistribute traffic flow and by that may cause its 
general decline. This phenomenon was observed among others in Stuttgart (1960), 
where closing the section of Kónigsstrafe street improved traffic flow in the area of 
Schlossplatz [61]. In general, such phenomena are not easily observable. Instead of 
experiment empirically on working roads to show particular behavior of the traffic, 


traffic engineers can use mathematical models. 
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For decades several approaches to traffic flow modelling were proposed to solve 
key transportation problems. The main author's interest during doctoral studies 
was macroscopic traffic flow modelling. The history of macroscopic modelling 
dates back to 1950s and arose from fluid dynamics concepts. The first macroscopic 
traffic flow model, proposed by Lighthill and Witham [68], is based on similarities 
between traffic flow and flood movements in rivers. Despite age and variety of 
approaches, this branch of science is still actively developing, raising new challenges 
and proposing new mathematical tools. 

The following thesis contains the state of art on selected issues related to macro- 


scopic traffic flow modelling and own contribution to the topic. 


1.2 Mathematical Models 


In this section, we briefly recall the main classification of mathematical models with 
a special focus on differences at the level of details. We show disparity between 
fluid and traffic flow. At last, we list the minimal requirements for physically 
reasonable traffic flow model. 


According to [82], mathematical models can be classified based on the following: 


e Level of details (microscopic, mesoscopic, macroscopic). 


Scale of the independent variables (continuous, discrete, semi-discrete). 


Representation of the processes (deterministic, stochastic). 


Operationalisation (analytical, simulation). 


Scale of application (networks, stretches, links, intersections). 


Let us have a deeper look into the classification with respect to the level of details. 

Microscopic models describe each vehicle independently, according to its 
speed and headings. It is often assumed that driver's decisions are made according 
to the traffic situation in front of his car. However, more complex models might 
take into account additional factors like interactions with other vehicles or driving 
style. The trajectory of each vehicle is given by an ordinary differential equation, 


whose solution gives the vehicle's trajectory. Due to the model's complexity and 
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the fact that the size of the system grows proportionally to the number of vehicles, 
online simulations are computationally very demanding. Moreover, solutions to 
system equations depend on each other, which makes parallel computation difficult 
to apply. The literature on microscopic models is broad and contains numerous 
number of models, like self-distance models [74,79], stimulus-response models [22, 
51] or cellular-automata models [57, 73]. 

Mesoscopic models fill the gap between microscopic and macroscopic models. 
The traffic flow is described in aggregate terms like in a probability distributions, 
however behavioral rules are considered for individual drivers. We distinguish 
three main kinds of mesoscopic models: headway distribution models [19,67], clus- 
ter models [63] and gas-kinetic models [76,80]. Mesoscopic models give a good 
compromise between the level of details and size of the area that can be analysed. 

Macroscopic models describe traffic flow by variables averaged over multiple 
vehicles: density, velocity and flow. The solution can be often expressed in a closed 
analytical form. Online simulations are capable, relatively easy to optimize and 
calibrate. 

Macroscopic modelling of traffic flow has roots in fluid dynamics. Though, it 
is intuitively clear that the dynamics of cars on a highway is different from that of 


flowing particles. We list below the main differences. 


e The driving behaviour differs among people. It depends on their experience, 
preferred style of driving but also on mood and tiredness at the particu- 
lar moment. Some behaviours can be also culturally conditioned and vary 


between countries. The fluid particles always follow physical laws. 
e The number of cars is far smaller than a number of particles in a fluid. 


e Fluid is isotropic, namely it has no directional preference and responds to 
stimuli from surrounding particles. Traffic flow is anisotropic, namely drivers 


move in one direction and they are influenced only by cars in front of them. 


Differences between traffic flow and fluid flow require a different approach to mod- 
elling. Aw and Rascle in [7] proposed the following minimal requirement for phys- 


ically reasonable macroscopic traffic flow model. 


(A.1) The system must be hyperbolic. 
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(A.2) The solution to arbitrary bounded nonnegative Riemann data in a suitable 


region of the plane must remain non-negative and bounded from above. 


(A.3) The speed of propagation of solution to any Riemann solver must be at 


most equal to the average speed. 


(A.4) The shock wave, whose propagation speed can be either negative or non- 
negative, are created due to breaking, whereas rarefaction waves are pro- 


duced due to acceleration and satisfy (A.3). 


(A.5) Near the vacuum, the solution to the Riemann problem must be very sen- 


sitive to the data. 


The condition (A.1) guarantees conservation of the number of cars, namely vehi- 
cles are not created nor destroyed. The condition (A.2) is a minimal requirement 
for any reasonable model. The condition (A.3) states that car travelling with a 
certain velocity receives no information from the rear. The condition (A.4) is nec- 
essary due to real life observation of traffic. The last condition (A.5) states that 


there is no continuous dependence with respect to initial data nearby vacuum. 


1.3 The macroscopic traffic variables 


We already briefly described different approaches to traffic flow modelling. Nev- 
ertheless, none of the models can be applied in a wide range if needed data are 
hard to collect. What kind of information are therefore useful and easy to reach? 
History of handling data of traffic flow dates back to Grenshield's studies in the 
1930s, what involved photographing the road section on fixed time intervals. 


According to [60] traffic can be observed from three perspectives: 


e Local(fixed position): camera, loop detector or other devices which captures 


changes at a certain point on the road. 


e Instantaneous(fixed time): camera or another device which capture longer 


road section at a certain time(e.g pictures from a helicopter). 


e Trajectory(moving with vehicle): in-car devices allowing to determine the 


position of the vehicle. 
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Which variables can be obtained using or mixing the above techniques? Let us 
start from studying the case in the microscopic framework. 
Consider a section of a multilane, unidirectional road without on-ramps and 


off-ramps. The position of the i-th vehicle at time t > 0 can be described by a 


unique pair X;(t) = (z;(t), (;(t)) € R x N, where the first component gives the 


position on the road section and the latter, number of the lane. Therefore we can 
give a graphical description of the traffic by providing the space-time diagram. 

For example, in Figure 1.1 the position of the vehicles is measured along the 
horizontal axis and time is represented on the vertical axis. The curved arrows 
correspond to trajectories of each car. We may observe that the driver starting 
from «5 is rapidly accelerating and overtaking three cars while the rest are rather 
timid. Two cars(starting from zy and xg) are moving along straight lines, namely 
their speed is constant and probably one following the other. The velocity of each 
car is given by v = z;(t). 

However, keeping track of every car might be problematic and very demanding 
for transportation systems. Instead of exact trajectories, we can study the velocity 
field v = v(t, x, £). Since for a given time and position on the specified lane, only 
one car may exists, the velocity vector field is uniquely determined. We observe 
also that the velocity field satisfies the equation v(t, X;(t)) = z;(t). 


da 


UZUAU5 XE L7L8 XQ X10 V11 


Figure 1.1: The example of vehicles trajectories on a road. 


In the macroscopic framework the velocity v is averaged over multiple vehicles. 
We can easily find two other macroscopic variables. The first one is called traffic 
flow, that we will denote by f, and describes the number of cars passing through 
a point in a time unit. The second is called traffic density, denoted by p, and 


describes the number of cars observed at a fixed time per unit space. 
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Example 1.1. Consider one lane, unidirectional road section of the length 1, 


parametrized by a coordinate x € R, with vehicles moving im the direction of 


increasing x. Let us study the simplest situation, namely with identical cars of 
length L, moving with the same constant speed v such that the distance between 
two consecutive vehicles is constant and equals d, see Figure 1.2. 

The traffic density, namely the number of vehicles per unit space, is constant 


in this case and expressed by 


pU 


The maximal density is achieved by taking d = 0, which is related in reality to the 


traffic jam with cars standing bumper to bumper, and therefore 


1 


Pmax — 


5l 


The traffic flow illustrates the number of vehicles per unit time. Let place an 
observer, who is counting vehicles in front of him. Each vehicle passes him every 


t=(L+d)/v, hence the traffic flow computed in a time unit is represented as 


U 


l=rzd 


Again, the above expression takes its biggest value for d = 0. Such a situation 
is however unrealistic if the cars are in motion. The real traffic flow observation 
shows that d — 0 is related rather to cars stuck in a traffic jam with maximal traffic 


density Pmax- 


observer 


&» S» €» €» & €» 


z 
L d vt 


Figure 1.2: Vehicles on a road section moving with the same speed 
v, having the same lenght Z and headings d. 
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1.4 The fundamental equations 


In this section, we give a mathematical description of macroscopic variables, intro- 
duce the fundamental equation and present the derivation of a scalar conservation 
law. Example 1.1 gave us a good intuition on traffic density and traffic flow con- 
cepts. One could also observe the simple relation between macroscopic variables, 


namely 


Jia) = plt, cyu(t, £). (1.1) 


Nonetheless, our example is oversimplified and we may wonder wheter it holds 
true in general. Consider the number of cars passing through the point x = x in 
a small time interval [to, to + At]. If we assume density and velocity functions to 
be continuous in time and space, then they might be approximated by constant 
values. Therefore, the number of cars passing through the point x, in a time At is 
approximately equal to v(t, x)p(t, x)At and (1.1) is still valid. This construction 
can be done only on a long road section in a short time interval or analogously 
on the short road section and longer time interval. However, it is not sufficient to 
derive the evolution of traffic in real life. 

Does formula (1.1) holds true for a more general framework? To answer this 
question, we need to define macroscopic variables at every point (x,t). Therefore 
let us assume the number of cars N (x,t) to be continuous and continuously differ- 
entiable. To find the local instantenous density we take the road section |x, «+ Ax] 


at given time t and let Az + 0, namely 


N(t,z + Az) — N(t, x) 
A0 Ar 


= ðN (x,t). 


Similarily, we may introduce a local instantenous traffic flow. Consider time inter- 


val [t,t + At] at given space x and let At + 0, namely 


Fa) = lim NEA?) - Na) 


AU At = ON (x, t). 


Velocity is a rate of change of vehicles position with respect to a frame of reference, 
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and is a function of time. By the above definitions we easily observe that for p 4 0 


f(t, x) 
p(t, z) 


v(t, x) = 


Thanks to the above considerations we can study the simple situation of traffic 
evolution on a road section [x1, £2] in a time interval |t, tg]. For the road without 


entrances and exits, the following formula conserving the number of cars holds true 


wo T2 t2 t2 
| p(x, tz) œ= f pa,t)da + | p(x1, t)v(z, t) a- | p (x2, t)v(z», t) dt. 

zı T1 ty ty 
(1.2) 


Assuming the functions p(x, t) and v(x, t) to be differentiable we may write 


t2 
He) = f Ojp(a,t) dt 
t 


1 


and 
pa, t)u(xe, t) — p(z1, t)v(z1,t) = JĄ On (p(x, tju(x,t)) dz. 


The equation (1.2) can be then written as 


| 2 ję {0;p(x, t) + O, (p(z,t)v(z,t))) dx dt = 0. 


Since a choice of a road section [x1,x2] and a time interval [t;, to] is arbitrary, we 


conclude that indegrand must be identically equal zero, namely 
Op + O,(pv) = 0. (1.3) 


This equation is a scalar conservation law in one space dimension. From the 
modelling point of view, it represents the conservation of the number of cars, 
namely vehicles are neither created nor destroyed. 

Since we have one equation with two variables, additional informations are 
needed. In this sense traffic flow models based on conservation laws could be 


divided into two main groups: 
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e Equilibrium traffic models assume that the traffic velocity in (1.3) is a 
function of density, namely v = v(p). Moreover, the dynamics of the traffic 
flow occurs along the equilibrium curve ((p, v(p)): p € [0, Pmax|}. They are 


also called first order models. 


e Non-equilibrium traffic models add additional partial differential equa- 


tions closing (1.3). They are also called higher order models. 


The first equilibrium traffic flow model was proposed by Lighthill and Witham 
[68] and, independently, Richards [81] (LWR). The traffic dynamics is expressed 


by the scalar conservation law 


Op + O,f = 0, f = plp). 


The function v : [0, pmax] — [0, Vmax] is such that v(0) = Vmax and v(Pmax) = 0. 
It is reasonable to assume also that v is C! and non-increasing. The LWR model 
will be presented in details in Section 2.2. 

In real life it has been observed that traffic is generally in non-equilibrium state. 
For small densities equilibrium models perform fairly good, however they do not 
describe correctly the dynamics of congested roads. This naturally leads to non- 
equilibrium models, which thanks to additional partial differential equation allow 
to consider non-equilibrium states. 

The first non-equilibrium traffic flow model was introduced by Payne and 
Witham [77]. They closed the scalar conservation law (1.3) by the equation ex- 
pressing traffic acceleretion 


z 2 
vU + vO,v = Ua ZY > X (1.4) 
T p 


where v(p) is an equilibrium velocity, c? is a diffusion parameter and 7 is a relax- 
ation time. In the 1995 Daganzo in [32] pointed out drawbacks of higher order 
models based on fluid dynamics. For example, vehicles in Payne-Whitham model 
can move backwards at upstream jam fronts, which turns out the model to be 
non-realistic. 

In 2000 Aw and Rascle in [7] and independently Zhang [85] proposed a new 


second-order model fixing flaws underlined by Daganzo. The Aw-Rascle-Zhang 
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model (ARZ) consists of the conservation of mass equation (1.2) and the equa- 
tion that mimics the momentum equation. More precisely, instead of unrealistic 
space derivative of pressure from Payne-Witham model, ARZ uses the convective 
derivative O; + v0,. By the fact that there is no conservation of momentum in the 
car traffic, the pressure term has been substituted by an "anticipation factor" p, 
which describes drivers” reaction to a variation of cars respect to the space. We 


may write the system in conservative form away from the vacuum as 


with 


Y = (ou) € R (0), ri)» (1-79) 
where the quantity y is a generalized momentum and p(p) is an anticipation factor. 
The ARZ model will be presented in details in Section 2.3. 

Nevertheless, none of the discussed construction is free of drawbacks. One of 
the main issue raised about the LWR model is infinite acceleration. In particular, 
when vehicles leave a congested road and enter an empty road, they immediately 
achieve their free flow speed. This can be solved by imposing additional term 
responsible for bounding acceleration, see [65]. Moreover, empirical studies show 
that the dynamics of traffic flow should be given on (p, f)-plane by a cloud of 
points rather than the equilibrium curve. 

The ARZ model fails to show continuous dependency of a solution with respect 
to initial datum near the vacuum. Furthermore, the maximal speed of the vehicles 
on an empty road depends on an initial datum. The full list of drawbacks is post- 
poned to another section, but the one mentioned here allow us to give a motivation 
for models with the phase transition. 

The idea behind phase transition models is to solve some drawbacks related 
both to LWR and ARZ models. The phase transition model treats differently traffic 
with low and high densities. According to Greenshield studies [55], the relation 
between traffic velocity and traffic density is (almost)linear. For this reason, the 


LWR model performs well in free flow regime, namely in the region with low density 


Nikodem Dymski 


Section 1.4 Page 11 


and high velocity. From experimental observations it appears that traffic in the 
congested regimes, namely sufficiently big densities, the fundamental diagram is 
2-dimensional and should contains a cloud of points. Therefore, in the congested 
region, the traffic flow can be modelled by a 2 x 2 system of conservation laws. 


Phase transition models can be written as follows 


Free flow Congested flow 
Y = (p,y) € Qr, Y = (p,y) e. 
(Y) = po(Y) 
Op + ôx f(Y) = 0, Op + sf (Y) — 0, 
v(Y) =v(p), dy + Or (y v(Y)) — 0, 


Above, p € [0, R] represents the density and y > 0 the generalized momentum 
of the vehicles. The sets Qs and Q. denote the invariant domains of the free and 
congested phases, respectively. Observe that in Qr the density p is the unique 
independent variable, so it is 1D while in (2, the independent variables are both p 
and y. 

The models briefly introduced above can predict and control with success the 
traffic flow on a simple road sections. However, real life traffic flow is much more 
complicated, thus more tools are needed to apply traffic flow models in a broad 
sense. For instance let us consider a road with traffic lights, toll gates or construc- 
tion sites. From the modelling point of view, all of them are "obstacles", reducing 


the traffic flow at fixed points on a road. We write such a constraint condition as 


p(t, $4) Vv (p(t, x;)) < qi(t), 


where x; are the "obstacles" positions, while q;(t) are the maximal flows allowed 
through them at time t 2 0. The concept of macroscopic traffic flow models with 
flux constraints has been popular in recent years. However, the first time this 
idea was proposed in a crowd dynamics framework to model the evacuation of a 
corridor through the exit door [30]. In the framework of traffic flow modelling, the 
problem was studied first in [29] and have been developed in [5, 20, 25, 26, 37, 38]. 
The idea of local point constraints can be generalized. Consider, instead of 


fixed obstacle, a slow vehicle moving like a bus or a truck reducing the traffic 
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flow at its position. This type of constraint is called a "moving bottleneck". The 
traffic evolution can be described by a strongly coupled PDE-ODE system. The 
coupling of conservation laws with ODE has been widely studied in [16, 37-39, 64]. 
More precisely, a system of PDE describes the traffic evolution, while an ODE 
describes the trajectory of the slow vehicle. In the PDE part, we might consider 
the already mentioned models like LWR, ARZ or PT. The constraint condition is 
constructed by X = z — y(t) change of coordinates, where y(t) is the slow vehicle 
position at time t > 0. In (X,t) coordinates, the velocity of the bus equals zero 


and conservation of cars equation (1.3) becomes 


Op + Ox (f(p) — Up) = 0. 


Therefore, the corresponding constraint condition might be written as 


f(p) — 9e € F(y(t)). 


The following PhD thesis contains the main results obtained during author's 
doctoral studies. His main interests in this period concern macroscopic modelling 
of traffic flow with constraints. The outcome of collaborative work with supervisors 
and other collaborators are three journal articles and two conference proceedings. 
The main ideas of such papers are stated below. 

In [42] we considered the ARZ model with fixed point constraint on the flow. 
We recall that the authors in [46] introduced two corresponding Riemann solvers, 
one fully conservative and one non-conservative. In our work we prove the existence 
of the weak solution, corresponding to the non-conservative Riemann solver, in the 
class of functions with bounded variation. The goal is obtained by showing the 
convergence of a sequence of approximate solution constructed via the Wave Front 
Tracking method [18,59]. A case study to describe the qualitative features of the 
solutions is also presented. 

The article [35] deals with a phase transition model with fixed local point 
constraint on the flow. We generalize the two PT models considered in [11, 12, 
53] and [14,15]. For more clarity, we consider two different PT models, both 
with metastable phases Q; N Qe # 0 and non-metastable phases Q; N Q. = (). 


The main result consists in the definitions of two new Riemann solvers and the 
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study of their properties. More precisely, we study the total variation estimates 


1 


ine Continuity. A case study 


and consistency of all Riemann solvers and their L 
describing the qualitative features of the solutions is also presented. 

In [10], we consider the constrained phase transition model with metastable 
phase. We prove the existence of a weak solutions in the class of function with 
bounded variation. The result is obtained via the Wave Front Tracking method. 
We point out that the main theorem distinguish two cases, namely when the con- 
straint q is higher or lower than the flux related to minimal density in the con- 
gested phase. However, the latter case requires supplementary conditions. It is 
worth mentioning that the solutions satisfy the entropy inequality with entropy 
pairs introduced in [11]. 

The two conference proceedings [34,43] concern traffic flow models on networks. 
The first is a generalization of the paper [35] to the case of the junction. The latter 
focuses on the LWR model with moving bottleneck. Both papers have a similar 
structure, namely, they define admissible solutions to the Riemann problem at the 
junction and introduce Riemann solvers generalized to the case of traffic networks. 


At last, we present a case study to give an intuition of solving the simple problems. 
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Macroscopic models 


2.1 [Introduction 


This chapter is devoted to macroscopic traffic flow models. Their purpose is to de- 
scribe the dynamics of traffic flow by using variables aggregated over road sections. 
We distinguish three main variables, namely (mean) density, (mean) velocity and 
(mean) flow, whose derivation is shown in Section 1.4. Macroscopic models are 
relatively simple and by that allow to real-time simulations of large traffic volumes 
and traffic on road networks. On the other hand, they are complex enough to cap- 
ture traffic flow phenomena. Moreover, macroscopic models are able to reproduce 
important features like formation and dissipation of traffic queues or appearing of 
shocks. With that being said, macroscopic models are very applicable to control- 
ling, predicting and optimizing traffic flow. 

The available literature extensively describes the different approaches to macro- 
scopic modelling. The traffic flow community actively develops new tools, both an- 
alytical and numerical, to solve problems motivated by real-life applications. In this 
section, we focus on issues related to the author's work during his doctoral studies. 
For information on different approaches we refer to state-of-art reports [9, 72, 78,84] 
and books [48, 60,82]. 
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2.2 LWR model 


The simplest and most widespread first order macroscopic model was derived by 
Lighthill and Witham [68] and independently Richards [81]. The dynamics of 
traffic is represented by a conservation of the number of cars equation, namely 
scalar conservation law 


ap + ðf (P) = 0, (2.1) 


where p = p(t, z) is the (mean) density at time t > 0 and position z € R, f = pv 


is the flux and v = v(p) is the (mean) velocity. Equation (2.1) represents the fact 
that vehicles are neither created nor destroyed on a the road without entrances 
and exits. We recall that roads with entrances or exists could be considered by 
adding a source term, see for instance [8] and the references therein. 

The equation is completed by a velocity-density relation v = v(p). It implies 
that a traffic moves along an equilibrium curve {(p,v(p)): p € [0, R]}. This as- 
sumption does not match to the reality. It has been shown that for congested 
roads traffic flow data is more likely represented on (p, f)-plane by a cloud of 
points rather than a simple curve. Furthermore, it implies that a small change of 
the density causes an instantaneous change of velocity. Non-equilibrium models 
overcome this drawback by considering velocity as a fundamental variable, see for 
instance Section 2.3. 

Whatsoever, in equilibrium models, only one averaged car population is consid- 
ered. In reality, the roads are occupied, among others, both by fast (almost racing) 
cars and slow ones dedicated to city driving. For this reason the assumption that 
they share, for instance, the same maximal velocity seems to be disputable. For 
generalization to multi-class we refer the reader to [13]. 

The last issue we mention about assumptions of LWR model is that the over- 
taking is not allowed. A generalization to multi-lane case can be found in [24]. 

In the theory of traffic low modelling it is convenient to consider the graph of 
the flux function f(p) = pv(p), which is called fundamental diagram. Throughout 


the thesis we assume that the velocity and flux functions satisfy: 


e |0, R] 2 pr v(p) € [0, Vmax] is a non-increasing Lipshitz function such 


that — = m) 
at v( R) = 0 and e(0) = Vmax- 
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e (0, R] > p> f(p) = pv(p) € [0, fmax] is a bell-shaped function, namely 
there exists Perit € (0, R) such that f'(p)(Pp—pet) 2 0 for a.e. p € [0, R]. 


(F) 


The flux function f is not necessarily continuous; for models with discontinuous 
flux function we refer the reader to [41,71]. By assumptions (V), (F) there exists 


maximal flow fmax achieved at some "critical" density perit € (0, R), that is 


rase = Ff (Perit) = Perit U(Perit); 


and moreover we have f(0) = f(R) = 0. The maximal flow fmax is also known 
as a capacity flow. The critical density perit divides the fundamental diagram into 
two particular regions - free flow regime for p € [0, perit] and congested flow regime 
for p € (Perit, R], see Figure 2.1. Observe that it is possible to have the same flow 


for different regimes and therefore for different velocities. 


f 
d max 


Perit R p 
Figure 2.1: Fundamental diagram for the LWR model. 
The expression of the velocity function v is chosen in order to fit best experi- 


mental data coming from the road under consideration. In the literature one may 


find the following expressions for the velocity functions satisfying (V), (F): 


Greenshields [54] : — v(p) = vua. — = p, (2.2) 
Vmax if p = [0, Perit); 
Daganzo [33] : v(p)= critVmax / R . 2.3 
ue Pho Mc Meme (3-1) spe (pai OP 
Perit p 
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Vmax — Ve 


Vmax — EE Ed if p € [0, Poit); 
Smulders [83] : v(p)= " > (2.4) 

Perit max (5 = 1) if p = [Perit, R], 

R hid Perit p 


for some Perit € (0, R). Observe that the Greenshields' velocity (2.2) is expressed by 
a single linear, strictly decreasing function. Daganzo's velocity (2.3) is represented 
by a piecewise function such that the corresponding flux f(p) — pv(p) is piecewise 
linear, with increasing part for low densities [0, pci) and decreasing part for high 
densities [pcrit, R]. The Smulders velocity function (2.4) is given by a decreasing 
linear function for low densities, while for high densities it has the same expression 


as Daganzo's velocity. 


2.2.1 The Riemann solver RSrwR 


In the theory of non-linear conservation laws has been shown, for example by the 
method of characteristics [18,31], that even for smooth initial datum the solu- 
tion may develop discontinuities(such as shocks) in finite time. Therefore it is 
convenient to consider problems with discontinuous initial data. 

Let us study Riemann problems for the LWR model (2.1), namely the simplest 
Cauchy problem with Heavyside-like initial datum 


(2.5) 


Pe if £ < 0, 
OP + O, zx. NARI 
tP Fo) p(0, x) n fr > 0. 


Below we give the detailed definition of the classical Riemann solver 


RS wr: [0, R]? + C° (R4; BV (R; (0, R])) 


corresponding to Riemann problem (2.5) for the case of C? strictly concave flux 
functions and for the case corresponding to the Daganzo's velocity (2.3). We defer 
the reader to [1,6, 18,31] for informations about more general flux functions. We 
recall that RSpwR associates to any initial datum py,p, € [0, R] a self-similar 


solution p(t, z) = RStwrlpe, PRI(T/t). 


Definition 2.1. Assume (V), (F) and that f is a C? strictly concave function. 
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Then the Riemann solver RSywn: [0, R]? + C° (R+; BV (R; (0, R])) is defined as 


follows: 


(R.1) If pe € pr, then RSrwnlpe, pr] is the shock wave 


pe if € < o(pe, Pr), 


| (2.6) 
pr if € > o(pe, Pr), 


RS wrlpe, Prl(€) = 


where the speed o(pe, pr) of propagation of the discontinuity satisfies the Rankine- 


Hugoniot condition 
(pe; Pr)(Pe — Pr) = Flpe) — F(Pr)- (2.7) 


(R.2) If pe 2 py, then RSiwelpe, pr] is the rarefaction wave 


pe if € < F'(pe), 
RSiwrlpe, Prl(ś) = 4 (FOE HH (po) € € < Flor), (2.8) 
Pr if E 2 f' (pr) - 


Observe that, from geometrical point of view, o(p, pr) is the slope of the 


straight line passing through (p, f(pe)) and (pr, f(p;.)). 


Definition 2.2. Assume (V), (F) and that v is the Daganzo's velocity (2.3). 
Then the Riemann solver RSiwr: [0, R]? + C° (R+; BV (R; [0, R])) is defined as 


follows: 


(R.1) If py, p. € [0, Perit], then RStwrlpec, Pr] is the contact discontinuity wave 


Pe if E < Vmax; 


pé if EÈ Vmax: > 


RSiwelpe, Prl(Ś) = i 


(R.2) If pe, pr € [Perit, R], then RSrwn|pu, pr] is the contact discontinuity wave 


pe if € ub PeritVmax 


R — Keri | 

RSnwR [Pe Pr] (£) = » pu. 
p WEZ Bie 

ou Perit 


(R.3) If pr < pau < pe, then RSrwn|pu, Pr] is the juxtaposition of two contact 
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discontinuity waves 


cri Vmax 


R — Perit” 
RS EE. „p PeritVmax 
LWR [Pe; Pr] (£) — $ Perit f CET < 3 < Vmax; 
R E Perit 


Pr if E > Vmax: 


(R.4) If pe < pais < pr then RStwrlpe, pr] is the shock wave given in (2.6), (2.7). 


Example 2.1. Consider a barrier on the grade crossing placed at x = 0, which is 
closed at initial time t = 0. For simplicity let f be the flux function corresponding to 
Greenshields velocity (2.2) with R = 1=Vmax, namely f(p) = p(1 — p). Initially 
the traffic in x < 0 has density p, € (0,1) and velocity v(pę) = 1 — pe. The 


corresponding initial condition takes the form 


Pe if x < 0, 
p(0, z) = | 
1 df 


The above choice for p, — 1 is to reproduce the effect in x « 0 of the closed barrier. 


Since pe < p, = 1, by (2.6) the solution is the shock wave 


KE n if x <o(pe,pr)t, (2.10) 


1 46 Sol pp pert, 


where o(pe, pr) < 0 is given in (2.7). 

The solution (2.10) represents a process of traffic jam formation with the dis- 
continuity line x = o(pe,p,)t separating still moving vehicles from those already 
caught in the traffic jam. 

A vehicle starting at the point x(0) = zo < 0 at time t = 0 has the following 
properties: 


position: x(t) = 


Xo + v(pe) t if t < ty, 
O(Pe; Pr) ty ift 2 ti, 


. ift «t4, 
velocity: x'( 
ift ti, 
0, 


acceleration: x” 


Nikodem Dymski 


Section 2.2 Page 20 


where by definition 


To 
O(Pe; Pr) = —pe, v(pe) = 1 — pe, t= = —Xo. 
v(pe) — o(pe, Pr) 


Example 2.2. Let us consider here the same setting of the previous example, 
however with the assumption that the barrier opens att =0. The traffic in x < 0 
has maximal density, while in x > 0 there are no cars, namely pe = 1 and p, = 0. 


In this case the initial condition takes the form 


1 eee 
p(0, x) = 
0 aa > 0. 


By (2.8) the solution is the rarefaction wave 


if x < —t, 
Ev > 
C=) if-t<a<t, 
ER 


p(t, z) = 


O N|= =| 


From the modelling point of view, a rarefaction wave describes traffic acceleration 
and simultaneously gradual decrease of traffic density. 

Again, consider the vehicle starting at t = 0 from xg < 0. Observe that this 
particular vehicle is stuck in the traffic jam until tı = xo/ f'(p;) = —1o. Then its 


trajectory t > z(t) is given by the solution of the Cauchy problem 


vie vp > (: 4 22 ; E) = z0. 


The properties of the vehicle are the following: 


ti (t) Xo ift < ti, 
position: z(t) = | 
t—2Vtit if tti, 
0 ift «t, 
velocity: x'(t) = " 
je: > ift> t, 
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0 ift < ty, 


acceleration: DS 4 . 


2.2.2 Drawbacks of the LWR model 


The LWR model performs fairly well in some basic road situations. However, in 


most of the real situations car traffic is complex and simple models might not be 


able to reproduce important effects. In the articles [32, 69, 75] were listed some 


drawbacks of LWR model, what we give below: 


Vehicles are supposed to reach new velocity immediately after the change of a 
road density. The lack of any delay related to drivers reaction implies an infinite 


acceleration. 


The LWR model (with passing allowed) does not recognize the distribution of 
desired velocities for light traffic but only the desired velocity of each vehicle. 
The velocity distribution across vehicles tends to spread a platoon linearly with 


time, while the variation within each vehicle with the square root of time. 


In contrast to reality, the transition from the free flow regime to the congested 
regime always occurs at the same density. This leads to the same outflow after 


breakdown. 


The velocity of the traffic depends only on the density what is empirically shown 


to be wrong. Moreover, only one velocity corresponds to a certain density. 


Equilibrium traffic flow models do not characterize the amplification of small 


disturbances in heavy traffic, called phantom jams [45,56]. 


Equilibrium models do not perform hysteresis phenomena, that is asymmetry 
between acceleration and deceleration behavior of driver-vehicle units. It is 


caused by retarded recovery of speed in deceleration-acceleration process. 


2.3 ARZ model 


To overcome drawbacks of LWR (2.1) and Payne-Whitham (1.3), (1.4) models, 
Aw and Rascle [7], and independently Zhang [85] proposed a new second order 
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model (ARZ). The model states the conservation of the total number cars and the 
conservation of the generalized momentum. 

In this section, we first present the ARZ model both in conservative coordinates 
Y and Riemann invariants coordinates W. The reason for providing W coordinates 
is due to indefiniteness of the system in conservation form at the vacuum. The 
Riemann invariant coordinates are also handy to calculate total variation estimates 
and provide easy geometrical interpretations for solutions to Riemann problems. 
Next, we give the corresponding Riemann solver in W-coordinates. At last, we 


define weak and entropy solutions. 
2.3.1 ARZ model in Riemann invariant coordinates 
'The ARZ model takes the form 


ıp + O,(pv) = 0, 


(2.11) 
dv + p(p)) + vd, (v + p(p)) = 0. 


Above p = p(t,x) and v = v(t, x) are density and velocity functions at time t > 0 


and position x € R, while p(p) represents the anticipation factor accounting for 


drivers! reaction to the state of traffic downstream. We assume that 
p(0)— 0 and p(p)>0, p(p)>0, 2p'(p)-pp'(p)» 0 for any p > 0. (2.12) 


Typical choice is p(p) = p", y > 0, see [7]. 

By assumption (A.1), reasonable traffic flow models have to be expressed by 
hyperbolic systems of conservation laws. We show now that (2.11) can be rewritten 
as a hyperbolic systems of conservation laws. By multiplying the first equation of 
(2.11) by v + p(p) and by adding the second equation in (2.11) multiplied by p we 
obtain 

Ory + Os (vy) = 0, 


where y = p(v + p(p)) is the generalized momentum. As a result, system (2.11) is 


equivalent to the system of conservation laws 
OY +0,F(Y)=0, (2.13) 
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with 


Y = (p,y)” € R2 \ {0}, rev) = no] 


The Jacobian matrix of the flux function F(Y) is given by 


—p(p) — ppp) 1 
A(Y) = DF(Y) = y? p) 2 (p) 
——, — YPP — — PP 
p? p 
The eigenvalues of A(Y) are 
AU m 
A(Y) = > — p(o) — PP (P), AY) = ; — p(p), 
with corresponding eigenvectors 
r(Y) = Y, r(Y) = (p, y + p?p (p))’- 


By the assumptions in (2.12) we have that Aı(Y) < A2(Y) for every p > 0, so 
system (2.13) is hyperbolic away from the vacuum. 


By direct computations we obtain 
VAY) r(Y) = —0(27 (0) + py'(p)) <0, VAY) r(Y) =0, 


namely the first characteristic field is genuinely non-linear, while the second char- 
acteristic field is linearly degenerate. 

We recall that the speed of propagation o = o(Y,,Y,) of any discontinuity 
(Ye, Y,) satisfies the Rankine-Hugoniot conditions 


(Y. -Y)e = F(Y.) — F(%). (2.14) 


In order to construct elementary waves, we use the Riemann invariant coordi- 
nates U = (v, w)”, which by definition [66] satisfies 


V yv : T3 = 0, Vyw-r, = 0. 
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Away from the vacuum p = 0 the Riemann invariant coordinates U are linked to 


the conservative coordinates Y by the identities 


We recall that w is a Lagrangian marker. The Riemann invariant coordinates U 
are a good choice because in this coordinates the total variation of the solution 
does not increase in time [44]. 

For any element U of the space U = {U € Ri:v < w) we can define the 
corresponding density p(U) = p "(w — v), generalized momentum y(U) = p "(w — 
v) w and flow f(U) = p '(w — v)v. The vacuum state p = 0 is described by a 
half-line Uy = {U € U : v = w} and non-vacuum state by US = UV. 


2.3.2 The Riemann solver RSARZ 


In this section we present the Riemann solver for ARZ model (2.11). Consider 


Riemann problem 


O,w + vO,w = 0, (2.15) 
U (0, z) = Uo(a), 


with Riemann initial datum 


U ifz«0, 
U, = 2.16 
olz) Ba eu) 


where U;, U, € U are constants. 

Before we introduce elementary waves, let consider the inverse function II € 
C? (R+; R,)C! ((0, +00); R+) of the function p > p(p)+pp'(p). By the assump- 
tions in (2.12) we have that II(0) — 0 and II is strictly increasing. The Riemann 
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solver RSARz for problem (2.15), (2.16) 


RJSARZ : Uxu > C? ((0, +00); Lig, (IR; 4)) , 
(Ur, U,) > RS arzlUe, Ur], 


can be described in terms of the following three elementary waves. 


e Take U; = (ve, we)? EU, U, = (vp, wr)? € uc, with w, = w, and v, < w. The 


elementary wave joining U; with U, is the shock wave 


U, if£ < o(U„U.), 


S[U,, U,]($) = i if £ > a (Uz, U,) 


where 


F(U.) Sz f (Ue) 
(Ur) — p(Ue) 


e Take U, € Uj, U, E U, with we = w, and v, < vr. The elementary wave joining 


o(U;, U,) = 


U, with U, is the rarefaction wave 


U; ie « X (U4), 
R[U+„U,J(E) = (" — p (II (w a PRONE ESN, 
U, if € > A(U,), 


where 
. |v=p(U)p'(p(U)) if U € ut, 
u(U) = AU IELU d 0 
m if U € M. 


e Take U;,U, € B. with w, Z w, and v, = v,. The elementary wave joining U; 


with U, is the contact discontiuity wave 


U, if E < Ver, 


qu o = fy if E> u, 


We stress that both shocks and contact discontinuities satisfy the Rankine- 


Hugoniot conditions (2.14). Indeed, any discontinuity (Už, Už) with speed of prop- 
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agation o* = o(U/,U;) satisfies the Rankine-Hugoniot conditions (2.14), which in 


the U-coordinates take the form 
p(U;)(o" — vr) = e(U7)(oc" —v;), p(U;)(o* —v,)w, = p(Uz)(o* — v;)w;. (2.17) 


In particular away from the vacuum, namely for Už, Už € Ug, by (2.17) either 
W; = W, or v; = o* = v;, which correspond to respectively shock and contact 


discontinuities waves. 


Definition 2.3. For any Up = (ve, we)? ,U, = (v,,w,)” € U with U £U, and 
(Uc, U,) EU xU, we define RS arzlUe, U,] as follows: 


1. IfU,€ n. we = w, and v, < ve, then RS Anz|U5 U,] = SIUS U,]. 
2. IfU, € uc, we = w, and v, < Up, then RS arz[Ue, U,] = RIU, U]. 
3. IfU„U, EU and ve = uy, then RS arz[Us, U,] = C[Us, U,]. 


4. If Ug Em and U, € us with wę = w,, then 


U; if € < Ur, 


RSarzlUe, U; (6) = Ę SEA 


5. If Ue, U, € UË and v, < v, < we, then 'RSanz|U;, U,| is the juxtaposition of 
S[U;, Um] and CIU, Ur], where Um = (vp, we)? € Us. 


6. If Up, U, € us and v, < v, < we, then RSARz(U,,U,] is the juxtaposition of 
R[U;, Um] and C[Um, Ur], where Um = (vp, we)? € UE. 


7. If Ue, U, € us and v, < we € v, < wp, then 


R[U;, Um (€) if E < Ur, 


RS arzlUe, U,| (6) = i if£ >v 


where Um = (we, we)? € My. 


8. If Up CUS, U, € Uo and we = w,, then RS arz[Us, U,] = R[Us, U,]. 
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9. IfU, € us. U, € Uo and we # wr, we can define RSarz[Us, U,] = RIU, Um], 


where Um = (wę, we)? € Ud. 
At last, we define RS jnz|U., U,] = U, for any U, EU. 


Remark 2.1. For any U, = (vs, Wx)? € US and U, = (Va, Vx)? € Uo then both 


Up usw 
U revu 


and 


RIU, Um|(x/t) à xt, 
| PIE) ra Um = (We, Ws)” € M, 


U, if £ > Wt, 


are weak solutions for the corresponding Riemann problem. However, as already 
observed by Aw and Rascle [7], only the first solution is physically reasonable. This 
justifies the choice of the item 9. in Definition 2.3. 


2.3.3 Weak and entropy solutions 
In this section we introduce the definition of solutions of the Cauchy problem 
p + O,(pv) = 0, 


Ow + v Ow = 0, (2.18) 
U(0, z) = Uo(z), 


where the initial datum Up is assumed to be in L^? (R;U/). Let Vo = ||Uol|ao. We 
extend the flux F to the whole of R, by taking F(0) =0. 


Definition 2.4. Let U, € L® (R; 4). We say that a function 


U € L*(R, x R;4) n C? (R4; Lig, (R:{)) 


is a weak solution of (2.18) if it satisfies the initial condition (2.18), for a.e. x € R 
and for any test function @ € Cz ((0, +20) x R; R) 


[[,.. 69 o2) G) dz dt = 0. 
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Proposition 2.1 (|4|). For any (U„U,)EUxU we have that 
U(t, x) = RSARZ(U:, U;] (x/t) 


is a weak solution of (2.18) with Riemann initial condition (2.16). 


Conservation laws admit in general more than one weak solution. This mo- 
tivates the introduction of an additional selection criterion, the so called entropy 
condition [66,70]. As suggested in [4] we consider entropy pairs (€, Q,) defined 
for any fixed k > 0 as 


0 ifv < k, 
Ex(U) = “ly — 
2 pg NY sep 
E D (2.19) 
0 ifv <S k, 
Q,(U) = f(U) 
k——————— if Ą 
AGE ifu>k 


Definition 2.5 (Entropy solution for the ARZ model). Fir Uo € LC (R;U). Let 
U € L” (R, x R;U)NC® (R4; LL, (R;U)) be a weak solution of (2.18) in the sense 


“loc 


of Definition 2.4. We say that U is an entropy solution if for any non-negative 


test function à € CX ((0, +20) x R; R) and for any entropy pair (2.19) we have 


If _ (EU) 26 + Q.(U) 0,0) dz dł > 0, k € (0, +00). 


2.4 PT models 


In this subsection we consider phase transition (PT) models of hyperbolic conser- 
vation laws for vehicular traffic. According to experimental data, the vehicular 
traffic flow acts differently depending on whether is free or congested. This leads 
to consider two regimes: free-flow phase Qf and congested phase Q.. The free-flow 
regime corresponds to low densities and can be approximated by a one-dimensional 
flux function, while the congested regime refers to high densities and the flow covers 
a two-dimensional domain, see [28, Figure 1.1] or [15, Figure 3.1]. Consequently, 


the traffic can be well described as a coupling of first order model (a scalar PDE) 
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in the free-flow regime and a second order model (a 2 x 2 system of PDEs) in the 
congested regime. 

This two-phase approach was first introduced by Colombo in [28] and after 
exploited by other authors in [11, 14, 15,53]. For instance, in [53] Goatin couples 
LWR model (2.1) for the free-flow phase Q; with ARZ model (2.13) for the con- 
gested phase Q.. Such model corrects drawbacks of LWR and ARZ models taken 
separately, see Subsection 2.2.2 and Remark 2.1. 

In the paper [35] we generalize two P'T models. The former one, denoted by 
P'T?, was studied in [11,53] and describes the coupling of the LWR model with the 
ARZ model. The latter, denoted by PT’, has been the subject of research in [14,15] 
where the authors considered Greenshields' type flux functions in congested regime. 
The congested regime is treated there as an extension of the free flow regime. In 
[14,15] the authors assume that the two phases have intersection, namely Or. Æ 
Ø, while in [11, 12,53] the authors assume that Q N Q. = Ø. In [35] we do not 
make any assumption on the intersection between free and congested phases. We 
suppose that the velocity in the free-flow regime is represented by a unique value, 
and coincides with the maximal velocity. This assumption is in order to avoid the 
loss of well-posedness, as observed in [28, Remark 2]. Therefore, in the free phase 
the velocity function coincides with Daganzo (2.3) velocity function. Moreover, 
in [14,15] the authors assume that the flux function vanishes at a maximal density, 
namely that the vehicles have (almost) the same length, while in |11, 12, 53] the 
authors do not impose this requirement. In [35] we consider both the cases, as 


both are motivated by practical reasons. 


2.4.1 The general PT models 


In this section we recall the PT models we introduced in [35]. Throughout the thesis 
we ignore the superscripts p, a only when they are not necessary. The common 
fundamental parameters for the PT? and the PT* models are: 

e R>0is the maximal value of the density; 

© Vmax > 0 is the maximal velocity; 


e v, € (0,Vmax] is the maximal velocity in the congested phase. 
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The PT model can be written as 


Free flow Congested flow 
Y = (py) € Oe, Y = (p,y) € Me, 
(p. y) € O (p, y) (2.20) 
Op + Onf(Y) = 0, Asp + 0, f(Y) — 0, 
v(Y je Vmax, Ory + Ox (y v(Y)) = (. 


Above, p € [0, R] is the density and y is the generalized momentum, while the 
domains €? and Qe refers to the free and congested flow phases, respectively. The 


corresponding velocity v(Y ) and the flow f(Y) are defined as follows: 


RAE v*(Y) = v (p) (+y) for PT*, us 
vP(Y) = ; — p(p) for PT”. á | 


The term (1 + y) in the definition of v^(Y) is a perturbation which corresponds 


to the thickness of fundamental diagram in the congested regime around the flux 


f(p,0) = pvz, (p). The function v2, : (0, R] > R+ is the equilibrium velocity for 


PT* model given by 


Above a € Rand r, € (0, R) are fixed parameters of the model. We point out that 


a 
eq 


has different physical meaning in the framework of PT^ and PT? models. 


v? (r) = Vmax by definition. At last, we stress that the generalized momentum y 


We recall the notation and the main assumptions on the parameters discussed 
in [11,15], see Figure 2.2. For PT? we require that p: (0, R] — R satisfies (2.12) 


p(0)=0 and p(p)>0, p(p)>0, w(p)+pp'(p)>0 for any p > 0. 


Fixy_ < y, and rf, r$ so that 


f 


(KAM < R, v(r3 ; ri y+/R) = Ve. 


0<r 


er < R, UU r4} y+/R) = Vmax, 
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pw.) R p pi) R p 


Figure 2.2: The first column refers to the case QN Qe = 0, namely 
Ve < Vmax, While the second column refers to the case Qr. z 0, 
ie. Ve = Umax. The first row refers to PT^ and the second row to 
zs Above e (rf, rfi Wh Yr = (tow als (X); 


= yf (Y) and Y? = v£(Y;) defined in Subsection 2.4.1. 


Notice that ri € r$, with the equality holding if and only if ve = Vmax. We can 


then explicitly characterize the free and congested regimes as 


Q = {Y € [ori] x Ry = QU). 
0, = [Y ege, R] x R:0<0(Y) «v, w- <2 <w. }, 
p 
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where w+ = y+/R and 


(o — ro) (a(R — p(R — ro) + Vmax) 
Q(p) = (=) (a (ro — p)(R — ro) + Vinaxto) 
P (Vmax + p(p)) for PT?. 


for PT 


Notice that v(Y) = Vmax for any Y € Qr and wa = Q(<)/r = Q(ri)/rí. More- 


over, we denote Q = 0; U (Y, and 


OF -2(YeO:pe(orf)) =Y EQ: pE frir] 
QD NOE Q = (Y € (0, R] x R :e(Y) € (0, vas], w(Y) € [w-,w4]}, 


where 


| |u/p if u € QF, 
w(Y) = p r - 
Wegen a ifu E Qr. 
r 


Under the conditions stated above, we point out that 


Voce => Mie, Dee, Ionen. 
waa S fe N=, (QU. 


c 


Remark 2.2. We underline that in the congested phase the variable w is a La- 
grangian marker, namely the equation O;w(Y )+v(Y ) Ozw(Y ) = 0 is satisfied, under 


the assumption that the weak solution Y to (2.20) takes values in Qe. 


We define functions that are useful to define the Riemann solvers given in the 


next subsections: 


R for RI”, 
b (w) for PTP, 
w(Y}) =w(%), 
Pi: (QUO) + (RUN), ed) = : 
v(Yi) — Vmax; 
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dy: (Q UN) + (RUN), =) = 


Y7: 0 (RUN), Y? = 03 (Yo) > i 


We stress that the maximal density is reached for o9(w,) = R both for the PT“ 
and PT? models. Since the considered PT models are given by 2 x 2 system of 
conservation laws in the congested phase, then they have two Lax curves corre- 
sponding to each characteristic field in Q., namely the curves along which the 
Riemann invariants are constant. We extend the Lax curves defined in Q. to Qg 


as the graphs of the maps 


[pi (w(Yo)), pt (w(Yo))] 3 0 9 Lycy(P) = f (o, w(Yo)p), 
[P> (v(Yo)) . ez (v(Y9))] 9 p = Lh lo) = v(Yo)p. 


for any Yo € Q7. Draw attention that the graphs of Lj, and £j, belong to the 
boundary of QS” in the (p, f)-plane. By using the Lax curves definition, we can 


give geometrical meaning of the functions introduced before, see Figure 2.2: 
e 1;(Y5) is the intersection of Lax curve Live) and (Y € 2: w(Y) = Vmax}. 


e i1(Yo) is the intersection of Lax curve Loy) and (Y € 2: v(Y) = ve}. 


e For any w € [w_,w,] the point (p?(w), p?(w) w) is the intersection of the Lax 
curve £l 


have p?(w) = p^!(w) for PT”, while p?(w) = R for PT". 


and (Y € Q. : v(Y) = 0}. More precisely, for any w € [w_,w+] we 


e y3 (Y) is the intersection of Lax curve LY, and (Y € 0: w(Y) = w4}. 


At last, we define the maps 


f be E and 3 : [0, Vmax] + [r5 p1(wa)] 


py : wo wy] 9 [rn 


so that py°(w(Yo)) and p£(v(Yo)) are the p-components of Y% (Yo) and v£ (Yo). 
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2.4.2 Main assumptions 


In this section we give eigenvalues, eigenvectors, Riemann invariants and main 


assumptions for the PT models in the congested phase. The eigenvalues in Q. are 
4 (Y) = v(Y) +Y :Vyv(Y), A (Y) =v(V), 


with corresponding eigenvectors 


PTY) = v, PTY) = Oy) ) | 


and Riemann invariants 
wi (Y) = w(Y), i Sy 


The above functions naturally extend to Q£”. We can see that APT(Y) > 0 for all 
Y € OC". We obtain by direct computation that 


VAPT(Y).rP'(y) 22Y - Vv(Y) + p? (VY) + 2 py0,0,0(Y) + y? Œv(Y), 
A Y): r$ (Y) = 0. 
As a consequence the second characteristic field is linearly degenerate. 

For modelling consistency, we assume that for any w € [w_,w4] the flux is 
decreasing and has at most one inflection point. For this reason, in Qe the waves 
of the first and second characteristic families have respectively non-positive and 
non-negative speed of propagations. Therefore, the following assumptions will be 


needed in order to construct general PT model: 
(H1) APT(Y) < 0 for all Y € NF 


(H2) the first characteristic field is genuinely non-linear in €£*, except for the 


PT? model. 


Notice that the first characteristic field for PT? model, namely PT* with a = 0, 


is linearly degenerate along the Lax curve £j and is genuinely non-linear along 
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any Lax curve £l, with w 4 0. It is interesting to observe that if Cl vanishes in 
Y € QS”, then (H1) and (H2) imply that v(Y) € (0, Vmax}. 


Remark 2.3. Regarding to the fact that for PT? model, the p function satis- 
fies 2p' (p) + pp'(p) > 0 for every p € [rr], (H1) reduces to the requirement 
r p'(r) > Vmax: Moreover, by (2.12) we have VA?T(Y) - rt? (Y) = —p (2p'(p) + 
py" (p)) < 0 for every p € (0, R] and (H2) easily follows. 

On the contrary, for P'T^ in general assumptions (H1) and (H2) cannot be 
easily derived in terms of the parameters of the model. We just mention that in 


the simplest case a = 0 condition (H1), as stated in [49], is guaranteed by 


-5 <w <0<w< p 
Remark 2.4. The assumptions (H1) and (H2) can be reformulated with regard 
to the first Lax curves. Since apap) = A?T(p, pw), then (H1) represents the fact 
that the first Lax curves are strictly decreasing. As a consequence the capacity drop 
in the passage from the free phase to the congested phase is ensured. Moreover, as 
a result of 35. £1, (0) = 5 VAT (p, pw): rPT (p, pw), we have that (H2) is equivalent 
to the fact that the first Lax curves are strictly concave or convez, except the case 
of P'T? corresponding to w — 0. Point out that by (2.12) we have ELL (0) <0 


for all p € [pl (w), p?(w)|, hence the first Lax curves for PT? are strictly concave. 


In the forthcoming subsections we recall the Riemann solvers corresponding 
to Riemann problems both with metastable phase Qf NA, Z 0, and without 
metastable phase, that is Qf N Qe = Ø. We introduce now the admissible solu- 
tions to the Riemann problem for both problems. Consider the Riemann problem 
for PT model (2.20) with the initial datum 


Y» HR, 
Y ae (2.21) 
Y, ifx >00. 


Let us recall the definition of solution to (2.20), (2.21) introduced in [28]. 


Definition 2.6. For any Y, Y, € Q a self-similar function Y = (p, y): Ry xR > 9 
is an admissible solution to (2.20), (2.21) if the following properties hold: 
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C.1 IfY„Y, € Qs, thenY=(p,y) with y = Q(p) and p=RSiwnlpe, pr] given by 


(2.9). The solution has no phase transition and attains values in Qs. 


C.2 IfY,,Y, € (X, then Y is the usual Lax solution to (2.20),, (2.21). The solution 


has no phase transition and attains values in Qe. 


C.3 If Y, € Q; and Y, € €, then there exists o € R such that: 


(a) Y (t, (7oo,ct)) € Qr and Y (t, (ct, +00)) € À for all t > 0. 


(b) The first Rankine-Hugoniot condition is satisfied for all t > 0, namely 
o (p(t, ot^) — p(t,ct-)) = f(Y(t.ot*)) — F(Y (t, ot7)). 
(c) The functions 


Y (t,x) Jd ae ot, Y(t,ott) <at, 
(t, c) > (t,£) > 
Noch dh Yes) iu > at, 


are respectively the usual Lax solutions to the Riemann problems 


op + Oxf (Y) =0, Op T osf (Y) =0, 
v(Y) — Vmax; dy RE O, (y v(Y)) = 0, 
Y, f x <0, Vote) ifx <0, 
Yoga u? Vibe d ER ae 
Vitor) 47 BU, Y, ia. 


C.4 In the case of Y; € (Y; and Y, € Q; analogous conditions to the previous case 


are required. 


Notice that for the PT? model, condition (C.2) states that Y = SAnz|Y; Y;], 
where RS Anz is given by Definition 2.3. 


2.4.3 The Riemann solvers RSR and RSs for PT models 


In this section we introduce Riemann solvers RSR and RSs for the Riemann 


problems with metastable phase and without metastable phase, respectively. We 
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point out that these Riemann solvers are defined below according to Definition 2.6, 
in the sense that (t,x) — RSR[Y, Y,|(1/t) and (t,x) > RSs|Y,,Y,|(z/t) are 
admissible solutions to the Riemann problem (2.20), (2.21). 

Besides the elementary 1-waves, namely rarefactions R and shocks S, and 2- 
waves, namely contact discontinuities C, which are defined analogously to those 
introduced in Subsection 2.3.2 for the ARZ model, we need to introduce phase 
transitions (PT). A PT wave is a shock-like wave between states Y; and Y,. belong- 
ing to different regimes. More precisely, we distinguish the following types of PT 


Waves: 
e YEO; \{(0,0)}, Y. = v5 (Y,) € Q \ Qs; 
e Y, = (0,0), Y, € Qe | Or; 
e Y, € OF, Y, € OF with f(Y.) = Cu (Pr) if Re N Qc = 0; 
e Y; = J£(Y;) € Qe, Y, = YEY) e OF if OLN 0; = 9. 
Therefore, the PT wave can be defined as 


Y, if € < o(Y5 Y»), 


PT[Y, Y,|(6) = E if £ > o(%, Y;), 


where 
FY) — fr) 
Pe—Pr 


Let us first consider a Riemann solver for the case with metastable phase, that 


(Ye Y.) = 


is Qf N Q. = Of 40. In this case the maximal velocities for both phases coincide, 


c 


namely Vmax = Ve. We recall that VÍ = WS, ri = r$ and for this reason we simply 


write i; and r+. Moreover, we write below f; for f(Y;), ve for v(Y;), we for w(Y;) 


and so on. 


Definition 2.7. The Riemann solver RSR: Q? — L®(R; Q) associated to Rie- 
mann problem (2.20), (2.21) is defined as follows. 


(R.1) If Y, Y, € Qg, then Y Z (p, y) with y = Q(p) and p = RSywn|pe, Pr] given 
by (2.9). 
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(R.2) If Yi, Y, € Qe then RSr|Y:, Y,] is the juxtaposition of possibly null 1-wave 
(Ye, Ym) and possibly null C[Y,,, Y.], where Wm = we and Um = Ur. 

(R.3) If Y; EQ; and Y, € Q; , then RSR[Y;, Y,] is the juxtaposition of a 1-wave 
(R.4) If Y, € OF, Y. EQ; and o (Ya, v5 (Y,)) > AT (u5 (Y,)), then RSn[Y,, Y;] is 
the juxtaposition of PT[Ye, va (Y,)] and possibly null Cya (Y,), Yr]. 

(R.5) If Y; € Og, Y, € Q and o(Yn (Y) < XP (5 (Y;)), then RSr[Yn Y. 
is the juxtaposition of PTY, Yp], R[Y;, v; (Y;)) and Clp; (Y), Y.], where 

Y, = Y; (Y?) is the state satisfying w, = w.. and o(Y4, Yp) = AFT (Yp). 


We point out that if = Ce) € 0, then o(%4,44 (Y,)) 2 APT (yz (Y,)) for all 
Yı € Q; and Y, € Q7; thus case (R.5) never occurs. We recall that by Remark 2.4 
the Lax curves of the first characteristic family for PT? model are always concave. 

Let consider now a Riemann solver for the case without metastable phase, that 
is Q N Qe = Ø. In this case the maximal velocities of each phase do not coincide, 


namely Ve < Vmax. We introduce the states 
Yr. = (r£, w- r£), += (r$, wy ri), 


which represents the points in 0, with minimal p-coordinate and maximal f- 


coordinate, respectively, see Figure 2.2. The states 
YE = (wii), Vi, = (wr), 


represents the points in QF with minimal p-coordinate and maximal f-coordinate, 


respectively, see Figure 2.2. At last we define 


po 


Definition 2.8. The Riemann solver RSs: 0? + L?*(R; Q) associated to (2.20), 
(2.21) is defined as follows. 
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(S.1) We let RSg|Y;, Y.] = RSnlY>, Y.] whenever 


2pl 


d°£ 
(Y, Y) ENR U OF U (v, Y) € Me x De: qr (PO > 0) 


u 1(%,Y,) EO; x0.: (Ve VE) > ATW) | 


diri 
U {M Y) € OF x Ne: gaz (Po) < 0). 
dA 
(S.2) If Y; € Qe, Y, € Qf and do (pe) < 0, then we let 


RSg[Ys v$(Y)])(£) for €< e(t (Yo), vi (Ye), 


RSs|Y,. Y,)(£) = 
RSa|vi(Y),Y.(£) for € > e(u£(Yo), yi (Yo). 


S.3) If Y, EQ, Y, EQ. and e(Y,, YE ) < APT (Y£ ), then we let 
f 


Y, for € < c (Yo, YE), 


RSs|Y,, Y.](£) = 
RSRIY_„Y.](6) for& 2o(Ys YE). 


21 


d^c 
(S.4) If Y; € OF, Y, € Qe and do (pe) > 0, then we let 


Y, for € < o(Ye, (Ve), 


RSs|Y, Y.(£) = 

Rón[v1(Y), VIE) for € 2 e(Yo vy (Ye). 
Remark 2.5. We point out that RSS differs from RSR only in the cases described 
in (S.2), (S.3) and (S.4), hence the Riemann solver RSs|Y;, Y,| differs from 
RSR|Y, Y,] if and only if (Ye, Y.) one of the following condition holds: 


dE 
Ye € 05 Y; E Or, a 00 < 0, (2.22) 
YEM, Y, E€ Qe, JOW YE) < APT (YE), (2.23) 
d? L1 
Ye OF, Y, EX, a (0 > 0. (2.24) 
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Observe also that for PT? model we have that RSg|Y+, Y,| differs from R.Sg|Y,, Y,] 
if and only if Y; € Qe and Y, € Qe. Moreover, this is also the case for PT* model 
2 pl 


if a (D) <0. 


At last we define a Riemann solver consistency and then give results on RSs 
and RSR. 


Definition 2.9. A Riemann solver RS: 0? — L*(R;Q) is consistent if the fol- 
lowing conditions are satisfied for any Ye, Ym, Y, € Q and £ € R: 


RSI, Y I(E) Af $ «6, 


RS|Y+, Ym] (E) em Y i LE 
RSV VID = Ya — ‘ ae 
RS VIE = À ^ Ra 
RS[Y,Y,|(T) f£26, 
() 
RS (Ya, Ya] (5) = Ym —. RS|[Y,Y,(5) = RS[Yo,Um|(§) af E < è (II) 
RS|Y, VIE) = Yin RS Ym, Yr (E) if ERE. 


The advantage of using consistent Riemann solver lies in the fact that it is necessary 
condition for the Lt-continuity of the semigroup associated to the Riemann solver. 


We present then the following propositions, with proofs deferred to Section A. 
Proposition 2.2. The Riemann solver RSR is Lt -continuous and consistent. 


loc 


Proposition 2.3. The Riemann solver RSs is Li -continuous and consistent. 


loc 
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Constrained LWR model 


3.1 Introduction 


In this chapter we describe LWR model (2.1), introduced in Section 2.2, with local 
point constraint on the flow. From the modelling point of view, the unilateral 
constraint on the flow can represent the reduction of car traffic caused, for instance, 
by toll gates, traffic lights or construction sites. The idea of conservation laws with 
unilateral constraint was introduced in the framework of traffic flow in [29] and 
then studied in [2,5,21]. 


3.1.1 The constrained Riemann solver RS;;vr 


We consider a Riemann problem (2.5) for the LWR model 


pe ifa <0, 


3.1 
Bo ifx 0) on 


op F Oxf (p) = 0, p(0, x) = 


with additional condition 
f(p(t,0)) <4, (3.2) 


where py, p, € |0, R] and the value of constraint q € [0, fmax] are constant. Observe 
that the classical Riemann solver RSrwa for the LWR model introduced in Defi- 
nition 2.1 does not necessarily satisfy condition (3.2). For this reason in [29] the 


authors introduced the constrained Riemann solver RSTywyą defined below. 
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Definition 3.1. Constrained LWR, Riemann solver RSiwr : [0, R? + C° (R+; BV(R; 


is defined as follows. 


1. If f(RSvwn|pe prl)(0) < q, then 


RStwrlpe; Prl(€) = RSiwelpe, Prl(Ś): 


2. If f(RSiwrlpe, Pr])(0) > q, then 


RSiwelpe, P(g) 1/6 <0, 


RSiwnlee p. (£) = s Siwr, PIE) if € > 0, 


where 
= minfp € [0,8]: f(p)=ql, p=max{p€ [0, R] : f(p) = a}. 


Notice that p and f coincide if and only if q = fax. On the contrary, if 
q € fmax, namely the constraint is effective, the solution performs at z = 0a 
stationary non-classical shock between f and p. Since f(5) = f(A) = q, the non- 
classical shock satisfies the Rankine-Hugoniot condition (2.7). It can be proved 
that [(t, x) > RS vrlpe, pr|(x/t)| is a weak solution, however it does not in general 
satisfy the Kruzhkov entropy condition [62]. 

The solution obtained by Riemann solver (3.1) can lead to significant increase 


of total variation, as it is shown in the following example: 


Example 3.1. Consider a constant initial datum p(0, x) = po € [0, R] such that 
f(po) > q, hence f(RStwrlpo, Pol(0)) > q. In this case RSiwrlpe, Pr] has a shock 
(Po; P) with negative speed of propagation, followed by a stationary non-classical 
shock (Pp, p) and a shock (5, po) with positive speed of propagation. Its total variation 
jumps from 0 to 2(p — p). 


We recall now the proposition from [29] with the main properties of RS{wr. 


Proposition 3.1. The constrained Riemann solver RS{wp has the following prop- 


erties for any pe, pr € [0, R]: 


^ 


e RSiwrlpe Pr]  RSywnlpe Pr] if and only if (pe, pr) € (p, R] x [0, 9); 
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[(t, c) > RStwrlpe, prl(1/t)] is a self similar weak solution to (3.1), (3.2); 


RS wrlPe. Pr] satisfies the constraint (3.2) in the sense that 


Jim | (RSpwnlpe Prl(r/t)) <q and lim f(RStwęlpe,p.](2/1)) < 4; 


RSTwRlp, Pr] € BV(R; [0, R]); 


RSiwr: [0, RP? + Li, .(R;R) is uniformly continuous. 


3.2 The Constrained Cauchy Problem 


We study now the constrained Cauchy problem 


Op + Or f(p) = 0, p(0, z) = po(z), (3.3) 


with constraint condition 


f (p(t,0)) < a(t), (3.4) 


where z = 0 is the position of the obstacle and q(t) is the maximal flow allowed 
at x = 0 at time t > 0. The definition of weak solution to constrained Cauchy 
problem (3.3), (3.4) is given as follows [29]. 


Definition 3.2. A function p € C? ((0, +20); LS, (IR; [0, R])) is a weak entropy 


loc 


solution to constrained Cauchy problem (3.3), (3.4) if for every test function @ € 
C? (R?; [0, +00)) and for every k € [0, R] we have 


0 < II (le Bl 210+ sm(o — E) (fp) - F(E)) 0,9) dedi 


+ [mto - t 6(0, a) asta fi (1- AL) f(k) ó(t,0) dt 


(3.5) 


and f (p(t,07)) = f (p(t,0*)) < q(t) for almost every t € (0, +00). 


Above the term p(t, 0^) expresses the measure theoretic trace implicitly defined 
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by 
. 1 
ins | i l(t, 2) — p(t,0*)| o(t,z)dedt =0, 
e 0+ € R4 x [0e] 
1 
Ra z lo(t, 2) — p(t,07)| O — 0, 
e>0t € R4 x[—e,0] 


for all ó € CX(R?; R). It has been shown in (5, Theorem 2.2] that both traces 


at x = 0 exist and they are finite. We also point out that Definition 3.2 chooses 


the solution with maximal possible traffic flow through x = 0, however due to 


non-classical shock (5, 6) at z = 0, the solution turns out to be non-entropic [62]. 


Remark 3.1. The first two addends in (3.5) correspond to the entropy condi- 
tion of Cauchy problem (3.3) in the sense of Kruzkov [62]. The last addend 
accounts for constraint (3.4), see [29, Definitions 3.1 and 3.2]. The condition 
f (p(t,07)) = f(p(t,0*)) < q(t) chooses the solution described in Example 3.1 
rather than the constant solution, which satisfies (3.5). Another equivalent formu- 


lations corresponding to unilateral constraints can be found in [5, Proposition 2.6]. 


At last, we introduce the uniqueness theorem for constrained Cauchy problem 
(3.3) (3.4), proved in [29]. 


Theorem 3.1. Assume that f € C9([0, R]; [0, fmax|) satisfy f(0) = f(R) = 0 
and that there exists per. € (0, R] such that (pax — p)f'(p) > 0 for almost every 
p € [0, R]. Assume that: 


1. po € L!(R;(0, R]) is such that sgn(po — pert) (F (Pert) — f(Po)) € BVR; R); 


2. q € BV ((0, +00); [0, Ff (Perit) ])- 


Then there exists unique entropy solution to constrained Cauchy problem (3.3), 
(3.4) in the sense of Definition 3.2. Moreover, p(t,-) € L*(R;[0, R]) is such that 


sgn(p(t, -)— Perit) (f (Perit) — f (p(t, ))) € BV( 
of the Cauchy problem (3.3), (3.4) corresponding to the pair (po, 4) which satisfies 


A 


R). Furthermore, if p is the solution 


the above assumptions, then the following Lipshitz estimates holds true for all 
t € (0, +00) 


lolt -) — P(E, Jen < loo) — Por) lacey + 2 lal) — 4C llos (to o 
(RR) (RR) ((0,41.8) 
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Constrained ARZ models 


4.1 Introduction 


In this chapter we consider the ARZ model described in Section 2.3 with point con- 
straint on the flow. In [46] Garavello and Goatin proposed two distinct constrained 
Riemann solvers. The former is fully conservative, namely both the Rankine- 
Hugoniot conditions (2.14) hold. The latter is non-conservative and only the first 
Rankine-Hugoniot condition (2.14), is satisfied, which ensures the conservation 
of the number of vehicles. In [4] the authors considered the fully conservative 
Riemann solver and proved the existence of the corresponding fully conservative 
solutions to general constrained Cauchy problems for ARZ model with initial data 
of bounded variation and piece-wise constant time depending constraint. The exis- 
tence of non-fully conservative solutions to constrained Cauchy problems for ARZ 
model with constant constraints has been studied in [50], under the assumption 
that the waves of the first family have only negative propagation speeds. However, 
the assumption on the propagation speeds does not allow to take into account 
vacuum states, what is a drawback from application point of view. 

In [42] we proved the existence of non-fully conservative solutions to general 
constrained Cauchy problems for ARZ model with initial data of bounded varia- 
tion, without any restriction on the propagation speeds. The proof is based on the 
Wave Front Tracking method (see [18,59] and the references therein). We provide 
a time decreasing functional Y which allows to give uniform bound for the total 


variation of a sequence of approximate solutions constructed via the Wave Front 
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Tracking method. Then by Helly’s theorem we get the existence of a convergent 
subsequence of the approximate solutions. At last we prove that the limit function 
is indeed a weak solution of the constrained Cauchy problem. These results are 


presented in Section 4.3.1. 


4.2 Description of the model and notations 


In this subsection we introduce notations and the definition of weak solution for 


the constrained Cauchy problem for ARZ model 


OY +0,F(Y) — 0, 
F(U (t, 07)) <q, (4.1) 
Y (0, z) = Y(Uo(x)). 


As in Section 2.3.1, above p = p(t,x) is the density, y = y(t, x) is the generalized 
momentum, p(p) satisfies (2.12) and 


Y = (p,y) € R2 \ {0}, pue (1 = 2) Y, 
U=W"(Y), f(U) =p"'(w— v)v. 


The term U(t,0*) expresses the measure theoretic traces along x = 0. The left 


and right measure theoretic trace are implicitly defined as 


. 1 
lim — 
e—>0+ € 


ff... Itm tol szda =o 
R+ x[—e,0] 


. 1 
lim = 
e30* € 


JI |U (t,£) — U(t,0*)| (é,x)drdt = 0 
R+ x[0,e] 


for all 6 € CX(R?; R). The existence of measure theoretic traces for weak solutions 


guarantees suitable BV bounds, see [31, Lemma 1.3.3]. 


Define the physical domain in the Riemann invariant coordinates 


D = {U = (v, w)" € [0, V] x [0, W]: v < w}, 
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where V and W are the maximal allowed speed and Lagrangian marker, respec- 
tively, with W > V. Denote by Dy = {U € D: v = w} the set of vacuum states 
and by D \ Do = {U € D: v < w} the set of non-vacuum states. 


We introduce the following notation, see Figure 4.1. 


Figure 4.1: Basic notations in (p, f)-coordinates. 


e For any w, € (0, W], the state U (w) = (Te, wę)” is such that p = p^! (we — Te) 


is the maximum of p +> [we — p(p)] p. 


e Let U(q) = (6(q),d(q))) € D be the unique state such that f(U(q)) = q and 
U(w(q)) = Ü(q). We take U(0) = (0,0)7. 


e The propagation speed of the discontinuity (Up, U,) for p, £ p, is 


f(U.) Fu f (Ue) 


ee =a): 


where p(U) = p-'(w — v), as introduced in Section 2.3.1. 
e We define 


max (v € (0, we): v + p(g/v) =w} if we 


> 
v(q) if we < w(q), 
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vods min (v € (0, we): v + p(g/v) = wi) if we > 
v(q) if we S ü(q), 


and denote U^(w;) = (v’(we), we)T, U(we) = (6(w;),w;)". Notice that w + 
: 


v" (w) is non-decreasing and w > 6(w 


> 0 is non-increasing. 
e If v, > 0, then we define U( (ur) = (blu v,)) = (Ur, Ur + p(q/vr)). 


Definition 4.1. Consider the constrained Cauchy problem (4.1) with initial datum 
U, € BV(R; D). A function 


U € BV( Ry x R; D) N C( R+; Lioc( R; D)) 


is a weak solution of (4.1) if Y = W(U) satisfies the initial condition Y (0, x) = 
Y(Uo(x)) for a.e. x € R, for any test function 9 € C2((0, +oo) x R; R) 


ui (p0, + vpO, o) dx dt = 0 (4.3a) 
RLxR 
and if @(-,0) = 0 then 
i (y&ó + vyO.) dz dt = 0. (4.3b) 
R+xR 
Futhermore, the traces of Y at constraint x = 0 satisfy 
f (U(t,0*)) = f (UG,0-)) <q for a.e. t > 0. (4.4) 


We point out that by the BV assumption the weak solutions admit traces in the 


classical sense. 


Remark 4.1. Since we assume that 6(:,0) = 0 in (4.3b) but not in (4.3a), along 
x = 0 weak solutions satisfy the first Rankine-Hugoniot condition (2.17), (with 
o* = 0, namely (4.4)), but not necessarily the second one (2.17),. We underline 
that in [4] the authors consider only weak solutions satisfying the second Rankine- 
Hugoniot condition (2.17), (also along x = 0). Therefore such solutions are weak 


solutions also in the sense specified in Definition 4.1. This is in the same spirit 
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of the solutions considered in [12, 35, 38, 42, 46] for traffic through locations with 


reduced capacity. 


Weak solutions given by the standard Riemann solver RSĄARz (2.3) does not 
satisfy in general constraint condition in (4.1). Thus, let us introduce the following 


sets 


04 = (U U.) € D x D: f (RS arz[Ue, U,]) (07) á q}, 
Q = { (Ue, U.) € D x D: f (RS arz[Ue, U,]) (0^) > q} ; 


We conclude this subsection by recalling from [46] the definition of the non- 


conservative constrained Riemann solver 


CRS Anz DxD E C? ((0, +00); Lie (R; D)) , 
(Us, U,) > CRS Anz|Us, U;]. 


Definition 4.2. For any U,,U, € D, we define CRS arzlUe, Ur] as follows: 


1. If (Ue, U,) € €, then 


CRS arzlUr, U,| (£) = RS arz [Ue, U,] (£). 


2. If(U„U,) € Qe, then 


. J RSanz|Ue. Ü (w;)](£) if € « 0, 
E M UNE) i£20. 


Notice that if (Uj, U.) € Qa, then U(t, x) = CRSarzlUe, U,|(r/t) has a stationary 


^ Z 


non-classical shock (U(wę),U(v,)) with 6(we) < d(v,) = vr. 


4.3 The main result 


In this subsection we present our main result from [42], that is we prove the 


existence of weak solution for constrained ARZ model (4.1). The proof is based 
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on Wave Front Tracking method, which is used to construct a globally defined 
approximate solution. Some technical proofs are postponed to Appendix B. 
We introduce Y: R} — R defined by 


Y(t) = TV (U (t, -);R) + 3TV (w(t, -); R_) (4.5) 


+ 2/TV(a(t);R_) + TV (v(t); R) + W — q(t) + TV. (ii(t, ))]. 


where 
vb) = &(w(t, )). Qm (w(t, ). 
and 
5,2) = w (v(t, z)) if u(t, x) € [6(W),6(q)], 
0 otherwise, 
_ | min (o? (w(t,07)) v(t, 07) } — v(t,07) if (U(t,07),U(t,0*)) € Oe, 
„a 


otherwise. 


Above, TV_ stands for the negative variation. By convention we assume that U 
is left continuous in time; therefore the maps t > TV(v(t,-)) + TV (w(t, -)) and 


t — Y(t) are also left continuous. 


Remark 4.2. Some comments on the definition (4.5) of T are in order. The term 
3TV (w(t, -); R_) +2 [TV (0(2);R_) + TV (v*(t); R_)] is sensitive only for changes 


in w-coordinate, thus it measures the strength of the contact discontinuities in R_. 


Moreover TV_(ń(t)) measures the strength of the rarefaction waves and of non- 
classical shocks (U(t,07),U(t,0*)). At last, y is a compensation term and takes 
positive value if and only if the solution has a non-classical shock. Notice also that 
Y) 0. 


We present now our main theorem from |42]. 


Theorem 4.1. Assume that initial datum Up € BV(R;D) is such that 'Y (0) is 
bounded. Then constrained Cauchy problem (4.1) admits a weak solution U € 
BV (R, xR; D)NC”(R.; L! (R; D)) in the sense of Definition 4.1 and for allt, s > 0 
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it satisfies 
TV (U) < T(0), IU (6, -) — U(s, za < Lit — sl, 


where L = Y(0) - max(V,p  (W)p'(p (W))}. 


Remark 4.3. We underline that TV(Ug) < +00 does not imply that T (0) < +00. 
Indeed, the functions 6 and v? are not Lipschitz continuous close to w = ü(q) and 


(t, x) is discontinuous if (t, x) > v(t, x) crosses the value d(q) or 6(W). 


The proof of Theorem 4.1 bases on Wave Front Tracking method(see [18, 59] 
and references therein) and is described in the two forthcoming subsections. In the 
former one we describe the key tools to construct piecewise constant approximate 
solutions U" to constrained Cauchy problem (4.1), namely a grid, an approxi- 
mate constrained Riemann solver CRS'Rz and the functional I". In the latter 


subsection we give the proof of Theorem 4.1. 


4.3.1 Wave Front Tracking 
Grid construction 


Fix q € (0, fmax(W)). We consider the grid Gn = (W, x Wn) n D, where W, = 


{wi j=—3-J,....J 


LC is a finite subset of [0, W] constructed as follows, see Figure 4.2. 


step 1 Consider the recursive sequence 

w(q) ij =0, 

wi! + p(q/wi*) ifj € N\ {0}, 

which corresponds to the continuous lines in Figure 4.2. In other words, w gives 
the w-coordinate of the intersection point between constraint curve w = v+ p(q/v) 


and v = w/—!. We point out that there exists J < (W — w(q))/p(q/W) such that 


wu «x W < w!*! and by monotonicity of p we have 


W — wit = p(q/w?*) > p(q/W). 
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step 1 
—-- step 2 


-—--— step 3 


vq) 


E idid step 4 
(wit! ) 


PAUSE step 5 


P(w!) wg) w(q) W v 


Figure 4.2: The grid G, with n — 3. 
step 2 Split [w°,w!] into n subintervals of equal length with endpoints 
ROM i € {0,... n}, 
and define recursive sequence 
w? uwiqaqussmq,U.'  ie{0,... n}, jE {1,..., J}, 


corresponding to dashed lines in Figure 4.2. We point out that wi = wit! = wit". 


step 3 In [ü(q), G(q)| we take 
w = w (w), i € {0,... n}, 


corresponding to dot-dashed lines in Figure 4.2. 
step 4 In [ô(w/+1),ü(q)] we let 


773 = 0(w] 3), ie {0,... n}, geli J}, 


corresponding to thin dotted lines in Figure 4.2. 
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step 5 Split (0, ü(w/*1)] into n subintervals of equal length with endpoints 


A EET 
—3— „V(w ; 
up = i ) i € {0,...,n}, 


corresponding to thick dotted lines in Figure 4.2. 
We show in the following lemma that the grid G, is well defined and the distance 
between two points of G, has a lower uniform bound; the proof is deferred to 


Appendix B.1. To simplify notation we write W, = {wi,...,ww} for Wi < wi. 


Lemma 4.1. Fiz n € N and q € (0, fmax(W)). The grid Gn is well defined and 


En = min (w1 — w;) > 0. 
We dee s vm 2 


The cases q = 0 and q = fmax(W) are the straightforward generalizations. 


Approximate Riemann solvers 


To properly define in Gn the approximate solutions U” constructed via Wave Front 
Tracking method, we split rarefactions and introduce the approximate Riemann 
solver S^; : Gn X Gn — PC(R;G,). In more detail, for any (U4,U,) € Gn x Gn 


with wę = w, and v, = w < v, = wy, we define 


U, FES 9 (Us, U1), 
RSARzIUwU;](6) = AU; ifo(Uj,Uj) < € < o(U 5 Uj), 1& j &k— 1, 
U, if> o0(UjaU.); 


where Up = Up, Uk = Up, U; = (wn4j, wę)” and o is defined in (4.2). The corre- 
sponding constrained approximate Riemann solver CRSĄpz: nX Gn > PC(R; Gn) 


is given as follows: 


1. If f(RSKpz[Ve, U,](0*)) < q, then 
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2. If f(RSĄgz[Ue U.](07)) > q, then 


^ 


CRS" no (Ue, U, (6) = RSAralDr. U(we)|(é) a < 0, 
RS'irzlU(vr), U;](Ś) if € 20. 


Interaction estimates 


We introduce now the Wave Front Tracking approximate solutions U” and define 
the map t — ["(t) = ["(U"(t)). To shorten notation we omit the dependence 
on n, that is we write for instance U in place of U” and e for en. Therefore, the 


non-negative map T(U”(t)) has the following form 


T(U(£)) = TV (U(t) + 3TV (w(t);R_) + 2TV (0(t); R-) + 2TV (v(t); R..) 


+2|W-7ft)+ D n(U(62),U(t2*)) 


Above, 5(t) = é(w(t)), v? (t) = v’(w(t)), J(t) C R is the finite set of discontinuity 
points of U(t,-) and 


S OE if w_ = d 6(W) € v- < v, < 9(q), 
n(U_,U,) = (v-) — ù (v4) if w = w, and 6(W) € v- < v, < t(g) 


D & 


otherwise. 


Remark 4.4. We underline that the estimate 


does not hold true in general. For example we can take Ug with a discontinuity 


(U,,U,) away from x = 0 such that the following conditions are satisfied 


We £ Wr, (W) < v, < v, «€ v(q). 
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However we have 


li > F « TV_(5 ; 

lim 9 n (U(t, z ), U(t, x )) < TV ((Uo)) 
xEI(t) 

These considerations show the difference between functionals T and Y defined in 

(4.5). Moreover t ++ Y(t) fails to satisfy the properties listed in Lemma 4.2 and 


thus giving motivation to introduce both functionals I and Y. 


Lemma 4.2. For any n € N let Uy € PC(R; Gn) be an approximated initial da- 


tum, and U” be a corresponding approximated solution constructed via Wave Front 
Tracking method. Then the map t > V"(U"(t)) is non-increasing and moreover 


decreases by at least e, if the number of waves increases. 


The proof of Lemma 4.2 is deferred to Appendix B.2. 


4.3.2 Proof of Theorem 4.1 


We approximate the initial datum Uo with U € PC(R; Gn) so that: 


m |U- Uglnam=0, TUS) < T0). (46) 


li 
n— +00 


lUo Msg S IUollr="(x) » 
By Lemma 4.2 the map t + l'(U"(t)) is non-increasing and as a consequence 
v(t) SPC) e re (Ute TO): 


By the standard application of Helly's Theorem, see [18, Theorem 2.4], we de- 
duce that only finitely many interactions can occur at finite time. Therefore, the 
construction of approximated solution U” can be extended globally in time. 


We begin by straightforward observation that 
|U" (t) Inc) < lUollres qo 
and the Lipschitz condition 


U^ (t, -) — U” (s, -)llt1R) < Lit — s] (4.7) 
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holds for L = Y(0) - max(V,p ! (W)p'(p !(W))). More precisely, if there is no 


interaction between t and s, then 


|| U^ (t, -) = UPS, Jig 
= 9 || — s)o (U (t+, £7), U(t*, £*)) (U(t*, £7) -U(t",z*)) || < Lit — sl. 


1€EJ(t+) 


Above, J(t*) C R is the finite set of discontinuity points of U(t*,-) and o is 
defined by (4.2). We point out that the map t — U"(t,-) is L!-continuous across 


interaction times. 


By Helly's Theorem U” converges (up to a subsequence) to some function U in 


LL.(R, x R; D) as n — +00 and moreover the limit function U satisfies: 


TV(U(.)) < Y(0), 
IU(t,-) - Us, kur < LE — sl. 
lU. ne) < llUollze o: 


It remains to prove that the limit function U is a weak solution of (4.1) in the 
sense of Definition 4.1. 

We deduce that the initial condition of (4.1) is satisfied by (4.6),, (4.7) and 
Li .-convergence of U” to U. To prove that U is a weak solution of Cauchy problem 
(4.1) it is sufficient to show that condition (4.3) holds true, namely for any test 
function ó € C2((0, +20) x R; R) 


ii (paip + vp Ord) dx dt = 0 
R+xR 
and if ọ(-,0) = 0 then 
Jl (y 0,6 + v y Ord) dx dt = 0. 
R, xR 


Clearly Y" = W(U”) and v”Y” are uniformly bounded and we are reduced to show 
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lim Jl (p" dih + v" p" 0,9) dx dt = 0, (4.8a) 
n— +00 R.xR 
and if 6(-,0) = 0 
lini n ri (y^ Ab + v^ y^ 0,6) dz dt = 0. (4.8b) 
n—+00 R4 xR 


We take T > 0 so that g(t, x) = O for every t > T. By the Green-Gauss formula 


we write double integrals in (4.8) as 


| Y GUO 27), U” (t, z*)) AY” (t, £) — AF"(t,2))o(t. x) dt, 
0 aet) 
where 
AY” (t,x) = Y” (t,£) — ACAR 
AF”(t, x) = oO" (U YU WCC Y Que 


Notice that for any classical discontinuity at the point (t,x) Rankine-Hugoniot 
conditions (2.14) is satisfied and 


o (U”(t, £7), U” (t, Y AY” (t,£) — AF” (t, x) = 0. 
If a discontinuity is a (stationary) non-classical shock, then 
c (U” (t, 07), U” (t,0+)) AY” (t, 0) — AF” (t,0) = —AF"(t,0). 


However, in such a case f (U (t,07)) = f(U(t,0*)) and then integral in (4.8a) equals 
zero. By the fact that the non-classical shocks occur only at x = 0 and thanks 
to assumption @(-,0) = 0, we conclude that also the integral in (4.8b) is equal to 
Zero. 


The proof is completed by showing that U satisfies the constraint condition 


f(U(t,07)) < q. 
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To reach this goal we use the fact that both U and U” are weak solutions of 


OY + O,F(Y) = 0, 
Y (0,2) = Y(Uo(x)), 


in R+ x R- and R+ x R+. By the Wave Front Tracking approximation of U” we 
have f(U”(t,07)) = f(U"(t, 0*)) < q for all t > 0. We apply Gauss-Green formula 
to the weak formulation of 0;p” + 0, f (U^) = 0 and then obtain 


af voa > f f (U”(t,07)) v(t) dt 


= ff (m motoeta) + r oso vod?) deat, 


where w € C£ (R+; R) is an arbitrary time-dependent test function with compact 


support and 6 € CX(R; 


R) is a space-dependent test function such that £(0) = 1. 


We pass then to the limit n  +oo in the term 
Ff (OŚW + reo WWE) ardt 
R+ xR_ 
and apply again the Green-Gauss formula. As a result we obtain 


q| w(t)dt2 lim f (U^(t,07)) v (t) at = f f (U(t, 07)) y(t) dt. 
Ry noo Jm, R+ 
Hence the trace f(U"(t,07)) weakly converges to the trace f(U(t,07)) and more- 
over f(U(t,0-)) € q for ae. t > 0. At last, we repeat analogously the above 
procedure to Ry x IR, and obtain /(U(t,0*)) = f(U(t,0*)) € q for ae. t > 0. 
Summarizing, U” can be extended globally in time and TV(U"(t)) < T(t) < 
F”(0) < T(0). 


4.4 A case study 


In this subsection we apply the Riemann solver CRS Anz to describe the evolution 
of traffic through a point constraint representing, for instance, a toll gate. Fix 


Tą X Tp < 0 and wy < wą. Consider in [r4,0) vehicles with zero velocity that are 
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bumper to bumper. We assume that they have different maximal densities p1 = 
p (wi) and pe = p~'(we) in [r4, zp) and [xp, 0), respectively. This corresponds to 


consider constrained Cauchy problem (4.1) with piecewise constant initial datum. 


Ay Wy) dÍd X4 
Qum). PGA TT 


,U2) ifzp<Xz<0, 


( 
Uo(x) = 
( 
( 


W3, We) if > £p. 


Po Py P1 Pr PaP1P2 p x 
Figure 4.3: The solution to Cauchy problem (4.1) with initial da- 


tum (4.9) constructed in Section 4.4. The shaded areas correspond 
to rarefactions. 


Define U; = U(w;) and let 0, = U;(U,, U2) be implicitly defined by 


0, = W — p(P1) = b2, 
see Figure 4.3. Assume that q belongs to (0, fmax(w1)). The solution is constructed 
as follows. We first apply RSarz at A(0, x4) and B(0, xp) and CRS Anz at (0,0). 
Two stationary contact discontinuities CD4 and CDp start from A and B. From 
(0,0) we have a backward rarefaction Ro, a stationary non-classical shock NS, 
and a forward contact discontinuity CDs. Denote by C and E the first and last 
interaction points between Ry and CDg. Denote by D and F the first and last 
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(x,t) > v(t, z) (x,t) > w(t, z) 


Figure 4.4: The solution to Cauchy problem (4.1) with initial 
datum (4.9) constructed in Section 4.4. Darker colors correspond 
to higher values. 


interaction points between Ry and CD4. Notice that Ro crosses CD g; on the other 
hand, Rọ does not cross CD, and expires. Both CD4 and CDg have speed of 
propagation 0, = Ô after their interaction with Ro. Once CDg reaches z = 0 
at G = (ta,0), a backward rarefaction Ra is created and the left state of NS; 
varies from Üo to Ü. Let H and I be the first and last interaction points between 
CD and Rg. Observe that Rg does not cross CD 4 and expires. CD4 moves with 


velocity ô; after time t = tz, reaches z = 0 at L = (tz,0) and then continues in 
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(0, +00). CDo and CD, do not interact because their speeds of propagation are w 
and 9», respectively. 

The solution above has the following physical interpretation. At time t = 0, 
the rightmost vehicle starts to move with constant maximal speed w». The other 
vehicles start as soon as the following distance in front of them is safe. This 
acceleration is related to the rarefaction Rọ. Due to the presence of the toll gate, 
which hinders the flow at z = 0, the vehicles initially in [rg, 0) stop to accelerate 
once they reach the velocity d and flow q, which is the maximal capacity of the 
toll gate. On the other hand, the vehicles initially in [r4, xg) stop to accelerate 
once they reach the velocity 0; = à of the other vehicles in the upstream of the toll 
gates. Once all the vehicles initially in [x8,0) have crossed the toll gate, namely 
at time t = tą, the upstream vehicles accelerate and reach the velocity 6; and flow 
q. This acceleration corresponds to the rarefaction Rg. Finally, after time tz all 
the vehicles have passed the toll gates. 

In Figure 4.3 we represent the initial datum and an overall overview of the 


solution corresponding to 
w=1, w=6/5,  p(p—p q-—v3/5  z4—-8, xp = -5. 


The quantitative evolution of the corresponding solution is represented in different 
coordinates in Figure 4.4. We finally observe that, once the overall picture of the 
solution is known, it is possible to express in a closed form the time at which the 


last vehicle passes through x — 0, indeed 


tę = [(zg — z4)p I (wi) — rep 1 (wz)]/q ~ 24.4716. 
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Constrained PT models 


5.1 Introduction 


In this chapter we consider the constrained versions of the PT* and PT? models 
introduced in Section 2.4.3 and study both the cases with and without metastable 
phase. We first introduce and describe in full details the corresponding con- 


strained Riemann solvers. We then study their main properties we obtained in [35], 


1 


ioc Continuity with respect to the initial data and their invariant domains. 


namely L 
Moreover, we show our exsitence result [10] for the constrained PT? model with 
metastable phase. 

As in Section 2.4.3 the Riemann solvers associated to the Riemann problem 
(2.20), (2.21) are denoted by RSR and RSs, respectively in the cases of intersecting 
and non-intersecting phases. 

Besides the initial condition (2.21), in this chapter we enforce the local point 


constraint on the flow at z — 0, namely 


f (U(t,07)) <q, (5.1) 


where q € (0,f;) is a given constant quantity. Admissible solutions to (2.20), 
(2.21) do not satisfy in general (5.1). Thus we introduce the following definition 
of admissible constrained solution to (2.20), (2.21), (5.1). 


Definition 5.1. Fix Up, U, € Q. A self-similar function U = (v,w): Ry xR > Q is 


an admissible constrained solution to constrained Riemann problem (2.20), (2.21), 
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(5.1) if for UT = U(t,07) and U* = U(t,0*) the following conditions hold true: 


e The maps 


U(t,x) ware UE if x « 0, 
(t,x) > (t, 2) > 
U- if x 0; UU ay ifx>o0, 


are admissible solutions to the Riemann problems for (2.20) with Riemann 


data respectively given by 


U, ifx <0, Ut ah gee dp 
U(0, z) = U(0, z) = 
U- Pec, U, PX 


in the sense of Definition 2.6. 


e The constrained condition f(U~) = f(U*) < q is satisfied. 


Remark 5.1. Notice that the linearized momentum y is in general not conserved 
across possible stationary discontinuities at x = 0. For this reason in the above 
definition we cannot impose y(U~)u~ = y(U*)v*, even if U-,U* € Q.. Therefore 
an admissible constrained solution to Riemann problem (2.20), (2.21), (5.1) taking 


values in €, is not necessarily a weak solution. 


We denote by CRSR and CRSs the constrained Riemann solvers associated 
to the Riemann problems (2.20), (2.21), (5.1), respectively for the cases with 
and without metastable phase. Their definitions given below are in accordance 
with Definition 5.1, in the sense that (t,x) > CRSg|U;, U,|(x/t) and (t,x) + 
CRSs|U;, U,|(x/t) are admissible constrained solutions. 

We let (with a slight abuse of notation) 


CRSR=RSR in =D, = {(Ue,U,) e X : f(RSR[Ue, U,|(t,0*)) < 
CRSs=RSS in D, = (Us, U,) € Q? : f (RSs[U;, U. tt, 0*)) < 


and we denote Dy = Q? \ Dj. 


We will also discuss their main properties, such as consistency, Lt 


joc Continuity 


and invariant domains. In this regard, we recall the following definition. 
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Definition 5.2. An invariant domain for a Riemann solver RS: Q? + L*(R; Q) 
is a set T C Q such that RS|Z,T|(R) CT. 


Proposition 5.1. Let RS: Q? > L*(R;Q) be a Riemann solver. If (t,£) + 
RS|U,, U,|(r/t) satisfies (5.1) for all U,,U, € Q, then RS does not satisfy (I) of 


Definition 2.9, hence it is not consistent. 


Proof. Fix Up = (Vmax; W+) = U,, see Figure 2.2. We have that f(U,) = ff >q. 
By the finite speed of propagation of the waves we have RS[U;,U,|(4) = U, for 
some € > 0. Let Um = U,. Assume by contradiction that RS[Um, U,](£) = Um 
for any € < €. Then f(RS[Um, U,](0*)) = /(U,) > q and this is contradictory to 
(5.1). We conclude that the property RS[Um, U,](x) = Um for any € < € in (I) is 


not satisfied. 


We define below the constrained Riemann solver for the case with metastable 


phase. 


Definition 5.3. Assume that Qs N Qe ź 0. The constrained Riemann solver 
CRSg: Q? — L® (R; Q) associated to (2.20), (2.21), (5.1) is defined as 


RSg|[Us, U,](£) if (Ua, U,) € Dy 
CRSr|U:, U,](E) = ^ ; 
RIVE RSR Ur, UJE) if E < 0, if (Us, U,) c D 


where Ü — Ü (we, q) E€ Qe and U = Ü(v,,q) € Q are uniquely selected by the 


conditions 


Vmax if f(y (U,)) > q, 
v. Y fz (U.)) Sq. 


fU)=f(l)=q w=max{w,w_}, ŭ= 


In Figure 5.1 we present all possible choices of U and U given by the above 
definition. It is worth noting that U and U satisfy the following properties: 


e if (U,,U,) € Do, then w > w and v, > Ô, 
e if (Us, U,) € Dz and U, € Q;, then ù = w_, 
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f 
Ur 
U; 
q t- 
Rp 
(Ue, Up) € 0? (Ur, Ur) e 02 
f £ 
q -Y Ur 
U 
q 
R p R p 
(Ua, Up) € Re x OF (Ue, Ur) € NF x A 


Figure 5.1: The selection criterion for U and U given in Defini- 
tion 5.3 in the case (U,,U,) € D». In the first picture U}, U? 
represent the left state in two different cases and Us. Û, are the 
corresponding Û. In the second picture qı, da are q in two dif- 
ferent cases, Uh, U> are the corresponding U and Ü,, U, are the 
corresponding U. 


e if (U,,U,) € D; and U, € Qg, then © = Vmax- 


We define below the constrained Riemann solver for the case without metastable 


phase. 


Definition 5.4. Assume that Qt N Qe = 0. The constrained Riemann solver 
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RP 
Case (5.2). Case (5.3). Case (5.4), v1 = 1(Uz). 


Figure 5.2: The selection criterion for U and U given in Defini- 
tion 5.4. 


CRSs: Q? + L*(R; Q) associated to (2.20), (2.21), (5.1) is defined as 


RSs[Us, U.](£) if (Ue, U.) € Di, 


CRSs[Ue, UAE) = À | mssU, ÔE if& <0, 


RSs[|U,U.(8) #67 0, 


if (U, Dx) e D», 


where U = U (we, q) EQ, and U = Ü(w;,v,,q) € Q are uniquely selected by the 


conditions 
f(U) = f(U) = max{ f(U) <q:UEMh, w= max{we, w_}}, 


vae WIO U 
w YFU <q. 


w = max{ w, w}, ù = 


Remark 5.2. In the case (U;, U.) € D», the choices of Ü and Ù in Definitions 5.3 


and 5.4 coincide if and only if one of the following conditions is satisfied 


(U, Ur) € OF x Op, q € (f, F(U), (5.2) 
(Ue, Up) € QF x Qr, q € (f(vi(Ue)), F(U), 5.3 
(Ue, Ur) € Qe x Qr, q € (F(YE(Ur)), F(V1(U2))), (5.4) 


and in this case f(U) = f(U) <q and ô = ve. The above cases are considered in 
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Figure 5.2. Thus CRSR(U,,U,]| Æ CRSs|U,,U,] if and only if (Uc, U,) belongs to 
Dı and satisfies one of conditions (2.22), (2.24), (2.23), or (Ue, Up) belongs to D» 
and satisfies one of conditions (5.2), (5.3), (5.4). 


In the following propositions we give the main properties of the constrained 
Riemann solvers CRSR and CRS». We point out that by Proposition 5.1 the 
constrained Riemann solvers CRSR and CRSg are not consistent; for this reason 


we consider in the next two proposition only condition (II) of Definition 2.9. 


Proposition 5.2. CRSp is LL, -continuous and satisfies (IT) of Definition 2.9. 


loc 


Proposition 5.3. CRSg is not LL,,-continuous and satisfies (IT) of Definition 2.9. 


loc 


The proofs of the above two propositions are deferred to Sections A.3 and A.4, 
respectively. 
We summarize in Table 5.1 the Li 


loc 


Riemann solvers considered in Section 2.4 and in the present chapter. We underline 


-continuity and consistency properties of 


that the lack of L} -continuity of a Riemann solver does not prevent the study of 


the general Cauchy problem, see for instance [3, 27]. 


metastable phase | Lt -continuity d 
RSR Yes Yes Yes | Yes 
RSs No Yes Yes | Yes 
CRS Yes Yes BLE 
CRS: No No „BASIĘ 


Table 5.1: The main properties of the Riemann solvers. 


We complete this subsection with some remarks on the invariant domains. By 
construction €) is an invariant domain for both CRS; and CRS». The spaces Qg 
and Q. are invariant domains for CRSs, however not for CRSg. Therefore we 
look for the minimal (w.r.t. inclusion) invariant domains for CRSp containing Q 


or Qe, see Figure 5.3. 


Proposition 5.4. Let CRSR be the constrained Riemann solver introduced in Def- 


inition 5.1. 
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(ICR.1) The minimal invariant domain including Q; is Te = Q¢ UL, UT, where 


T = (UEQ.: f(U) £q< FZ (U))], 


je fu E Qe: F(U) > a, LEE (PW) > o}. 


(ICR.2) The minimal invariant domain including Q is 
Te = Qe U $ (Gua (Ginx) F- 


The proof is deferred to Subsection A.4. 


R p R p 


Figure 5.3: Invariant domains Zp (left) and Z, (right) described in 
(ICR.1) and (ICR.2) of Proposition 5.4. 


We recall that the first Lax curves are strictly concave for the PT? model and for 


Ly 
the PT* model under the assumption 7 e (re) < 0. This implies Z = 0. At 
p 


last, we give the minimal invariant domains for CR s containing Qg or Qe. 


Proposition 5.5. Let CRSs be the constrained Riemann solver introduced in Def- 


inition 5.4. 
(ICS.1) The minimal invariant domain containing €» is 


Ty defined in (ICR.1) of Proposition 5.4 ifq < ft 
QUEUE: v = ve} if q 2ft. 


? 


f= 
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(ICS.2) The minimal invariant domain containing Qe is 


Qe U La vaa Q(q/Vmax)) } if q < tes 
Qe ifq2 fy. 


Lo 


Since the proof is analogous to that of Proposition 5.4, we skip it. 


5.2 Existence result for constrained PT? model 


with metastable phase 


In this section we give our existence result [10] for the constrained PT? model with 
metastable phase QF = Q; NA. For convenience we use the Riemann invariant 
coordinates U = (v,w). We adapt the notation introduced Subsection 2.4 to the 


U-coordinates, see Figure 2.2. 


5.2.1 Notations, definitions and main result 


In this subsection we reformulate in the U-coordinates the main assumptions on 
the parameters, adapt accordingly useful notations introduced in Subsection 2.4, 
see Figure 5.4, give the definition of solutions, state the main result in Theorem 5.2 


and at last introduce the Riemann solvers. 


Figure 5.4: Notations introduced in Subsection 2.4 adapted to 
U-coordinates. 
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We consider the constrained Cauchy problem for PT? model (2.20) 


Free flow Congested flow 
U = (v,w) € Qr, U = (v,w) € (X, (5.8) 
OU) + sf (U) = 0, Ou(U) + ôf (U) = 0, 
ie e du + vO,w = 0, 
U (0, x) = Uo(x), (5.6) 
f (U(t.07)) <a, (5.7) 


where p(U) = p !(w — v), f(U) = vp(U), and q € [0,f*] is a given constant 
quantity. Above, see Figure 5.4, in analogy with the notations introduced in Sub- 


section 2.4, we define 


Qg = {U = (v, w) : V = Vmax, W- — 1 S w S w4}, 
OF = {U = (v, w) : V = Vmax, W- —1<w< w}, 
OF = {U = (v, w) : V = Vmax, W- S w < w+}, 

Q. = {U = (v, w): 0 sus Vmax, W- S w S w4} 


We assume that requirements (H1) and (H2) are satisfied, therefore the first 


Lax curves are strictly decreasing and concave. The p function satisfies (2.12) 


p(0)=0, pp)>0,  2p(p)-p'(pp»0  foreveryp>0. (5.8) 
For later use we introduce the map W: 2 + |w.., w,] defined as 


W(U) = max{w_,w}. 


Moreover, see Figure 5.5, let vy € [0,Vmax] and w, € [w_ — 1,w,] be defined as 
follows: 
ifq=f": Vo = Vmax; Vo = Vmax, Wq = W4, 


EPEITI W -VN (1) — W4, W =p ) + Vmax 


q 
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OSN WNE 
Wom |. 
> 
p 
Figure 5.5: Representation of wq, vy and hq in the case q € (0,f ). 
The curve in the figure on the left is the graph of h, and corre- 
sponds to the horizontal solid segment in the figure on the right. 
In particular w, = h,(v.). 
1 = T q = q : q 
ifq e 3. v exe) on V; e(t) =w wQ— w--l-r- 
ifq=0: v =0, v; =0, w, — w- — 1. 


For any q € (0,f*), let hy: |v; , v] — [w-,w4] be given by h,(v) = v + p(q/v), see 
Figure 5.5. Notice that h, is strictly decreasing by (H1) and is strictly convex by 
(5.8). 


5.2.1.1 The constrained Cauchy problem 


The notion of solution to Cauchy problem (5.5) is obtained by combining the 
notion of solutions of the LWR and ARZ models, and by choosing which phase 
transitions are admissible, see |11]. We recall that a discontinuity between states in 
Q is entropic if and only if its speed of propagation is Vmax. Therefore we consider 


only the generalized entropy-entropy flux pair introduced in [4] for the ARZ model, 


0 if v È k, 
E'(U) = p(U) 


—1 ifv<k, 
p! (W(U) — k) 
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where U € Q and k € [0, Vmaxl. 
We introduce now the definition of solution to Cauchy problem (5.5), (5.6). 


Definition 5.5. Fir U, e BV(R;Q). We say that a function 


U € L™((0, +00); BV (R; 2)) n C?(R,; Li (R; 0)) 


is a weak entropy solution to Cauchy problem (5.5), (5.6) if the following conditions 
hold: 


p 


(S.1) Initial condition (5.6) holds for a.e. x € 


(S.2) For any test function ó € CX ((0, +00) x R; R) we have 


JL mme 102.0 (Jaen - (7). 


(5.3) For any k € [0, Vmax] and $ € C?*((0, +20) x R; R) such that $ > 0 we have 


ik | (E'(U)8,0 + Q*(U)8,9) dx dt > 0. 
R+xR 


In [11, Theorem 2.8] the authors proved the existence of a weak entropy solution 


satisfying the above definition. We recall their result below. 


Theorem 5.1. Assume that initial datum Ug belongs to L! nBV(R;Q). Then 
Cauchy problem (5.5), (5.6) admits a weak entropy solution U in the sense of 


Definition 5.5. Furthermore there exist two constants C and L such that for any 
t,s 20 


TV(U(t)) < TV(Uo), lU(tllu-o; € C. UO — U()llanquo < LÉ — s]. 


In the following definition we introduce the notion of solution to constrained 
Cauchy problem (5.5), (5.6), (5.7). 
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Definition 5.6. Fiz U, e BV(R; Q). We say that U € L®((0, +00);BV(R;0))n 
C? (R,;LL,(R;O)) is a weak entropy solution to constrained Cauchy problem 
(5.5), (5.6), (5.7) if the following holds: 


(CS.1) Initial condition (5.6) holds for a.e. x € R. 


(CS.2) For any ó € C%((0, +oo) x R; R) we have 


[oo r3.) dzat = o (5.9) 
— 
and if o(-,0) = 0 then 


nac : U)8.ó + f(U)0,9)W(U) dadt = 0. (5.10) 


(CS.3) For any k € [0, Vmax] and 6 € CX((0, +20) x R;IR) such that ¢(-,0) = 0 


and à > 0 we have 


i (E*(U)2,6 + Q*(U)O,@) dxdt > 0. (5.11) 


(CS.4) Constraint condition (5.7) holds for a.e. t > 0. 


In the following proposition we state which discontinuities are admissible for 
the solutions to (5.5), (5.6), (5.7). 


Proposition 5.6. Let U be a weak entropy solution to constrained Cauchy prob- 
lem (5.5), (5.6), (5.7) in the sense of Definition 5.6. Then U has the following 


properties: 


e Any discontinuity ofU has speed of propagation o satisfying the first Rankine- 


Hugoniot condition 


e(t) (p(U(t.5(0*)) — p(UlŁ,607))) = FUE O) — KOO), 
(5.12) 


and away from x — 0 it satisfies also the second Rankine-Hugoniot condition 
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e(t) (e(U(t.5()*)) W(U(Ł,6(2)*)) — p(U(Ł,5(4)7)) W(U(Ł,óC2)”)) ) 
= f(Ut.5()*)) w(ut.à()*)) — FUE.) W(U(L,6(2)7)). (5.13) 


e Any discontinuity of U away from x = 0 is classical, namely it satisfies the 


Lax entropy inequalities. 


e Non-classical discontinuities of U may occur only at x — 0, and in this case 
the (density) flux at x — 0 satisfies (5.7). 


Proof. These properties follow from (CS.2), (CS.3) and (CS.4). We underline 
that (5.12), (5.13) are equivalent to 
(v(t,0*) — e())e(U(t,0*)) = (v, 0) - oft))p(U(Ł,07)). 
(w(U(t,0*)) - w(U(t,07)) ) (v(t, 0") = o(t))p(U(Ł,07)) = 0. 
The above equalities are satisfied by phase transitions and shocks because for 


them W(U(ż,0*)) = W(U(,07)), and by contact discontinuities because for them 
e(t) = v(t, 07. 


Remark 5.3. Since we assume that 6(:,0) = 0 in (5.10) and in (5.11) but not 
in (5.9), along x — 0 weak entropy solutions satisfy the first Rankine-Hugoniot 
condition (5.12) (with o(t) = 0), but not necessarily the second one (5.13). This 


is in the same spirit of Remark 4.1. 


Let [w_ — 1,w,] > w > ü(w,q) = ($(w,q), $(w,q)) € Qe and [0, Vmax] > v > 
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U(v,q) = (v(v, q), x(v, q)) € Q be defined by, see Figures 5.6 and 5.7, 


(&(w, q), (w, q)) = 


+ . 
vi,W_) if w < w S w, 


(5.14) 
Vmax, wg) ifv 7 vi, 


(šv, q) «(v q)) = 4 (v, hą(v)) ifv € z, vi], 


( 
( 
(Vmax: Wa) if w < wg, 
( 
( 
( 


v7, w4) if v < v7, 
where w= wo U and Ww = w o U. As a consequence 


=p" ) - $(w, q)), 
—*(w,q)) ifv <vq, 


> ifv > V 


Remark 5.4. By definition we have f (Ü(w, q)) = f (U(v, q)) — q. Furthermore, 
w > U(w,q) and v + U(v,q) are continuous if and only if q > f-, and in this case 
they are Lipschitz continuous. On the other hand, if q < f~, then w > Ü(w, > 
and v ++ U(v,q) are only left-continuous. Moreover &(w,q) > w and (v, q) > 

At last, w > w(w,q) and v > v(v,q) are non-decreasing, while w +> (w,q) and 


ut W(v, q) are non-increasing. 


Let TV, and TV_ be the positive and negative total variations, respectively. 


For any U: R > Q we let 


TU) = TV4 ($(w, q); (—0,0)) + TV. (&(w, q); (—00,0)), (5.15) 


= TV. (v(v, q); (0, +00)) + TV- (&(v, q); (0, +oo)). 


d 
S 
\|- 


For all U € Q and k € [0, Vmax] we define 


q F3], 
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Va Vmax V 


Figure 5.6: The functions U and U in the case q € (f~, f*) defined 
in (5.14). 


where 


w ifw>0, 


0 otherwise. 
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Vmax U 


Figure 5.7: The functions U and U in the case q € (0,f~) defined 
in (5.14). 
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In the following theorem we give our main result obtained in [10]. 


Theorem 5.2. Let Uj) € L'NBV(R;Q) and q € [0,f*] satisfy one of the following 


conditions: 


(H.1) qe [f-,f*]; 
(H.2) q € [0,f-) and T(Uo) + Y (Uo) is bounded. 


Then the approximate solutions U, constructed in Section 5.2.4 converge to a weak 
entropy solution U € C9(R,; BV(R;Q)) of constrained Cauchy problem (5.5), 
(5.6), (5.7) in the sense of Definition 5.6. Moreover there exist constants C, and 
L; that depend on Uo and q, such that for all t, s > 0 


TV(U()) < Cy, OE) - Ulli € LE-s| IU Olli) < Roa 

(5.16) 
Furthermore, non-classical discontinuities of U can occur only at x = 0, and if for 
any k € [0, Vinax] and ó € CĘ((0, +20) x R;R) such that ¢ > 0 we have 


T 


lim NK (U, (t, 07 )) b(t, 0)dt = ye NË (U(ż,07))o(t, 0)dt, (5.17) 


n— +00 0 
then the (density) flow at x — 0 is the maximal flow q allowed by the constraint. 
The proof is deferred to Section 5.2.4. 


Remark 5.5. If q € [f~,f*], then both w > Ü(w, q) and v > U(v,q) are Lipschitz 
continuous, hence T (Uo) + T(Uo) is bounded for Up with bounded total variation. 


5.2.2 The constrained Riemann problem 


This subsection deals with constrained Riemann problem (5.5), (5.7) with initial 


condition 


U, ifr«0, 
UO x) = (5.18) 
U, dig o0. 


Below we recall the definition of the constrained Riemann solver CRS that we 
introduced in [35]. 
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Let RS? be the Riemann solver for PT? model given in Definition 2.7. We 
point out that in general (t,x) — RS%[U,,U,](x/t) does not satisfy constraint 


condition (5.7). This motivates the introduction of the sets 


Dı = 4(ULU,)EQNxQ: f(RSZ(U„U,]|(0*)) < q} 
= {(Up,U,) € % x Qg: F(Ue) <q) 
U { (Ue, Ur) € Qe x Q : f((ur, (W(UZ))) < a} 
U{ (Us, Ur) € AF x Qg : mint f(Ue), f(v (U.))] < a); 
D= Q Dr 


and the constrained Riemann solver CRS}, given below. 


Definition 5.7. The constrained Riemann solver CRS: Q + L®(R; Q) associ- 
ated to constrained Riemann problem (5.5), (5.7), (5.18) is defined as 


RSplU;, Ee) if (Us, U,) € Dı, 
CRSR[Ur, U.](£) = RS^[U;, Ue] (6) if E <0, 


: if (Ue, Ur) € Do, 
RSU UE if € > 0, 


where U, = U(wę,q) € Qe and U, = U(v,,q) € Q are given in (5.14). 


Clearly, the above definition corresponds to Definition 5.3 for the PT? model. 
In Figure 5.8 we clarify the selection criterion (5.14) for Up and U, in the case 


(Up, U,) € Dz. We point out that 0, and Ù, satisfy the following general properties. 
If (Ue, U,) € D», then we > w, and v, > dy. 
If (Uz, U,) € Ds and U, € Q;, then & = wL. 
If (Ue, U,) € Dz and U, € Qr, then V, = Vmax- 


We recall that the maps (t, 7) — SU; U,](x/t) and (t, x) > CRS%[U,,U,](z/t) 
are solutions to Riemann problems (5.5), (5.18) and (5.5), (5.7), (5.18) in the sense 
of Definitions 5.5 and 5.6, respectively. Moreover both S^, and CRS% are Li, 


loc” 


continuous, see Propositions 2.2 and 5.2. 
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Figure 5.8: The selection criterion for Ue = U(we,q) and U, = 
U(v,,q) in Definition 5.7 for (U,,U,) € Dz and q € (0,f^). In the 
first picture U}, U? represent the left state in two different cases 
and 01, U? are the corresponding Uy. Analogously in the second 
and third pictures for U}, U? and U!, Ù. 


5.2.3 A case study 


In this subsection we simulate the traffic across a toll gate located at z = 0 and 
with capacity q by applying model (5.5), (5.6), (5.7). Let w_ and w, be the 
Lagrangian markers of vehicles that are at time t = 0 at rest in |ra4, rg) and 


[z p, 0), respectively. The corresponding initial condition is 


U, ifr€[ra,rp), 
Uo(z)=4U, if x € [xp,0), 
Uo if x € R\fza,0), 
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with Uo = (0, Vmax), Us 
fa 


(p !(w..),0) and U, = (p~'(w,), 0), see Figure 5.9. 


uo U; U,p 


Figure 5.9: The solution constructed in Section 5.2.3. 


The solution is constructed below by solving Riemann problems at the discon- 
tinuities of Up and by considering the interactions of the waves between themselves 


or with the point constraint x = 0. Put 


A Z 


Ue = U(w_, q), U, = Û(w4, q), U = U(Vaas; q), U, = (r, W(U,)). 


At x = 0 we apply CRS% and obtain a backward rarefaction Ro(U,, D, ), a station- 
ary non-classical shock NSo(U,, U) and a forward contact discontinuity CDo(U, uo), 
which moves with speed Vmax. At x = xg we apply RS? and obtain a stationary 
contact discontinuity CDp(U;,U,). Let C and E be the starting and final interac- 
tion points between CDg and Ry. During such interaction CD accelerates, while 
Ro crosses CDp and changes its values. After time t = tg, CDg has speed Ÿ, > 0 
and interacts with NS at G. At G we apply CRS% and obtain a backward rar- 
efaction Ra(u., U;) and a stationary non-classical shock NSa(Uy, U). 

At x = z4 we apply RS% and obtain a stationary phase transition PT 4(uo, Ue). 
Let D and q be the starting and final interaction points between PT and Ro. 
During the time interval (tp,t,) we have that PT4 accelerates and Ro starts to 
disappear. After time t = tp we have that PT 4 has speed v, > 0. Let H and I be 
the starting and final interaction points between PT4 and Rg. During the time 
interval (£j, tr), P'T 4 accelerates and Rg starts to disappear. After time t = ty, 


PT4 moves with speed we > 0. Finally, PT 4 interacts with NSg at L and then 
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moves with speed Vmax- 


In Figure 5.10 we show the solution corresponding to p(p) = p? and to the data 
ta=—8, zpg——5, W=l, w,-0/5, Vmax=3/5, q=V3/5. 


Such solution is obtained by the explicit analysis of the wave-fronts interactions 


with computer-assisted computation of the interaction times and front slopes. 


td 


L 


my | 


LZ 
SĘ 


ZZ SĘ j 

E 

RZ j 
DU 


(t,£) > v(t, z) (t, z) > w(t, z) 


Figure 5.10: The solution constructed in Section 5.2.3. 


Once the overall picture of the solution is known, it is possible to express in a 
closed form the time at which the last vehicle passes through x = 0, indeed 


f, = EB TAPE TBPL ogg. 


q 
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5.2.4 Proof of Theorem 5.2 


In this subsection we prove Theorem 5.2. The proof consists of the following steps. 
We start by constructing grid G, and defining approximate Riemann solvers S7", 
CRSĘ” in Ga. Therefore, we construct approximate solutions U, = (vn, wn) via 
wave-front tracking method [18,59] to constrained Cauchy problem (5.5), (5.6), 
(5.7) in the space PC of piecewise constant functions. We show that U, can be 
defined globally in time; this goal is achieved thanks to a non-increasing Temple 
functional T,, which decreases each time the number of discontinuities of U, in- 
creases. Next, we prove that U, converges to function U, which is a solution to 
(5.5), (5.6), (5.7) and satisfies the estimates listed in (5.16). At last we consider 
the flux density of the non-classical shocks. 


Fix n € N sufficiently large. We simplify the notation by letting 
‘We = We, Ü, = U(we, q), U; = Ü(we, q) 
and so on, where à and à are defined in (5.14). 


The grid 


Define G, = 0. NP, see Figure 5.11, with P given in the (v, w)-coordinates by 


(UET iv) x (US {w}, 
where M, N, v! and w”, are defined as follows: 


e If q = 0, then we let M = 1, N = 2, 
v = i2 "Wmas ifi e T0 62" > 


w.-ilqi27^ ESET Oe; 
w_+(é—2")2-"(w, —w_) ifiEf27+1,...,27+1). 
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e If q € (0,f”), then we let M 23, N =3, 


12 7v, if i € {0,...,2"}, 
z —1/,,,4:27—i ie à n n+1 
v = $ Rh C ) ifi € 427 +1,..., 2T}, 


esan a n cui) appe aae, 


w_-1+i2 (wą — w- +1) if i € T0 ssa. 
w= 4 w, + (i — 2”) 27"(w.. — wy) ifie (2^ c 1,...,2^H), 
we (1 =2-27)2%(w4,=wo) ifeef2" 1 +1,...,3: 2°}. 


e Ifq € [f7,ft], then we let M 22, N=3, 


52 y E1640: 5. h 
hz p ifie{2 ALL 2t}, 


w_—1+i12™ ifi € {0,...,2"}, 
w = $w-+(i—27%)2(w,=w-) if ie {2" +1,..., 2°}, 
Wy + (i — 27%11)27*(w+=w,) ifae {2741 +1,...,3-27%). 


Notice that if q € (f^, f^], then not necessarily wt Z witt. 


The approximate Riemann solvers 


To properly define in G, the approximate solutions U” to (5.5), (5.6), (5.7) con- 
structed via wave-front tracking method, we introduce the approximate Riemann 
solvers RSR”, CRS” : Gn x Gn — PC(R; Gan) obtained by discretizing the rar- 


efactions. In more detail, for any (Uj,U,) € Gn x Gn such that w, = w, and 
h+k 


ve =v" < v, = v"t*, we define 


U, if € X c (Ug, U1), 
RSR [Us U(E) = AU; if o(U,_,U;) «€£&o(UsUja), 1&3 &k- 1, 
U, MEE ouis] 
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Figure 5.11: The grid G, for q € (0,f^) and n = 2. The curve in 
the figure on the left is the support of h,, which corresponds to (a 


portion of) the ho 


rizontal line in the figure on the right. 


where Up = Ur, Ui = U, and U; € Gn is such that v; = v^? and w; = we. The 


approximate constrained Riemann solver CRS%” is defined as follows: 


1. If f (RS? [U,, U, (0 


) < q then 


CRS?” Up, U,] = RSZ” [Us U,]. 


2. If f (RSR [Ue, U,|(07 


CRS (Ue, Ur|(E) = 


The approximate solut 


We give below the construction of the approximate solution U, € PC( 


-)) > q then 


RSR [Ue UIE) if E< 0, 
RSP"U. U, (E) if £ 20. 


ion 


Ry x 


to (5.5), (5.6), (5.7). As a first step we approximate U, with Up, € PC( 


R; Gn) 


R; Gn) 
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such that 
|on||= < [v||L=, [wat < [w||r=, 
TV(v,) <S TV(v), TV(w,) < TV (w), 
| (5.19) 
Hm feel. = 0, Jim fus — wil, = 0 
T(Uss) € CT (Un), T(Uon) < CT (Un), 


for a constant C. The approximate solution U, is then obtained by a wave-front 
tracking method, which exploits CS" at z = 0 and RS” away from x = 0. 
We prove that only finitely many interactions may occur in finite time; this ensure 
the global (in time) existence of U,. Then we show that U, € PC(R, x R; Gn). 


At last we demonstrate that U, converges (up to a subsequence) in Lie 


to a limit 
U and show that U is a constrained weak entropy solution to (5.5), (5.6), (5.7) in 
the sense of Definition 5.6. 


The Temple functional 


Notice that any contact discontinuity (CD) has non-negative speed (of propaga- 
tion), any shock (S) or rarefaction shock (RS) has negative speed, all the non- 
classical shocks (NSs) are stationary and the speed of all the possible phase tran- 
sitions (PTs) ranges in the interval (—f~ /(p~'(w_) — r), Vmax). Below we say that 
(Ue, U,) is a null wave if Up = U,. Notice that if (U,,U,) is a PT then U, € Qy and 
U, € QZ, moreover if (Ue, U,) is a PT with w, > w_ then p, = 0. 


Denote by (t) the number of discontinuities of U, (t, -) and introduce Th: Ry > 
R, defined as 


A x 


T„(t) = TV (vn (t, -)) + TV(wn(t,')) + 2T(t) + 2T,(t), 


where T,(t) = T(U,(t,-)) and T„(t) = T(U,(t,-)). Conventionally, we assume 
that U, is left continuous in time, i.e. U,(t,:) = U,(t.,-). Then also T, is left 
continuous in time. By the monotonicity of w > $(w), w > (w), v > v(v), 


v > &(v), see Remark 5.4, and the definitions of T and T given in (5.15), we have 
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that 
IO AWA (+(wnlt. jh (00, 0)) +TV- («(w,(t. )); (-0,0)) 
= M [#(wr(t,2*)) — t(w,(t,z7))], 
+ Y tee) (ter), 
T(t) = TV. (#(on(t,-), q); (0, +00)) + TV- (avs. ), a): (0, +90) ) 
= M [polt £t), q) — v(v.(t27),a)], 
cERSn 
+ M [(vn(t,x ),q) — ñ (vnlt, £"), a)l] 
TERSn 
where 


CD, =JreR: (U,(ż,z”),Un(t,x")) is a CD in z < 0 such | 
that wn(t, £7) > maxfwn(t, £t), wą) 


(U, (t, £7), U.(t, x*)) is a RS in z > 0 such 
RS,=<xeER: E 
that v,(t, z*) > max(v,(t, £7), v, } 


Let e, > 0 be the minimal (v, w)-distance between two points in the grid G,, 
namely 
En = min max |v" — v? |w! — w?|}. 
UL,U?EGn 
U!zU? 
The next proposition ensures that the number of discontinuities of U, is uni- 
formly bounded in time. Moreover, it gives uniform bounds on the total variation 


of the approximate solution, which allows us to use Helly's Theorem. 


Proposition 5.7. For any fired n € N sufficiently large and Uon € PC(R;G,), 


we have that: 


(a) the map t > T„(t) is non-increasing and decreases by at least €n any time 


the number of waves increases; 


(b) U,(t,-) € PC(R;G,) for all t > 0. 
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The proof of Proposition 5.7 is deferred to Appendix A.5. 
Beside the bound on the number of wave-fronts proved in Proposition 5.7, 
we need to bound also the number of interactions. This is the aim of the next 


proposition, which together with Proposition 5.7 ensure the global existence of Un. 


Proposition 5.8. For any fixed n € N sufficiently large and Uon € PC(IR; Gn), 
we have that the number of interactions in (0, -+oo) is bounded. In particular U, 
is globally defined. 


'The proof is deferred to Appendix A.6. 

We prove in Proposition 5.7 that U, takes values in G, and we estimate 
TV(U,(t,-)) uniformly in n and t. This together with Proposition 5.8 guaran- 
tee that the number of interactions and the number of the discontinuities of U, 


are both bounded globally in time. 


Convergence 


We first observe that 
lp: — pr] S L(lve = v,| + [we — w,|) 
where L = max{r_, ||1/p'||L<(jp-1(w_),p-1(w,)];R) | because 


p” "(wę Le Ver) if Wer € [w_, w4], 
Per = 
(wer +l—w_)r_ if wer E€ |w_— 1,w_). 


As a consequence TV(p(U)) < L(TV(v) + TV(w)), hence 
TV(U) < (1 + L)(TV(v) + TV(w)). 
Moreover, by Proposition 5.7 and (5.19) we have that for any t > 0 


TV (vn (t, -)) + TV (w,(t,-)) < Ta(é) < T„(0) 


< 
< TV(v) + TV(w) + 2C (tU) i: Ÿ(Uo)). 
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As a consequence TV(U,) is bounded by 
C2027) (TYG) +TV(w) + 2C (T(Uo) + Ÿ(Uo)) ). (5.20) 


Since Uo,n takes values in Q, for any t > 0 we have that also U,(£, -) takes values 


in Q, hence 


[Un (t, -Jlr i) < R+ Vmax- 


Moreover 
Un (t, ) =< VAS, Js e; < Lt = s|, (5.21) 


with Ly = C, max(vmax, Rp'(R)). Indeed, if no interaction occurs for times be- 


tween £ and s, then 


lU. (t. +) — Us ss lireo 
< Y e = 380 (0.5007) - PUC 5,0)) | 


icJ(t) 
+ S5 [e - 8.0 (vnlt, E) — vnlt, 409) 
icJ(t) 


< Lilt sy s|, 


where ôi (t) € R, i € J(t) C N, are the positions of the discontinuities of U, (t, -). 
The case when one or more interactions take place for times between t and s is 
similar, because by the finite speed of propagation of the waves the map t +> U,,(t, -) 
is Lt-continuous across interaction times. 


Thus, by applying Helly’s Theorem, the approximate solutions (U,)n converge 


(up to a subsequence) in LL,(R, x R;Q) to a function 


U € L” (R4; BV(R; Q)) n C° (R4; L,(R; Q)) 


and the limit satisfies the estimates in (5.16). 


Proposition 5.9. Let Up € L''OBV(R; 9) and q € |0,f*] satisfy (H.1) or (H.2). 
If U is a limit of the approximate solutions (Un), constructed in Subsection 5.2.4, 
then U is a solution to constrained Cauchy problem (5.5), (5.6), (5.7) in the sense 
of Definition 5.6. 
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Proof. We consider separately the conditions listed in Definition 5.6. 


(CS.1) 


(CS.2) 


(CS.3) 


Initial condition (5.6) holds by (5.16), (5.21) and the Li 
U, to U. 


io" convergence of 


We prove now (5.9), that is for any test function 9 € C2((0, +oo) x R; R) 


we have 
E U)8,ó + f(U) 8,9) dxdt = 0. 
R+xR 


Choose T > 0 such that $(t, x) = 0 whenever t > T. Since U, is uniformly 


bounded and f is uniformly continuous on bounded sets, it is sufficient to 


prove that 
[ neo Jno + f (Un) 8,9) dxdt — 0. (5.22) 
By the Green-Gauss formula the double integral above can be written as 
jA D (ao) - AjQ(t))o(t.6,(0)at. 
icJ(t) 


where 


Ap(U, (0) = (Un) E, LE) — pn) (950^). 
Af;(t)  F(U.(6,,00*)) — f (Galt, 5,07): 


By construction any discontinuity of the approximate solution U,(t,-) 


satisfies the first Rankine-Hugoniot condition (5.12), therefore 
Snl Ap(U) (t) — AF) 2-0, ie I(t), 


and (5.22) is trivial. 


The proof of (5.10) is analogous because by construction any discontinuity 
of U,(t,-) away from x = 0 satisfies also the second Rankine-Hugoniot 
condition (5.13). 


We prove now (5.11), namely that for any k € [0, Vmax] and test function 
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$ € CX((0, +20) x R; R) such that ¢(-,0) = 0 and ¢ > 0 we have 


JJ, „E U)0,6 + Q*(U)0,6 )dzdt > 0 


where 


0 if v 2 k, 
E'(U) = 
U) PU) —1 ifv<k, 
p-! (W(U) — k) 
0 if v È k, 
Q(U) FU) —k ifv<k. 
p (W(U) — k) 


Choose T > 0 such that (t,x) = 0 whenever t > T. By the a.e. conver- 
gence of U, to U and the uniform continuity of E* and Q^, it is sufficient 


to prove that 


lim inf f g n (E'Uó + Q(U,)2.6)drdt > 0. (5.23) 
0 R 


n— +00 
By the Green-Gauss formula the double integral above can be written as 
[ Y (8.0 (t) AE (£) — AQh(t DEC Si (£))dt, 
icJ(t) 


where 


AER (E) = E'(U.(t,5,()*)) —E*(U,(4,8, (07). 
AQ; (2) = Që (U.(65,(9*)) — a (U.(t5,0)). 


To estimate the above integral we have to distinguish the following cases. 


e If the ith discontinuity is a PT, then we let x = o! (f) and observe that 


p(U, (t, c7) < min(o(U,(1.7*)) p(w- = R), 


DE 23) = Vmax > ma ); 
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hence 


—1 if DE £7) <k < Vmax 


AEĘ:(t) Pr 
0 if k < vnlt, £t); 
D 
| Je f(U.(5,27)) if Un(t, £") <k < Vmax; 
-AQk (t) = Ft 
0 if k < Un (t, De) 


where pk, =p" *(w(U,(t,2")) —k). If vs (fu) < k < Vmax, then 


= o US £7), Un(t, £™)) a a ) SE 02) 


Bre Pr 
F En ) - ) (o(Un(t.27), Unlt,2*)) — e (ok: Ul 2) 
a ee = na n——— 
> 0 


e If the ith discontinuity is a CD, then we let x = 6° (t) and observe that 
ói (t) = v, (t, £7) = v„(t, xt) implies that à? (t) AES*(t) — AQĘ(t) = 0. 


e If the ith discontinuity is a S, then we let z = ô$ (t) and observe that 


p(Us (t, )) < p(U(t,a”)), 
~ c ) > Un(t, x ^ 
f(Us(t a) > f(Us (t7), 
Ši (t) = ec A £7), Un(t,£*)) < 0, 
wa = WA ta = wn (t, £") 2 w, 
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hence 
p(U, (t, TO )) m p(U, t, x )) ifv xt v "m 
P f vnlt, £t) < vnlt, £7) < k, 
AR es p(Un(t,2*)) —1 Uh) < k S vnlt, £7), 
p (ws — k) 
0 if k < vnlt, £t) eo £7), 
FUER IAU) if vnlt, £+) < vnlt, £7) < k 
(wi — k) SA ma | 
-AW (t) = 44 — FUn(t,2*)) if vnlt, £t) < k < ud £7), 
p (ws — k) 
0 if k < wt £") < vnlt, £7). 


If k > v„(t, £7) or k < v, (t, ct), then obviously à? (£t) AE&*(t) — Agt*(t) = 
0. Furthermore, if v(t, z^) < k < v, (t, £7), then 


à, (t) AE" (t) — AQ" (E) 


20 
> 0. 


e If the ith discontinuity is a RS, then we let x = 0% (t) and observe that 
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Wa = w,(t,z )=w,(t,z')>w.-, 
hence 
p(Un(t,2*)) — p(Uu(t, 27) if un(t, £7) < uwadost5oes 
Pw - E) MIU AP AS 
AER (t) = p(Un(t,27)) il gt) < k S ost, 
p^ (ws — k) 
0 if k < au Ju m) 
UCZE x) B f(Us(t z )) if un (t, £7) < Wo £t) < k 
p (wa x TVs d | 
-AQ;(t) = 4 We) |, if v, (ta) < k < v (tx), 
p- (wi cz k) 
0 Ru £7) < RET je 
If k > v (t, £t) or k € v (t, x^), then obviously 
b(t) AES (t) — Agk*(t) = 0. 
Furthermore, if v,(t, x^) < k € v. (t, z^), then 
0, (0 AE (t) — Aq" (t) 
s R p(Un(t, £7)) | T(U, (5 27) 
no E Gz D ) rain | 
" EE )) -1) 
p^ (wę — k) 
>0 
x CUS) U,(t,a+)) + o (Ust. 27), (p^ (wa. — k), 1)) 
| —————————— K, 
<0 


> — pue) (o^ (wo) [p(Uxlt,7)) — (Us (t, 2*))] 


because p(U,(t,2~)) > p^ (ws — k) > p(U,(t,z*)) > r_ and because by 
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the concavity of Li, we have 


0 > o (Un(t, £7), Un(t, £*)) > o (U,(t, a7), (p^ (w.. — k), k)) 


dcl 

> dre (p(Un(t, s 

= wa — p(p(Un(t,27))) — p(U,(t, x ))p' (o(Un(t, x ))) 
d£ —1 —1 If —1 


The above case by case study shows that 


n— +00 


lont | : | (EN(U,.)0.6 + 0*(U„)0.6)drdt 
0 R 
T 


= liminf n 2 (E EAER E) — A E) Jol, 6, (0))at 
nT JO i Rs. (t) 
2 


>, == max / 
an = mo (o) 


xlimint f 25 (P (Unt, E) — p(Un(t, d4(4)*)) Jolt 95 (0)at 


n— +00 


2T : 
> — —||ollqsC max |pp'(p)| = —M, 
r pelr-,R] 


where o? (t) € R, i € RS„(t) C N, are the positions of the RSs of U,(t, -) 
and C, is defined in (5.20). 


We claim that for any fixed h > 0, there exists a dense set K, of values 
of k in [0, Vmax} such that 


T 
iar M 


iERS, (t 


Sr 


(5, (mt (E) - NADU 5 (Jdt > — 


To prove it we fix a,b € [0, Vmax| with a < b and show that there exists 


k € (a,b) such that the above estimate is satisfied. Let 


l= E + 2 
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and introduce the set 


. 2N+1 
zza 


Kn N (a, b). 


Let €, > 0 be the maximal (v,w)-distance between two “consecutive” 
points in the grid G, having the same w-coordinate, namely, with a slight 
abuse of notations, we let 
Eu. mex (vt! — vf). 
(vt, w), (vitt w)EGn 
viz vtl 

Let n, € N be sufficiently large so that £y, < 2/l. Take n > np. We claim 
that for any i € RS„(t) we have 


Ka n (vnlt, OT), t ROH) 


has at most one element. Indeed, if K, has more than one element then 


for any i € RS,(t) we have 


^ 


Un(t, 65 (£)*) — v(t, nlt) ) < En 


As a consequence the sum 


>. (EAER — Ad‘) 


kEK y 


has at most one nonzero element; moreover 


-m(pW)(t, EOT) - o(Un)(t, da(t)*)) < D> (Kg: (0 - Ad"), 
keKy 
=? max lo = aa — 
MT elena | TC, élire 
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Therefore we find 


$5 Y (RAE | - Aq) > -mo,. 


iERS, (t) kc Kj, 


By exchanging the sums, multiplying by the test function and integrating 


in time we get 
` a > ( (öc AES (t) — AQ; (0))e(t,9(0)t >=M. 
kEKh 1ERSn 


Moreover, by construction we have that Kj; is a non-empty set with a 


finite number of elements (it has at most ^ M elements), hence 
h M max [ A (öt AEF*(t) — ARR E) JACE, 6, (£))at >=M. 


In conclusion we proved that there exists k € Kj, C (a,b) such that the 
above estimate is satisfied for any n > np; therefore, since K, has a finite 


number of elements, we have 


mc 


mir f al (c AEF (t (1) — AQ$'(t) ) A(t, 6; (£))at > — 


n— +00 
iERS, (t) 


Since a and b are arbitrary, the above estimate holds true for a dense set 


of values of k in [0, Vmax]. 


Actually, the above estimate holds for any k in [0, Vmax] because the term 
in brackets in the above formula is continuous with respect to k. Finally, 


for the arbitrariness of h, we have that 


lim inf | d 9 (5, ©AEĘ (e) — AQ$(t))ó(t,6, (£))dt >0 


n— +00 
iERSn(t) 


and this concludes the proof of (5.23). 
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(CS.4) We prove now that (5.7) holds for a.e. t > 0, namely 


f(U(t,07)) <q for a.e. t > 0. 


By construction f(U,(t,0^)) < q for any t > 0, namely the approximate 


solutions satisfy (5.7). Since weak convergence preserves pointwise in- 


equalities, it is sufficient to prove that f(U,(t,0*)) weakly converges to 


f (Utt, 63): If o is a smooth test function of time with compact support 
in (0, +oo) and g is a smooth test function of space with compact support 
and such that (0) = 1, then 


À f (U.(t,0-))e(t)dt 


E^ 


The right-hand side passes to the limit, yielding the analogous expression 


] (ENa) + (0. 2) ele) aoa 


oo 


with U, replaced by U. By using again the Green-Gauss formula, one 
finally finds that 


im | f(U.(,0)é(0at = f f(U(t,07))ó(t)at. 
noo JR, R. 
As a consequence f (U.(t, 07)) weakly converges to f(U(t, 07)) and there- 
fore f(U(t, 07)) <q for a.e. t > 0. At last, since we already proved that 
U satisfies the first Rankine-Hugoniot condition, we have f(U (t, 07)) = 
f(U(t,0*)), hence f(U(t,07)) < q for ae. t>0. 


The density flow through x — 0 


Let U be the solution of constrained Cauchy problem (5.5), (5.6), (5.7) constructed 
in the previous section. By Propositions 5.6 and 5.9 we have that non-classical 
shocks of U can occur only at the constraint location x = 0, and in this case 
the (density) flow at z = 0 does not exceed the maximal flow q allowed by the 
constraint. 


In the case of a constrained Riemann problem (5.5), (5.7), (5.18), we know 
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that U coincides with (t, x) > CRS%|[U,,U,|(x/t), moreover if (Us, U+) € Dz then 
the flow of the non-classical shock of U coincides with q. In the next proposition 
we show that also for a general constrained Cauchy problem the flow of the non- 
classical shocks of U coincides with q if the traces at a = 0 of the approximate 


solutions (U,), satisfy a technical condition. 


Proposition 5.10. Let U, € L! n BV(R; Q), F € [0,f*] satisfy (H.1) or (H.2) 
and U be a limit of the approximate solutions (Un)n. Assume that the traces at x = 
0 of (Un)n and U satisfy (5.17), that is for any k € [0, Vmax] and @ € CX((0, +00) x 
R; R) such that ¢ > 0 


T T 
lim x (Un(t,0-)) a(t, at = | N^ (U(t,07))ó(t, 0)dt, 
n— +00 0 0 
with 
o f= | vere 
NE(U) = a WOSE) NT 
k HY qe. 


If the limit U has a non-classical discontinuity at time to > 0 then 


fU (to, 0*)) = q. 


Proof. We first prove that for any k € [0, Vmax} and ¢ € CX((0, +00) x R; R) such 
that ¢ > 0 we have 


Í (f (EU): + Q*(U)0,6 )dz + Nę(U(t, 07))ó(t, 0)) dt > 0. (5.24) 


We stress that (5.24) differs from (5.11) not only for an extra term involving 
Nabigi(U (t, 0*)), but also because here we do not require that @(-,0) = 0. 

Choose T > 0 such that g(t, x) = 0 whenever t > T. By (5.17), the a.e. conver- 
gence of U, to U and the uniform continuity of E* and Q*, it is sufficient to prove 
that 


lim inf i (J (EX(Un)axe + a (Un) Axe) dx + Nk (U„(ż,07))o(t, o) dt > 0. 
(5.25) 
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As already observed in the proof of Proposition 5.9, by the Green-Gauss formula 


the double integral above can be written as 


f 


(EOAR — Aat E Jolt, Edt, 
(t) 


icJ 
where 

AEK (0) = E+ (U(t, Ot) — E* (Ut 8,07). 

AQE (t) = a*(U,(t,5,(¢)*)) — a (Uat, 5,0"). 
To estimate the above integral we can proceed as in the proof of Proposition 5.9, 
with the exception that here the ith discontinuity could also be a NS. In this case, 


ją Var k f 


Figure 5.12: q € (f,f*), vgp = F/p"'(W(Un(t, 0*)) — k) and 
vy = Un(t,0*). The first two pictures show that if vy < k < ug, 
then up p < k. In the last picture we consider the case vg < ug < k 
and show that ug p < Uj p < k. 


that is, if the ith discontinuity is a NS, then 


€ US 07)) eg, s w-0)<me0"): 
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fa Uo =Vmax k 


Figure 5.13: Above q € (0,f-), vy = Un(t,0*) and vg, = 
q/p "'(W(U,(t,0*)) — k). With the first two pictures we show 
that if uy < k < ug, then up, < k. In the last picture we consider 
the case vy < ug < k and show that upp < vj p < k. 


hence 
q B q 
p”'(wn(t,07)—-k) pi (w(u,t. 0+)) — k) 
AQE'(t) = do 05:077) Su te 0t) < k, 
1 -k if v„(t,07) < k < v, (t, 0+), 


if k < vn(t, 07) < (07) 


zd if q £0, 


Ni (U, (t, 07)) 


I 
Tr 
— 
| 
+ 
a 
= 
TS 
S 
ŻE 
TH 
© 
cL 
x 
| 
— 
NE 


ig =. 


Notice that if q = 0, then U,(t,0*) = (0,Vmax) and U,(t,07) € [p !(w.), R] x 
{0}. We observe, see Figure 5.12 and Figure 5.13, that —AQ*'(t) < 0 and that 
—AQh"(t) + NE(U,(t, 07)) > 0 and therefore 


(E) AER E) — Aak*(t)) 9(t,5,(0) + N;(U(Ł,07))ó(t, 0) 
= ETAT + N*(U, (t, 07)}) (0) > 0. 


Thus, by proceeding as in the proof of Proposition 5.9 it is easy to see that (5.25) 
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holds true. Let us just underline that beside the NSs, the only possible stationary 
discontinuities at z = 0 are PTs and CDs, however in both of these cases we have 
f (U, (t, 07)) = 0 and therefore NF(U, (t,07)) = 0. 

We can now prove that if U has a non-classical discontinuity then f (U (t, 0*)) =q. 
This is of course obvious if q = 0, due to (CS.4) and the fact that f(U) > 0. 
We can therefore assume that q > 0 and that z > U(to,x) has a (stationary) 
non-classical shock (U;, U,), with v, < v, and f(U;) = f(U,) = f <q. We want to 
prove that f = q. Consider the test function 


ao: (cow) (LE 0) 


where ps is a smooth approximation of the Dirac mass centred at 07, 62, namely 


Pe € C? (R; R+), E> 0, supp(ó.) C [0,e], PATES = 1, Qe + o. 


Observe that as € goes to zero 


olto, x) = 0 — 0, 


t—to+2e 
do = | ^ ees o. 


—tote 


dolt, x) = ( / d 209 CARE pitt +2) E 


a|—e 


Xp, (00:9(62) > FOR) (D. 


Then by (5.24) for all k belonging to the interval (6(we, q), 6(v,,¢q)) we have 


q*(U) —Q*(U,) + f E E : | 


5 SA -k) J i E a RC - 5) 5 [s R ) a 


Since f < q, the above estimate implies that f = q and this concludes the proof. 
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We underline that the entropy condition (5.11) “becomes” (5.24) if we do not 
require that the test function 6 satisfy the condition $(:,0) = 0. Even if it is not 
necessary for the proof of Theorem 5.2, we conclude this section by considering in 
(5.10) a test function $ which may not satisfy the condition ¢(-,0) = 0. 


Proposition 5.11. Let U, € Lt n BV(R;Q 


and U be a limit of the approximate solutions 


„qe [0,f*] satisfy (H.1) or (H.2) 
Un)n. If the traces at z = 0 of (Un)n 
and U satisfy for any à € C2((0, +oo) x R; R) 


PK — 


T 


„ia | (0,07) [W(U. (t, 07)) —W(Un(t, 0*))] (t, ojat 
a [i f(U(t.07)) [W(U(t,07)) - W(U(1,0*))]. o(t.0)dt (5.26) 


then U satisfies the following integral condition for any ó € C? ((0, +00) x R; R) 


nanc * U)0:6 + f(U)0.6) W(U) dx dt 
= f(U(£,07)) [w(u 07) = W(U(t.0*))]. 40,0) dp 


Proof. Choose T > 0 such that g(t, z) = 0 whenever t 2 T. By (5.26), since U, is 
uniformly bounded and f is uniformly continuous on bounded sets, it is sufficient 


to prove that 


[ [ee 1):6-+ [(U,)0,6) W(U,) dadt 
- [ f (Un(t,07)) [WU (4,07) - WU, (6,0*))]. é(0)0£— 0. (5.27) 


By the Green-Gauss formula the double integrals above can be written as 


[ Y (B (AY; E) — AQi(t)) 6(t, 5) (0))4t, 


icJ(t 


where 


AY; (t) = p(Un(t, 5, (0)*)) W(Us(t,05(0)*)) — p(Un(t,ó,(t) )) W(Un(t, 6, (0) )); 
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AQ; (t) = F(Un(t,67(t)7)) W(Us(t ó(t)*)) — f (Unt, 5, (0 )) W(Un(t, 95 (t) )). 


If U,(t,-) does not have a non-classical shock at 6/(t), then by the Rankine- 
Hugoniot conditions 
& (t) AYG) — AQi (t) = 0; 


moreover, if 6, (t) = 0 and U,,(t,-) has a stationary discontinuity at z = 0, namely 
a phase transition or a contact discontinuity, then v,(t,0*) = v„(t,0) = 0 and 
therefore sign(v,(t, 0*)) = 0. 

On the other hand, if 0? (t) = 0 and U, (t,-) has a stationary non-classical shock at 
x = 0, then 


à; (t) = 0, f (Un(t,0*)) =q, w(U„(ż,07)) > W(Un(t, 0*)), 
and therefore 


& (AY; — AQ) = -q (W(UG(t,0*)) — W(U,(t.0-)) 


Nz 


As a consequence (5.27) is trivial. 
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Networks 


6.1 Introduction 


This chapter is devoted to our traffic flow models on networks presented in [34,43]. 
A network is a directed graph (Z, 7), where Z is a finite set of unidirectional roads 
and J is a finite set of junctions. 

Here we consider the basic case 7 = {J}, and assume that the node J is placed 
at x = 0 with n incoming roads I; = (—0,0) € T, i € | = {1,...,n}, and m 
outgoing roads I, = (0, +oo) € Z, j E J= {n+1,...,n +m}; the general case is 
analogous. 

The evolution of traffic along the roads is described by the PT model in [43] 
and by the LWR model with moving constraint in [34]. At the junction z = 0 we 
propose accordingly two Riemann solvers. We recall that a Riemann problem at 
a junction is a Cauchy problem with constant initial data on each road. Several 
Riemann solvers at junction are available in the literature, see [23,36,58,59]. For 
instance, in [23] the authors propose a Riemann solver for the case of a node having 
more exiting than entering roads, which is determined by applying the following 


rules: 
I) The percentages of drivers who come from the i-th incoming road and take 


the j-th outgoing road is given by fixed coefficients o;j; € (0, 1) such that 


n+m 


y» Qj = 1. 


j=n+l 
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II) Drivers maximize the flow through the junction. 


We underline that rules I) and II) are not sufficient for nodes having more entering 
than exiting arcs. A way out of this problem is proposed in [36,49] for telecom- 
munication networks: the authors design a Riemann solver that inverts the order 
of rules I) and II). 

Differently from the previous chapters, for notational simplicity, each Riemann 
solver associates to any Riemann initial condition the traces at the junction of the 


corresponding weak solutions rather than the weak solution itself. 


6.2 LWR with moving bottleneck on networks 


In this section we consider the LWR model coupled with moving bottleneck [38] 
in the case of networks. A moving bottleneck models a slow vehicle, e.g a bus or a 
truck, which reduces the road capacity at its position. The road traffic with a slow- 
moving vehicle is described by the strongly coupled PDE-ODE system introduced 
in [38,52]. More precisely, the PDE is a scalar conservation law (6.1) which models 
the evolution of traffic, while ODE (6.2) describes the trajectory of the slow- 
moving vehicle. For more informations on different approaches of coupled PDE- 
ODE systems, we refer the reader to [16, 37-39, 42, 64]. 


6.2.1 A single unidirectional road 
We take the Greenshields’ velocity (2.2) with R = 1, namely 
U(p) = Vmax(1 — p). 


Notice that perit = 1/2. Let V, € [0, Vmax) be the maximal speed of the bus and 
y — y(t) € R be its position. Then, y € R is given by the function 


w(p) = min{v(p), Ve. 


The LWR model with moving bottleneck takes the form 


Op + O,f(p) — 0, t>0,xeR, (6.1) 
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g(t) = w(p(t, y(t)*)), t>0, (6.2) 
f (p(t, y(t))) — y(t)p(t, y(t) < ——(Vmax — 90), t20. (6.3) 


AVmax 


Clearly, the trajectory of the bus is determined by the function y : IR, > R 
satisfying (6.2). We underline that if v(p(t,y(t)*)) < V, namely the road is 
highly congested at the bus position, then its speed coincides with the speed of 


surrounding cars. 


Remark 6.1. Condition (6.3) can be obtained by considering the bus reference 
frame, where the velocity of the bus is equal to zero. Indeed, by the change of 


coordinates X = x — y(t), conservation law (6.1) becomes 
op F Ox F(p, y) = 0, 


where the flux function F is defined by 


F(p,y) = Fle) — UP, 


see Figure 6.1. Fix a € (0,1) and consider at X = 0 the reduced flux function 


Ja(p) = ZĘ = 2). 


a 


Notice that a is the reduction rate of the road capacity at the bus position X = 0, 
indeed max, fa = a: max, f. We correspondingly consider at X = 0 the reduced 


flux function 
Falp, Y) = fale) — YP- 


Since 
2 
Fa ; J) = —— max — , 
max Falp, b) = zy — (Vmax — 9) 
it is reasonable to impose that at X = 0 
i A2 
F(p,y) € T (max = ÿ) , 
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namely (6.3). We point out that 


peł, (Vmax m V,) , 


.\2 
dus ex — 9) - ÁVmax 


AVinax 
Let us consider initial datum with jump at z — 0, which is also assumed to be 


the initial point of the bus trajectory, namely 


pt ifm 0, 
p(0,x) = y(0) — 0. (6.4) 
ps Abo, 


We look for self-similar solutions of problem (6.1), (6.2), (6.3), (6.4), thus we let 
y(t) to be piecewise-constant. The standard Riemann solver RSrwa for the LWR, 
model introduced in Section 2.2, does not necessarily satisfy constraint condition 


(6.3). For this reason in [38] the authors introduced the constrained Riemann 
solver BRS wr: [0,1]? > LL; [0, 1]). 


4 
f(p)-PV»o 
Fa 
0 p 
0 Pa 6 Da 1 p PV; 


Figure 6.1: Fundamental diagram with constraint. Left: Fixed 
reference frame. Right: Bus reference frame. 


Definition 6.1. Fix F = 
BRSywn: [0,1]? + LL; [0, 1]) is defined as follows: 


1. If f ((RStwrlee, Prl(W)) < W RSuwelpe; P](Ve) +F, then 


2 . . 
—— (Vmax — Vi) . The constrained Riemann solver 
AVmax 


BRS wrlpe, Pr] = RSiwelpe, prl and y(t) = w(p,)t. 
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2. If f((RStwn[pe PV) > W RSiwelee, pr](V) +F, then 


RSıwrlPe, Pal(é) if E < Vo, 


Rs. bene, SR eM 


BRS twrlpe, Pr] (€) = 
where 


a = min{p € [0,1]: f(p) = pV +F}, 
Pa = max{p € [0,1]: f(p) = pV, + FH. 


We point out that if constraint condition (6.3) is not satisfied for classical Riemann 
solver RS wr, then the solution has a non-classical shock (Pa, a) moving with 


speed of propagation equals Vj. 


6.2.2 Networks 


Let assign to each road of the junction Ip, for h € H=|IUJ, the corresponding 
flux function f,(p) = v^..p(1— p). The corresponding Riemann problem at the 


junction J is then 


0 Oz e d$ xe. 
s fion) di ee hEH, (6.5) 


p(O, x) — Ph LE In, 


where p? € [0, 1] is the constant initial density along road I,. If we further consider 
a bus initially at z = 0 and moving towards road 14, k € J, then beside (6.5) we 


consider also 


Me = w(pe(t, y(t)*)), t» 0, 
y(0) = 0, (6.6) 


fi (ost, y(t))) — (pt. y(t) < o (ax — D), > 0. 


max 


Above, y(t) € 14 is the position of the bus at time t > 0 and o; € [0,1] gives the 
reduction of the road capacity due to the presence of bus. 
Before we introduce in Definition 6.5 our constrained Riemann solver B RS; yg 


for (6.5), (6.6), we give the following definition of admissible constrained Riemann 
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solvers. We denote by RSfwp the classical Riemann solver for the LWR model 
corresponding to the flux fa, h € H. Moreover, let BRStwr be the constrained 


Riemann solver defined in Definition 6.1 and corresponding to road 14. 


Definition 6.2. An admissible constrained Riemann solver at the junction J € J 
for (6.5), (6.6) is a map BRS” : [0, 1]7*7 + [0,1]"*™ such that for any initial da- 
tum (p, ..., Pim) € [0,1]7*% we have that (f1,...,Pn+m) = BRS pl, ..., p... 


satisfies the following conditions: 
(A.1) For every i € |, RSiwr[p?, pi] consists of waves propagating in x < 0. 


A.2) For every j € J\ {k}, RS2 yp [p;, p] and BRSpwn|px, pt] consists of waves 
LWRIPj: Pj k 
propagating in x > 0. 


(A.3) The number of vehicles is conserved at the junction, namely 


(A.4) BRS’ is consistent, that is 


BRS" |p, wę » Pen] = (P1, ot , Pam): 


By the last condition of the above definition, we can observe that the n-tuple traces 
at the junction is a fixed point for BRSi wn. 

Since we wish to maximize the flow through the junction, we propose the pos- 
sible traces and their ranges corresponding to a given initial datum (p4,..., Pm). 


Let us first define the following function. 
Definition 6.3. For any h € H, the function Ty: [0,1] — [0, 1] is defined as follows: 
e fu(m(p)) = (p) for every p € [0,1]; 


e r„(p) £ p for every p € [0,1] \ {1/2}. 


The function 7, is well defined, continuous and satisfies 


On pe 1/2 1/2 < ail) <1, LO e p< 1 4 0 Say <S 1/2. 
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Proposition 6.1. Let p? be the initial datum on the incoming road I; for i € |. 


The set of reachable fluxes f;(p;) for an admissible constrained Riemann solver is 


[0, fi) if p? = [0, 1/2], 


9t (p?) = 
[0, (1/2) ifo € (1/2,1]. 


Proof. We apply the construction done in 47, Proposition 4.3.3]. We consider 
the case p? € [0,1/2]; the case o? € (1/2, 1] is analogous. By (A.1) we have that 
RSiwslp?, pi] must have only waves propagating in z < 0. If p; € {p?} U(r;(p9),1] 
then RSiwp[p?, Pi] is either constant or has a single shock with negative speed. 
On the contrary, if p; € [0, 7:(p?)| V (99) then RStynlp, Pil is either a rarefaction 


or a single shock, however both of them propagate in z € 0, which concludes the 


proof. 


From the above proposition the following conclusion is clear. 


Corollary 6.1. The maximal flow of the incoming road I; at the junction J is 


ymax( 5 - Ue if py € [0, 1/2], 
Mk 1/2) if o? € (1/2,1]. 


We observe that there exists a unique p; € [0,1] so that f,(pi) = *P**(p9) and 
RSiwrlp?, Pi] has waves propagating in x < 0. 


Proposition 6.2. Let j € J and p; be the initial datum on the outgoing road I;. 


The set of reachable fluxes f;(p;) for an admissible constrained Riemann solver is 


In f5(0/2)] à pf € [0, 1/2] and j # k, 
suo) = 4 UL FiCS) ifo) € (1/2,1] and j # k, 
[0, fe(Ba,)] if P € (0, fol, 
th (| df 0$ € (Pa 1]. 


Proof. Except for the case j Z k, the proof is analogous to proof of Proposition 6.1. 
Consider j = k and p? € [0, 94,]. We observe that p; € [0,5„,| can be joined with 
py by a classical wave. If py € (Pa,,Tk(Pa,)) U {Pa,} then BRSpwn|pk, py] has 
a possibly null shock (pk, Pa,), followed by a non-classical shock (0a,, Pa.) and a 
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shock (5a,,pę). We point out that py € [T(fo,), Pa.) cannot be connected with 
fa, by a wave propagating in x < 0. At last, the case j = k and py € (Pap, 1] is 


analogous to the situation when j Æ k and p? € (1/2,1]. 


Corollary 6.2. The maximal flow of the outgoing road I; at the junction J is 
(1/2) if p? € 


I" 

BA 5 =. Fis) if p; = 

: i TS if py € 
fp) if pk E 


— 


0, 1/2] and j źk, 
1/2,1|andjźk, 
0, Pax]; 
Dan» 1]. 


Additionally, for any j # k there exists a unique p; € [0,1] so that f;(p;) = 


max 


W 
exists a unique py € [0,1] so that fy(Px) = V%*(p}) and BRS welpr, p] has 


(p°) and RS! yvrlP5, Ph has waves propagating in x > 0. Analogously there 


waves propagating in x > 0. 


For each junction we define a traffic distribution matrix, namely the matrix 
representing the percentage of cars from incoming road J; moving towards outgoing 


road [;. 


Definition 6.4. A distribution matrix Amxn for the junction J € J is given by 


Qn+1,1 tt Qn+ln 
As z , 
Qn+m,1 ^^ An+tm,n 
n+m 
where o; > 0 for every i,j and D aj; = 1 for every i. 
j=n+l 


More precisely, if C is the number of cars coming from road J;, then Cay, is the 
number of cars moving from J; towards J;. 
We conclude this subsection by defining our admissible constrained Riemann 


solver at the junction J. 


Definition 6.5. Fix a distribution matrix Amxn. The associated constrained Rie- 
mann solver BRSwr : [0, 1] — [0, 1]^*" is defined as follows for any initial 


datum (pf. .., pou) € [0,14 
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1. Consider the closed, convex and non-empty set 
Mat Ain) eue KOO AGE aie EMail: 


where 9t; = Ni(o;) = [0,a*(e})] and 9t, = (pi) = [0,97*(pj)] are 
respectively defined in Propositions 6.1 and 6.2, see also Corollaries 6.1 
and 6.2. 


2. Let (%1,...,n) EN be such that 


25 E max „) 5i (6.7) 


and for any i € | consider the unique p; € [0,1] so that f;(p;) = % and 
RSiwrlp?, pi] has waves propagating in x < 0. 
3. Define 
(abe Ray PEM (as E 
Then for any j € J consider the unique p; € [0,1] so that f;(p;) = "yj and 
RSI we Dy, p\] has waves propagating in x > 0. 


4. At last, define BRStyn|P, --- Pam] = (pu... Pam). 


Remark 6.2. We underline that problem (6.7) may have more than one solution. 
In such a case, we can impose additional requirements on the distribution matrix 


A as in [47, Section 5.1], or introduce the priority vector, see [40]. 


6.2.3 A case study 


We study a junction having two incoming (n = 2) and two outgoing (m = 2) roads 
with corresponding flux functions f,(p) = 4p(l — p), h € {1,...,4}. Consider 


constant initial density (p, ...,p4) € [0, 1]*, see Figure 6.2, center, as follows: 


0 < p < 1/2, 1/2 < <1, 1/2 < p <1, 1/2< ed 
F(p?) = 1/2, f(P2) = 2/5, f(e3)=7/10, =f (ef) = 1/2. 
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Let us choose the parameter a3 appropriately to obtain f(a) = 7/20 and consider 


że K A 
1/2 2/3 


the distrution matrix 


max 
Y 


Figure 6.2: Left: the set 9t. Center: the fundamental diagram 
with initial datum. Right: the set N. 


We consider two cases. The former one describes a junction without bus, while 
the latter describes a bus with maximal velocity V, = 1/6 entering road /3. By 


Propositions 6.1 and 6.2, the sets of admissible fluxes at the junction are 

Nı = Ma = [0,1/2], N= [0,1], N = [0,7/10], N = [0, 7/20]. 
Above, N? and Nz correspond to the cases when the bus is present or not, respec- 
tively. 

In the case without the bus at the junction we have 


N = Uu %2) EM x M: A (41,72) € Ne x NU). 


The maximal admissible flow through junction max(4, „s)en(71 + 72) is achieved at 
the point Q = (1/2,3/8) € N, see Figure 6.2, left, therefore 


The solution of Riemann problem at the junction for (6.5) is then completely 
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determined. 


In the case with the bus, the set of reachable fluxes is given by 
W = [(71,72) € 9 x Ma: A (71, 99)" E W x Mu}. 


The maximal admissible flow through the junction max, „s)ese (7172) is achieved 
at the point Q^ = (2/5,9/20) € NÈ, see Figure 6.2, right, thus 


1/2 <p, « 1, 1/2 < p < T, P3 = fos, Pa = Pa; 
f(p1) = 2/5, J(p2) = 9/20, f(ps) = 7/20, f(p4) = 1/2. 


The solution of Riemann problem at the junction (6.5), (6.6) is completely deter- 
mined. 

In Figure 6.3 we display both solutions at time t — 1/5. The blue line describes 
the solution for the problem (6.5) without the bus. We point out that in this case 
the solution is constant along roads J, and /4, while along each road I> and J; it 
has a single shock. The red line in Figure 6.3 describes the solution of problem 
(6.5), (6.6) in presence of the bus. Notice that in this case the solution has a 
single shock along road 14, a single rarefaction along road /5, a non-classical shock 


followed by a shock along road /3, and it is constant along road T4. 


6.3 PT models on networks 


In this section we adapt the Riemann solvers for PT models described in Section 2.4 
to the case of road networks. More precisely, we introduce a Riemann solver at 
the junction SZ, obtained by adapting the Riemann solver RSpr. Here we 
use RSpr to indicate the Riemann solvers RS%, RSP, RSS and RS} defined in 
Subsection 2.4.3. 


The Riemann problem at J is 


PT model iUe 


(6.8) 
Y,(0, z) = Y? z € I, 


where (Y,..., Y? m) € Q"*" and “PT model” stands for phase transition model 
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Road J; Road I2 


p(t=1/5) p(t=1/5) 


>| 0.887 


iJ 0.146 


, Road I3 Road I4 
plt=1/5) plt=1/5) 
0.903}, 


0.853 
0.773 


> > 


0 x 0 x 


Figure 6.3: Above we represent the densities at time t = 1/5 along 
the roads /,,...,74. In red the case with the bus and in blue the 
case without the bus. 


(2.20). In particular the traffic evolution along each road Ip, h € {1,...,n +m}, 
is described by PT model (2.20). 
As in the previous section, before stating our Riemann solver at the junction 


RSpr for (6.8) we give the definition of admissible Riemann solvers. 


Definition 6.6. An admissible Riemann solver at the junction J € J for (6.8) 
is a map RS”: QF" — Qn. such that for any initial datum (Y9,...,Y9,,) € 


^^ nm 


Q^" we have that (Yy,..., You) = RS|Y9,..., Y? n] satisfies the following 


^-^ nm 


conditions: 
e For every i € | we have that RSpylY/, Y;] consists of waves propagating in x < 0. 


e For every j € J we have that RS Y? consists of waves propagating in 
x > 0. 
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e The number of vehicles is conserved at the junction, namely 


2,f%)= 2, f(y. 


e RS” is consistent, that is 


RSW wad sa ck o8): 


As in the previous section we study the possible traces and their ranges corre- 


sponding to a given initial datum (Y?,..., Y? m). 


Proposition 6.3. Fir initial datum (Y$,...,Y9,,,.) e Qm, 


^^ naim 


e If Or OQ, z 0, then the set of reachable fluxes is given by 


O(Y?) = 


e If Q NQ. — 0, then the set of reachable fluxes is given by 


0, f(Y?)] if Y? ENT and f(Y2) < fr, 
oy) = 4 Oke | if Y? € OF and f(Y?) > fe, 

0, f(vt(Y?))] u C/)] if Y? EQF, 

0, f(v$(Y?))] u £12). i£ Y? € Qs. 


Proof. It suffices to make the observation that waves with negative speed can be 


of three types: 


e waves of the first characteristic family, 


e phase transitions (Y,, Y.) € QF x O7 with f(Y;) > f(Y,) and we < w, = w, 


e phase transitions (Y, Y.) € (Of x QZ) U (Oz x QF) with w, = w, (in the 


case NN = ()). 
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The set O;(Y;/) is non-convex, therefore we do not consider the metastable states 


in the definition of the corresponding maximum flows. Hence, we let 


y if Y? € Qr, 


( if Q4 0, Z 0, 
FWUYP)) ifY? e. 
WYO = À (709) EYL EQ and FY?) «t. 
ts: if Y? EQ; and f(Y9) 2 fj, ENAR. — 9. 


FH) EYP € Qf UN, 


Corollary 6.3. For any initial datum Y? € Q and 3, € O;(Y/), there exists a 
unique Y; € Q such that RSpr[Y!, Y;] consists of waves propagating in x < 0. 


Proposition 6.4. Fix initial datum Y? € Q forj € J. The set of reachable flows 
is given by 
[0, ff] if Y? EM, 


[0, f( (PM) if YP € Me. 


Proof. It suffices to make the observation that waves with positive speed can be 


O;(Y?) = 


of two types: 


e contact discontinuities, 


e phase transitions (Yy, Y,) € Q; x Oz with fe < f, and w(Y;) < w, = w. 


We point out that in Proposition 6.4 we do not have to distinguish between the 
cases of intersecting and non-intersecting phases. The set O;(Y?) is convex and 


we denote its maximum as 
jn if Y? € Qr, 
fi) if YP ER. 


Corollary 6.4. For any initial datum Y? € Q and 7, € Oj(Y?), there exists a 
unique Y; € Q such that R Spr [Y?, Y;] consists of waves propagating in x > 0. 


We introduce now the definition of our admissible Riemann solver at the junc- 
tion J. 
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Definition 6.7. Fix an (n4-m)-tuple of priority parameters (01,...,Onim) € IR"*" 
such that 0, > O for all h € H and 


n+m 


A= Dot 


j=n+1 


The associated admissible Riemann solver RSpq: O+ + (Q"*" is defined as 
follows for any initial datum (Y£, ..., YP m) e (m, 


^^ ntm 


1. Define the maximal possible through-flow at the junction, namely 


n+m 
De min} pa Y Jr op 


j=n+l1 


2. Consider the closed, convex and non-empty sets 


Cine = [Q9 € Floor al nd 
n+m n+m 
Cow = cz ee news) € lI ©. 4 (Y?)] : 9 = B . 


j=n+1 j=n+l1 


3. Denote by (%1,...,7n) the projection of the point (T6,,...,1'0,) on the conver 
set Cine and by (Fn+1,---, ^nm) the projection of (T0541, ..., D054,4) on the 


convex set Cout. 


4. For every h € H, set 


Yno if f(VYao) = Fh, 


Y, otherwise, 


Y= 


where Y, is the unique solution to the equation f(Y) = ^y, given by Corol- 


lary 6.3 or Corollary 6.4. 
5. At last, let RSprÓ D ns Y Bede, Yat: 


We point out that RS}, is admissible in the sense of Definition 6.6. 


Nikodem Dymski 


Section 6.3 Page 120 


6.3.1 A case study 


Consider a simple network made of two incoming roads 4 = 1,2 and one outgoing 
road j = 3 that meet at z = 0. Assume that the traffic on each road is described 
by either PT* or PT?. Once the initial data (Y, Y2, YP) € Q’ is chosen, we can 
determine the maximum flows 4T^*(Y9), 42ex(y3), 4mex(Y92). Then, we define 
I = min(4f** + 47e, 4m9*Y and fix (01,02, 03) = (0,1 — 0,1), for some 6 c (0, 1). 
We have 


Cic = { (1572) € [0, Ne (YP)] x x [0, Ao AMA Yı F c= D Cout = {T}. 


Then, we take 73 =I. We have to distinguish two possible cases: either ^T'?*(Y 9) 4- 

ges y n = T OK or Des Re Le): In the first 
case, we take (74,72) = (yP9*(Y9),yme*(Y2)). In the second case, let (71, 72) 
be the point of intersection of the sets {(71,72) € IR* x Rt: 44 4- 59 = T] and 
(1,92) € R* x R*: y2 = ((1—0)/0)m). If 051,52) € [0,475 (YP) x [0,9277 (Y2)], 
then we take (74,72) = (71,72). Otherwise, (74,72) is the projection of (71,72) on 


Cine. Finally, we can find the traces at the junction (Y;, Y2, Y3) by Corollary 6.3 
and Corollary 6.4. 
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Technical proofs of Subsection 2.4.3 


In this section we collect the proofs concerning the properties of the Riemann 
solvers RSR and RSs. 


A.1 Proofs of the main properties of RSR 


Proposition 2.2 follows from the following two lemmas. 


Lemma A.1. RSR is Li,.-continuous in Q*. 


Proof. Assume that (YF, Y) converges to (Yp, Y.) in Q? as e > 0. We claim that 
RSg[Y?, Y] > RSRIY, Y,] in Li... We consider only the case where RSg[Y;, Y,] 
is just one PT, namely the case described in (R.4) of Definition 2.7; for the remain- 
ing cases it is sufficient to apply the continuity of RSpwr and RS Anz. Observe 
that Sg |Y7, Y] is the juxtaposition of RSg[Y7, Y3 (Y?)| and RSg|v; (Y), Y]. 


Therefore we are have to consider only the following cases. 


e Assume that p; = 0. If v4 (Y?) > v4 (Y), then of YF, v4 (Y7)) > v(Y,). Then 
we can conclude by observing that u(Y£) + v(Y,). 

e Assume that p, Z 0 and Y, z Y,(Y;), with Y, defined in (R.5). In this case 
w, = w, therefore v; (Y?) — Y,, which implies that o (YF, v; (Y*)) > o(Y,,Y,). 
Then we can conclude by observing that v(Y/) + v(Y,). 

e Assume that p, Z 0 and Y, = Y,(Y;). In this case w, = w_ and RSR[Y/, Y] is 
described by either (R.4) or (R.5). The former case is analogous to the previous 
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one. For the latter case, by Y (YF) + Y, and v; (Y?) + Y, we have that both 
c(Y?, uj(Y£)) and APT (v5 (Y?)) converge to o(Y4, Y.) = APT (Y,), then it is easy 


to conclude. 


Lemma A.2. The Riemann solver RSR is consistent. 


Proof. Since RSn|z, z] = z for all z € Q, we can assume that Y;, Y and Y, are 
all distinct. The Lax Riemann solver is consistent, hence it is enough to consider 
RSR|Y, Y,] with at least one PT. In this case by definition RSR has exactly one PT 
and it goes from Q; to Q7. Moreover, &Sg|Y;, Ym] does not contain any C. Hence, 
to prove (I) or (II) of Definition 2.9 it is sufficient to consider the cases described 
in (R.4) or (R.5), p, # 0 and wm = w_. The result then easily follows. 


A.2 Proofs of the main properties of RSs 
Proposition 2.3 follows from the following two lemmas. 


Lemma A.3. The Riemann solver RSs is LL, continuous. 


Proof. By the analysis done in Lemma A.1 and Remark 2.5, it is sufficient to 
consider the case in which &Sg|Y,Y,] = PT[Y;, Y,] and one of the conditions 
(2.22),(2.24),(2.23) holds true. Let (Y£,Y?) + (¥,Y,) in Q? and consider the 


following cases. 


e Assume that (Y; Y,) satisfies (2.22) with Y, = yf (Yp) and v(Y;) = ve. In this 
case it is sufficient to observe that v£(Y7) — Y; and vi(Y£?) + Y,, thus both 
AP (V7) and AP" (s(V7)) converge to APT(Y;) and a(i (YF), WO) > o(Y Y). 


e Assume that (Y;,Y,) satisfies (2.23) with Y, = Y°. In this case, it is sufficient 
to observe that v; (Y?) > Y., AP! (ug (Y7)) + M (Y) and o(¥F, Y.) > o (Ye, Y.). 


e Assume that (Y;, Y,) satisfies (2.24) with Y; = y£ (Y,) and v(Y,) = ve. Let Y? be 


the intersection point of Lu and Ls In this case it is sufficient to observe that 


w°(Yf) + Y, and Y? + Y, to obtain that both APT(u£(Y7)) and APT(YE) converge 
to AT (Y), and e(Y7, v1(Y7)) > o (Ya, Y.). 


Lemma A.4. The Rieman solver RSS is consistent. 
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Proof. Since RSg|z, z| = z for all z € Q, it is not restrictive to assume that Y, 
Ym and Y; are all distinct. Furthermore, RSg|Yę, Ym] does not contain any C. By 
Lemma A.2 and Remark 2.5, to prove (I) or (II) of Definition 2.9 we are left to 


consider the following cases: 
e (Y;,Y,) satisfies (2.22) with Ym such that wm = we and Um > vy; 
e (Y,,Y,) satisfies (2.23) with Ym such that wm = w.. and um > w; 


e (Y, Y,) satisfies (2.24) with Ym such that wm = w, and Um > v. 


The analysis of the above cases is straightforward, and the result easily follows. 


A.3 Proofs of the main properties of CRSs 


In this final section, we accomplish the proof of Proposition 5.3. Recall that the 


same constrained Riemann solver for PT? has already been studied in [12]. 


1 
loc 


Indeed, take F > f(uż) and consider U,,U,,UĘ € Qs with f(U,) = q < f(UF) 
and U; — U,. In this case CRSs|U$,U,| does not converge to CRSs|U,, U,] in 
Lipe, since CRSs|U,, U,| = U, in R_ and the restriction of CRSs|U$,U,] to R 


converges to 


Example A.1. The constrained Riemann solver CRSg is not Li .-continuous. 


1 
loc’ 
Ur ax 9 (Us, 24), 
Ze tf o(Ue, 2) < x < 0, 
where zy = u* if U; EQ; and zz = v£(U;) if U, € OF. 


Lemma A.5. The constrained Riemann solver CRSs satisfies (II) of Defini- 
tion 2.9. 


Proof. Assume that CRSg|U;, U,,](z) = Um = CRSs|Um, U,.](x), for some x € R. 
Since by Lemma A.4 and Lemma A.7 the Riemann solvers RSs and CRSg al- 
ready satisfy (II), by Remark 5.2 we have to consider the cases where at least one 
among (Ur, U,), (Uc, Um) and (Um,U,) belong to D» and satisfy one of the condi- 
tions (5.2),(5.3),(5.4). We observe that C(RSs|U;, Um] cannot present any contact 
discontinuity, otherwise it would not be possible to juxtapose CS s[U;, Um] and 
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CRSs|Um, Ur]. Hence, we are left to consider (Up, U,) € D» satisfying (5.2) with 
Um € {U, Ù}, or (5.3) with Um € {U, Ù}, or (5.4) with Um € {Ü}UCRSs[U, U](R_). 


In all these cases, it is easy to see that (II) holds true. 


A.4 Proofs of the main properties of CRSp 


In the next lemmas we prove Proposition 5.2. 


1 


loc Continuous. 


Lemma A.6. The constrained Riemann solver CRSp is L 


Proof. By the Li, -continuity of RSR, it suffices to consider (U$, UE) > (Us, U») 
with (UF, U£) € Dz and to prove that RSR[US, 0] — RSp(U„U,] in R- and 
RSRIU*, UE] > RSn|[U;, U,.] in Ry, where 0* = U(w(U$),F) and U* = U(v(U$),F). 
For notational simplicity, below we denote U = U(w(U;), F), U = U(v(U,), F) and 
ty i (Ug, Up). 


We first consider the cases with (U;, U,) € Di. 


e Assume U, U, € Qf and f(U,) = q. Then, either U, € Q; or U, € Of. In 
the first case, o(U$,U*) — 0 and RSp(U$,U*) — U in R_, while in the latter 
Ü* — U, and RSR[US, 0]  RSR[U;, Ue] = Ue. Moreover, in both cases Ù — U; 
and RSg[U5, UE] — RSg(U,, U,]. 

e Assume U;,U, € Qe and f(u,) = q. Then, Ue > u, and RSg[|U;,0*] — 
RSg|[U;, u.]. Moreover, it is sufficient to consider the cases Ù € Q; and UE € Q. 
In the first case Wz (US) + us, o (UF, dy (U£)) + 0 and RSg[U5, U] > RSn[u., U,] 
in R}, while in the latter U* > u, and RSg[U5, UE] — RSqg[u,, U,]. 


e Assume U, € Q7, U, € NF and f(v4(U;)) = q. Then, 0* — v4(U;) and UF = 
V1(U,). As a consequence, RSR[US, 0] 4 RSR[U;, Y(U] and RSg[U5, UE] = 
RSr[Y: (Ue), Ur]. 


e Assume U, € O7, U, € OF and f(U;) > f(y (U,)) = q. Then, Ü = v; (U,) 
and therefore RSR[U;, 0°] — RSR[U;, 43 (U,)]. Moreover, o(U*, #3 (UZ)) — 0 and 
therefore RSę[U*, UF] > RSn[v; (U,), U,] in Ry. 


e Assume U, € Qf, U, € OF and f(44(U,)) > f(U) = q. In this case, 


o(Uę,Ó*) > 0 and RSp[U$,6*] + U, in R_. Moreover, Ü° = U; and RSR[Ü, UE] + 
RSqg|U,, U,]. 
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We finally consider the cases with (U;,U,) € Do. 


e Assume U;, U, € Qr and f(U;) > q. Then, Ô + Ü and Ü° > U. Asa consequence, 
RSqn[U;, Ù]  RSR(U, 0] and RSR[Ü°, UE] + RSR(U,U,.]. 


e Assume U;,U, € (X and f(u.(U;,,U,)) > q. Then, U + U and therefore 
RSR|U$,U*] > RSn|[U;, U]. Moreover, it is sufficient to consider the cases U* € Qe 
and Ü € Or. In the first case U^ > U and RSga[U,U:] + RSn[U, U,], while 
in the latter (whether U = v4 (U,) or U = U*) we have yz (U£) — 44 (U.) and 
RSg[U5, UE] + RSn[U, U,] in Ry. 


e Assume U, € QJ, U, € OF and f(V4(U;)) >q. Then, Ue + Ü and Ue =U. Asa 
consequence, we have ASg[U;, U*] 4 RSn[U;, 0] and RSn[U5, UE] — RSRIU,U,.]. 

e Assume U, € QF, U, € OF and min(f(Uj), f(v3 (U,))) > q. Then, 0* = 0 
and U* = U. As a consequence, RSp(U;,U*] — RSR[Ur, 0| and RSg[U5,Uz] ^ 
RSalD U,]. 


This concludes the proof. 


Lemma A.7. The constrained Riemann solver CRSR satisfies condition (II) of 
Definition 2.9. 


Proof. We assume 


CRS p[Ue, Un](£) = Um = CRSn[Ua, U;](2), for £ER. (Zl) 


Since by Lemma A.2 the Riemann solver RSR satisfies (II), it is not restrictive to 
assume that 
{(U, Us): (U, Um), (Um; U,)} 1 D2 X 0. (Z2) 


We also observe that (Z1) implies 
CRS n|U,, Un] does not contain any contact discontinuity, (Z3) 


because otherwise it is not possible to juxtapose CRS g|U;, Um| and CRS n|U,,, Ur]. 


We are then left to consider the following cases. 
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e Assume U;, Um € Qg. In this case, by (Z3) we have (Up, Um) € Dz and U = 
U(v(Um), F), namely f(U,) > q = f(Um). Then, by (Z1) we have that either 
Um € Q; and f(%; (U,)) > q or Un € QF and U, € Oy. In both cases it is easy to 


conclude. 


e Assume U,,U,, € Qe. In this case, by (Z3) we have either (Ue, Um) € Do or 
(Uc, Um) € Dı and w(U;) = wm. In the first case, whether U(v(Um),F) € Qe, 
Um = U(v(Um), F) and Wm < w(U;) or ((v(U4), F) € OF, Wm = w- € w(U;) and 
f(Um) > q, by (Z1) we have that v(U,) = v(Um). In the latter case, by (Z1) and 
(Z2) we have that f(U,,) = q, v(U,) > v(U,,) and (Um, U+) € Dy. In both cases it 
is easy to conclude. 

e Assume U, € Qz and Um € 0; . In this case, by (Z3) we have (Ur, Um) € D» and 
Um = U(u(Um), F). Then, by (Z1) we have that U, has to satisfy f(v; (U,)) >q. 
Hence, it is easy to conclude. 

e Assume U, € Q; and Um € Q7. In this case, by (Z3) we have w, = w. 
Moreover, by (Z2) and (Z3) we have either f(Um) =q < f(U>) and v(U,) > v(U,,) 
or f(Um) > q and v(U,) = v(Um). In both cases it is easy to conclude. 


This concludes the proof. 


Finally, we accomplish the proof of Proposition 5.4 on the minimal invariant 
domains for CRSR. We remark that (q/V, Q(g/V)) is the point of intersection 
between the lines {U € Q : f(U) = q} and Qg: if q > Vo_, this point belongs to 


the region Qe, otherwise it is in Qf. 


Proof of Proposition 5.4. Let us first prove (ICR.1). The invariance of 7; is an 
easy consequence of Definition 5.7, hence we are left to prove the minimality of Zr. 


Let Z be an invariant domain for CS g containing Qr. Then, Z has to contain 


CRSpl((U U.) € À : F(U) > q))(R) = AU {U E€ Re: F(U) = q} UD. 


As a consequence, Z has to contain also 


CRSp[((ULU,) € À : f(Ue) = f(U-) = q, v(Ue) > v(U,)]](R) 
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Ti fq Vox, 
{U EL: JO) Za} ate Vas 


Finally, if q > Vo-_, then Z has to contain also 


CRS pl{(Ue, Un) € OF x Ne: f(U) < a = F(U,)HR) = (U € : f(Y1(U)) € a}. 


In conclusion we proved that Z 2 Ty. 

Now, to prove (ICR.2) it is sufficient to observe that by Definition 5.7 we have 
that there exist U;, U, € Qe such that the values attained by CRS;ę[U;,U,.] exit Qe 
if and only if q < Vo_, and in this case 


CRSg[Q2)(R) \ Q = CRS R[{(Ue, Ur) € 02 : fog (U,)) > GR) \ Re 
= ((g/V, Q(a/V))) C M. 


This concludes the proof. 


A.5 Proof of Proposition 5.7 


By construction, for any t > 0 sufficiently small U,(t,-) belongs to the space 


PC(R; Gn), more precisely it is piecewise constant with jumps along a finite number 


of straight lines. If at time t > 0 an interaction occurs, namely two waves meet 


or a wave reaches z = 0, then the involved waves may change speed or strength, 


while new waves may be created. To prove that U,(t, :) belongs to PC(R; Gn) we 
have to provide an a priori upper bound for the number of waves, which follows 
from (a). 

Clearly, if at time t > 0 no interaction occurs then m(t) = 7,(t+). For this 
reason we consider below all the possible interactions and distinguish the following 


main cases: 
e a single wave reaches z = 0 and no NS is involved; 
e a single wave reaches x = 0 and a NS is involved; 


e two waves interact away from z = 0; 
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e two waves interact at z = 0 and no NS is involved; 
e two waves interact at z = 0 and a NS is involved. 
For completeness we estimate 


ATV, = TV (v.(t4,-)) — TV(wn(t,)), AT = Ta(wn(ts,-)) — Ta(walt,-)), 
ATV, = TV(wa(t+,:)) — TV(wn(t,-)), AT, = Ta(valts,-)) ^ Ta(vn(t,-)), 


and 
Al = t(t+) — HE), AT, = Fults) - Talt). 


For simplicity in the exposition, whenever a NS is involved we consider sepa- 
rately the cases q € [f~,f*) and q € [0,f ^). Notice that w” > w, if and only if 
q < f^, or equivalently Vmax Æ v$. Notice also that if q = f* then Dz = Ø, while 
if q = 0 then D, = Ÿ. At last, notice that if q € [f ^, f^) and (U,,U,) € Do, then 
We = w, and V, = vr. 

In Table A.1 we collect the most relevant possible interactions considered in 


this section and list the corresponding possible results in terms of wave types, Af 


and AJr. 


We start with the interaction estimates. 


e If a wave (U;,U,) reaches x = 0, Un(t,07) = Un(t,07) and (U,,U,) € Di, then 
the constraint has no influence on the wave and 0 = ATV, = ATV, = At. Since 
any CD has non-negative speed, we have that AT, < 0. Since any RS has negative 
speed, we have that AT, < 0. As a consequence A7, < 0. 


e Assume that a wave (Us, U,) reaches x = 0, U,(t,07) = U,(t, 0*) and (U„U,) € 
Dy. 


If q € [f-, f^), then one of the following cases occurs: 


CDj (U2, U,) is a CD. In this case Ÿ, > v; = v, = V, > We, We = We > Ùr 2 Ür > Wy 
and f(U;) > q > f(U,). CRS%|U,,U,] has at most three waves (Up, Up), (Us, Ù) 
and (U,,U,) that are a S, a NS and a possibly null CD, respectively. As a 
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Interaction | Result At AT, 
CD} (S,NS,CD), (S,NS) € {1,2} <0 
RS; (NS,CD) =1 <0 
CD, | (S,NS,CD), (SNS), (PT,NS,CD), (PT,NS) efl2 |<0 
RS, (NS,PT), (NS,CD) =]: <0 
CD-S | (S,CD), (PT,CD), PT <0 =0 
CD-RS (RS,CD) =0 =0 
CD-PT | (PT,CD), PT, (CD), 8 <0 <0 
S-S S <0 =0 
S-RS S <0 <0 
RS-S S <0 <0 
PT-S PT <0 =0 
PT-RS PT, CD <0 <0 
CD-S, | (S,CD), PT,CD), PT <0 <0 
CD-RS, | (RS,CD) =0 <0 
CD-NS, | (RSs,CD), RSs e [-1,2"-1] | «0 
OD-PT, | (LCD), PT, (S,CD), 8 <0 <0 
S-So S < 0 = 0 
S-RSo S <0 < 0 
RS-So S «0 «0 
NS-So (S,CD), CD «0 «0 
N&PT, | (S,CD), 5, CD <0 <0 
PT-So PT <0 =0 
PT-RSo PT, CD < 0 < 0 
PT-NS CD, PT < 0 < 0 
CD-S7 (S,NS,CD =. <0 
CD-RS; | (5,NS,CD), (NS,CD), (S,NS), NS € {-1,0,1} | «0 
CD-NS; | (RSs,NS), (S,NS) € [0,27 -2] | < 
CD-PTĘ (S,NS,CD = 1 < 
NS-S* (NS,CD) =0 <0 
NS-RS; | (NS,CD) =0 <0 
CD-S7 (S,NS,CD =] < 
CD-RS, | (S,NS,CD), (S,NS,PT), (NS,CD), (NS,PT), (SNS), NS | €{-—1,0,1} | < 
CD-NS, | (RSs,NS), (S,NS) e [0,27 —2] | «0 
CD-PT, | (S,NS,CD), (SNS) exor. 7) 40 
NS-S7 (NS,CD) =0 < 
NS-RS, | (NS,PT), (NS,CD) =0 < 
NS-PT; | (NS,CD) =0 <0 


Table A.1: Overview of the interactions considered in the proof of 


Proposition 5.7. 


consequence 


ATV, = 2(ve = 14) > 0, 


AV, 
AT. 
AŤ, 


therefore Aj € {1,2} and AT, < —2(¥, — ve) < 0. 
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RS; (U„U,) is a RS. In this case ve = v, < Up = ©, We < We = We = Wr, 
f(U) =q < f(U,) and Uz, U, € Qe. CRS%[U;,U,] has two waves (U;,U,) and 


(U,, U,) that are a NS and a CD, respectively. As a consequence 


ATV, =0, 
ATV, = 2(w =) > 0, 

AT, — 0, 

AT,--[W,.-v4].-][& 9]. = — (v, — ve) = (w = 9.) < 0, 


therefore Aj = 1 and A7, = —2(v, — vj) < 0. 


If q € |0, f7), then one of the following cases occurs: 
CD, (Us, U,)is a CD. In this case $, > v; = v, = V, > Ve, We > We > Ù, 2 8,2 wu, 
and f(U;) > q > f(U,). CRS%|U,,U,] has at most three waves (Ug, Up), (Ue, Ù.) 
and (U,,U,) that are a S or a PT, a NS and a possibly null CD, respectively. 


Asa consequence 


ATV, = 2(ve = $,) > 0, 
ATV, = 2(%e = we) > 0, 


therefore Aj € {1,2} and A7, = —2(¥, — ve) —2(w, — $,) < 0 


RS, (U,,U,) is a RS. In this case v; < v, S V,, Wa S W, < W = We = We, 
F(U) = q < f(U,) and Us, U, € Qe. CRSZ [Us U,] has two waves (Uz, Ü,) and 
(U,, U,) that are a NS and a PT or a CD, respectively. As a consequence 
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therefore Af = 1 and A7, = —2(v, — ve) < 0. Notice that v, > v, if and only 


if we = w, = W and v, > v; = v4. 


Assume that two waves (U;,U,,) and (U,,,U,) interact at time t > 0. Let 
U, =U,(U;,U,). Notice that U, = U, if and only if (Uz, Un) is a S or a RS, while 
U, = U; if and only if (Ue, Um) is a CD. 


e If the interaction occurs at x Z 0, then one of the following cases occurs: 


CD-S (Ur, Un) is a CD and (U4,U,) is a S. In this case vg = Um > Up = v, 
Wm = Wr, w, belongs to the closed interval between w, and w,, W(U;) = W(U.), 
Wm = Wr, f (Un) > f(U,) and U, U, € Qe. RSR[Ue, U,] has at most two waves 
(U,,U,) and (U,,U,) that are respectively either a S and a CD, or a PT and 
a possibly null CD. As a consequence 0 = ATV, = ATV, = AU. = AŤ, 
therefore Aj < 0 and A7, = 0. 


CD-RS (U;,Um) is a CD and (Um,U,) is a RS. In this case v; = Um < Ur = w, 
We = W, Wm = Wr, F(U) < FU), f(Um) < TW) and Us, Um, Us, Up € 
Qe. RSR|U:, U,] has two waves (Ur, U.) and (U.,U,) that are a RS and a CD, 
respectively. As a consequence 0 = ATV, = ATV, = AT, = AT», therefore 
At = 0 and A7, = 0. 

CD-PT (Ur, Um) is a CD and (Um, Ur) is a PT. In this case v; = Um = Vmax > 
Up = Us, Wm < W X w, W, belongs to the closed interval between w, and wy, 
U, € Qf, Un € OF and U,, U, € Qe. RSĘĄ|U,,U,] has at most two waves (Us, U.) 
and (U+, U,) that are either a PT or a S and a possibly null CD, respectively. 


Asa consequence 


ATV, = 0, AT, < 
ATW, = |we — w,| — (Iw; — Wml + |Win — v,|) < 0, AŤ, = 


therefore Aj < 0 and A7, < 0. 


S-S (Ur, Um) and (Um,U,) are Ss. In this case vg > Um > Ur, We = Wm = Wr ŻW 
and Uy, Um, Ur € Qe. RSĘ(U,U,| has one wave (Up, U,), which is a S. As 
a consequence 0 = ATV, = ATV, = AT, = AT,, therefore Aj = —1 and 
AT, = 0. 
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S-RS (U;, Um) is a S and (Um, Ur) is a RS. In this case v; > v, > Um, WŁ = Wm = 
wr > w` and Uy, Um, Ur € Qe. RSĘ[U,, U,] has one wave (Up, Up), which is a S. 


Asa consequence 


ATV, = —2(v, — Um) « 0, AT e 
ATV =0, AY. e 


therefore Aj = —1 and AT, < 0. 


RS-S (Ur, Um) is a RS and (Um, U,) is a S. In this case Um > v; > Up, We = Wm = 
wr > w` and Uy, Un, U, € Qe. RSR[Ur, U,] has one wave (Us, U,), which is a S. 


Asa consequence 


ATV, = —2(um — ve) « 0, AM = 
ATV, — 0, A < 


therefore Aj = —1 and A7, < 0. 


PT-S (Ur,Un) is a PT and (U4,U,) is a S. In this case vg = Vmax > Um > Ur, 
We < W7 X Wm = Wr, U, € XQ; and Um, U, € Qe. RS%|U,,U,| has one wave 
(U2, U,), which is a PT. As a consequence 0 = ATV, = ATV, = AT, = AŸ,, 
therefore Af = —1 and A7, = 0. 


PT-RS (U, Un) is a PT and (Um,U,) is a RS. In this case v; = Vmax Z Ur > Um; 
We < W7 X Wm = Wr, U, € OF and Um, U, € Qe. RS%|U,,U,| has one wave 
(U,,U,), which is either a PT or a CD. As a consequence 


ATV, = —2(v, — Um) < 0, AT, = 0 
ATV), =0: AT, < 0, 


therefore Aj = —1 and AT, < 0. 


e If the interaction occurs at z = 0 and (U;, U,) € D, then one of the following 


cases Occurs: 


CD-So (Ur, Un) is a CD and (Um, U,) is a S. In this case v; = Um > Up = Us, Wm = 


Wp, w, belongs to the closed interval between w, and w,, W(U;) = W(U.), um = Wr, 
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f(U,) < (Um) € q, min( f(U2), f(U.)) < q and Un, U, € Qe. CRSR[Uz, Ur] has 
at most two waves (U;,U,) and (U.,U,) that are respectively either a S and a 
CD, or a PT and a possibly null CD. As a consequence AT, <0= ATV, = 
ATV, = AT,, therefore Af < 0 and AT, < 0. 


CD-RSv (Ur, Un) is a CD and (Um,U,) is a RS. In this case ve = Um < w, 


Wm = Wr, Wx = We, f (Ue) < FU), f(Um) < TAS): max{ f (Um), f (U.)) < q and 
Ue, Um, Us, U, € Qc. CRSĘ(U„U,] has two waves (Ug, Ux) and (U.,U,) that are 


a RS and a CD, respectively. As a consequence 


> 
R> 
3 


ATV, =0, 
ATV, = 0, 


eo © 


> 

re 

3 

A IX 


therefore Aj = 0 and A7, x 0. 

CD-NSo (Ue, Um) is a CD and (Um, Ur) is a NS. In this case v; < ve = Un < Up € 
Vos WT S we S Wr < Wm, J (U) < JU) < q = f(Um) = f(Urs) and Us, Um, Ux, 
U, € Qe. CRSĘ|U,,U,] has a fan of RSs ranging from U; to U, and a possibly 
null CD (U,,U,). As a consequence ATV, = —2(wm — w) < 0 = ATV, = 
AT, = AT, therefore At € [-1,2" — 1] and A7, < 0. 

CD-PT, (U Un) is a CD and (Um, U,) is a PT. In this case v; = Um = Vmax > 
Up = Us, Wm < W X W, w, belongs to the closed interval between w, and w,, 
mint f (U,), f(U,) < q, max{f(Um), f(U,)) < a, Ue € Qr, Un € Q; and UA, 
U, € Qe. CRS[U, U,] has at most two waves (Ur, Ux) and (U.,U,) that are 


either a PT or a S and a possibly null CD, respectively. As a consequence 


ATV, ext AT, « 0, 
ATW, = |we — w,| — (Iw; — Wml + |Wm — Wl) < 0, AT, = 0, 


therefore Aj € (—1,0) and A7, < 0. 

S-So (Ur, Um) and (Un, U,) are Ss. In this case v; > Um > Up, We = Wm = Wp > W7, 
f(U,) < f(U;) € q and Ur, Um, U, € XX. CRS%[U,,U,] has one wave (U;,U,), 
which is a S. As a consequence 0 = ATV, = ATV, = AT, = AT,, therefore 
At = —1 and A7, = 0. 


S-RSo (Ur, Un) is a S and (Um, U,) is a RS. In this case v, > v, > Um, We = Wm = 
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wr 2 w`, f(U,) < f(U) < q and UZ, Um, Ur € Qe. CRSĘ(U,,U,] has one wave 


(U,,U,), which is a S. As a consequence 


ATV, S 2G, — 95) <0, AV y= 
ATV, = 0, AT, < 


therefore Aj = —1 and AT, < 0. 


RS-So (Uz, Um) is a RS and (Um, U,) is a S. In this case Um > v; > Up, We = Wm = 
w, > w7, f(U,) < f(Uc) < q and Ug, Um, U, € Qe. CRS%|U,,U,]| has one wave 


(U,,U,), which is a S. As a consequence 


ATV eae e AY cs 
ATV, — 0, ANY 


therefore Aj = —1 and AT, < 0. 

NS-So (Ur, Un) is a NS and (U4,U,) is a S. In this case um > v; 2 v; = w, 
Wi = W > Wm = Wr > w`, TIU = f(Um) = q > f(U,) and Us Um, Ur € Qe: 
CRSĘ|U,, U,] has at most two waves (Up, U.) and (U,,U,) that are a possibly 
null S and CD, respectively. As a consequence ATV, = —-2(u, — v) < 0 = 
ATV, = AT, = AT,, therefore At € {—1,0} and A7, < 0. 

NS-PTs (Ur, Un) is a NS and (U4,U,) is a PT. In this case vm = Vmax > YZ 
Up = Vy, We = Ws > WQSEW. > Wm, f > fU = f(Um) =q > TUT UT <q, 
Uy, U, € Qe and Um € OF. CRS%|U,,U,| has at most two waves (Up, U+) and 
(U., U,) that are a possibly null S and a possibly null CD (but not both null), 


respectively. As a consequence 


ATV, = —2(Vmax — ve) < 0, AT, = 
ATV, = —2(w — Wm) < 0, AT, = 


therefore Aj € (—1,0) and A7, < 0. 


PT-So (Ur, Um) is a PT and (Um,U,) is a S. In this case v; = Vmax > Um > Ur, 
We < w < Wm = Wr, f(U-;) < F(Um) < max1 f (U;), f(Um)} < q, U; e Or 
and Um, U, € QZ. CRS[U;, U,| has one wave (U,,U,), which is a PT. As a 
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consequence 0 = ATV, = ATV, = AT, = ANI therefore Af = —1 and 
AT, =0. 

PT-RSo (U, Um) is a PT and (Um,U,) is a RS. In this case v; = Vmax 2 Ur > Um; 
We < W7 X Wm = Wr, Up € OF and Um, U, € Qe. CRSZIU,,U,| has one wave 
(U,,U,), which is either a PT or a CD. As a consequence 


ATV, = 20. 94) < 0, AY 
ATV,„ =0, AŤ, < 


therefore Aj = —1 and AT, < 0. 

PT-NSo (Ur, Um) is a PT and (Um, U,) is a NS. In this case v; = Vmax 2 Ur > Um; 
WZW X Wr < Wm, f (Ue) < f(Um) = f(U,) = q, UŁ E OF and Um, U, € 
0. CRS'S|U, U,] has one wave (Up, U,), which is either a CD or a PT. As a 


consequence 


ATV, = —2(v, — Um) < 0, AM 
ATV y = — {um — w,) « 0, AT. 


therefore Aj = —1 and A7, < 0. 


e Assume that two waves (Up, Um) and (Um,U,) interact at z = 0 and (U,,U,) € 
Dy. 


If q € [f-, f^), then one of the following cases occurs: 


CD-S4 (Uę, Um) isa CD and (Um, U,) is a S. In this case Ÿ, 2 v; = Um > Up = Vp > 
Ve, We — We > Wr > Wm — Wr, W > Ür, /(U,) > (U) >q JG) > f(U,) and 
Uz, Um, Ur € Qe. CRS%[U,,U,] has three waves (Us, Uz), (Ue, Ur) and (U,,U,), 


which are a S, a NS and a CD, respectively. As a consequence 


v V 


ATV, = 2(v, — $,) > 0, AT, = —($,, — te) — (wę — tim) < 0, 
ATV, — 0, AT, = 0, 


therefore Aj = 1 and AT, = —2(¥, — v,) — 2(we — $,) < 0. 
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CD-RS4 (Uz, Um) is a CD and (U4,U,) is a RS. In this case fẹ € v; = Um = 
Öm < Vr = Ýr, Vp < Üm, We = Ñe > Üm Z Wm = Wr, Ùr < ùm, f(Um) < f(U,) S 
(Ur) < f(U.), f(Um) € q < f(U.) and Ue, Um, Ur € Qe. CRSĘ(U„U,] has at 
most three waves (Us, Üs), (05, 0.) and (U,,U,), which are a possibly null S, a 


NS and a possibly null CD, respectively. As a consequence 


ATV, = 2(ve = Ÿe) > 0, 


2(w = wr) if (Um) = i " 


ATV, = 
p if f(Um) < q 


and therefore 


Abe(-1,0,1, AT, « —2($,, — vp) + p if {(Um) = à « 0. 
-hüm — à) if (Un) <q 
CD-NSq (Ur, Um) is a CD and (Um, U,) is a NS. In this case v; = Um < w, 
Wr < MINĄ Wg, Wm}, WŚ Wm, f (Um) = JU) = q # f(Us) and Us, Um Ur E€ Qo- 
If f(U;) < q, then CRS%[U;,U,] has a fan of RSs from U; to Uy and a NS (Uy, U,); 
as a consequence ATV, = —2(wm — we) < 0 = ATV, = AT, = AŸ,, therefore 
At € [0,27 — 2] and A7, = —2(wm — we) < 0. 
If f(Uj) > q, then e < v, = Um = Wm, We = We > Wm = Wm, CRSRÎUz, Ur] 


has a two waves (U,,U,) and (U,,U,), that are a S and a NS, respectively; as a 


consequence 
ATV, = 2(ve = Ve) > 0, AT, = — (ve = vo) = (we = Wm) < 0, 
ATV, — 0, AŤ, = 0, 


therefore Aj = 0 and AT, = —2(w — Wm) < 0. 
CD-PT; (Up, Um) is a CD and (U4,U,) is a PT. In this case wm <w <w = 


Ùm < We, Dr = Um exu > Ur = M, > Fp, Wm < W7 S W, < W < We = Ñy, 


M 
F(U) > FU) > q > max{f(Um), f(U,)}, Ue € OF, Un € OF and U, € Qe. 
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CRSĘ|U,, U,] has three waves (Ue, Uz), (Ue, Ur) and (U,,U,) that are a S, a NS 


and a CD, respectively. As a consequence 


ATV, = Av, — Ÿe) > 0, AT, = — (max — 9) — (we — Wy) < 0, 
ATV,-—-—2(w.—w.)«0, AT, =0; 


therefore Aj = 1 and A7, < —2(Vmax — Ur) — Z(wę — w) < 0. 


NS-S4 (Uz, Um) is a NS and (Um, Ur) is a S. In this case Um > v, = V, > w, 
We > ùr > Wm = Wr > w 2 w7, J(U) > q = TU = TU: > f(U,) and 
Ur, Um, U, € Qe. CRSR|U:, U,] has two waves (U;,U,) and (U,,U,) that are 
a NS and a CD, respectively. As a consequence ATV, = —2(Um — vy) < 0 = 
ATV, = AT, = AT,, therefore At = 0 and AT, = —2(um — v,) < 0. 

NS-RS4 (Uz, Um) is a NS and (Um,U,) is a RS. In this case v; < Um = V < Up = 
Ür, We > Wm = Ùm = Wr > ùp, f(u) TL) > q = Fy = f(Um) and Ux Um, 
U, € Qe. CRSR|U;, U,] has two waves (Up, ,) and (U,,U,) that are a NS and a 


CD, respectively. As a consequence 


ATV, — 0, AT, = 0, 
ATV, = 2(w, — w,) > 0, AT. = — (v, — Um) we < 0, 


therefore Aj = 0 and A7, = —2(v, — Um) < 0. 
If q € [0, f7), then one of the following cases occurs: 


CD-Są (Ur, Um) is a CD and (Um, U,) is a S. In this case $4, 2 Ve = Um 
ŻA We = WYZN > Wm = m = W, > w7, TOO US) >q È f(Um) > 
f(U,) and Ur, Um, Ur € Qc. CRS%[U,, U,] has three waves (Us, Uz), (Ue, Ü,) and 
(U,, U,), which are a S, a NS and a CD, respectively. As a consequence 


ATV, = 2(v, — %) > 0, AT = — (fm wes (we — Wm) < 0, 
ATV, — 0, AT, = 0, 


therefore Aj = 1 and A7, = —2(%m — v.) — 2(w; — Wm) < 0. 
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CD-RS, (Ur, Um) is a CD and (Um, U») is a RS. In this case $; € ve = Um S Vm, 
Um Van Ur S Veg We = e Z Ùm Z Wm = Wr > WT, We > Wm, f(Um) < f(U,) S 
F(U) < f(Ux), f(Um) € q < f(U.) and Ue, Um, Ur € Qe. CRSĘIU„U,] has at 
most three waves (U,,U,), (05, 0.) and (U,,U,), which are a possibly null S, a 
NS and a possibly null CD or PT, respectively. As a consequence 


ATV, = 2(v, — Ve) + 2(¥, — vr) 20, 


(w, — wr) if f(Um) = si 


ATV, = 
| if f(Um) <q 


and therefore 


At € {-1,0,1}, AT, = —2($4, — w) — 2(v, — Um) — 2(We — Wm) 
E f nae (CP à 
2,8) if (Un) <q 

CD-NS, (Uz, Um) is a CD and (Um, Ur) is a NS. In this case vg = Um = Vm < Ur, 

w, < min[w;, Wm}, We = Wy, Wm = Am, f(Um) = f(U-) = q # f(Ue) and UW, 

Um € Qe. 

If f(Uc) < q, then v, > Ÿe > w = Um, Wr < We < Wm and CRSZ[U„U,] 

has a fan of RSs from U; to Ü and a NS (Ge, U,); as a consequence ATV, = 

—2(w, — we) < 0 = ATV, = AT, = AYn, therefore At € [0,2” — 2] and 

AT = —2(Wm — we) < 0. 


If f(Uz) > q, then v, > v; = Um = tm > Ve, We > Wm > Wr and CRSR[Uz, Ur] 


has a two waves (U;,U,) and (U;, U,), that are a S and a NS, respectively; as a 


consequence 
ATV, = 2(ve = 44) > 0, AT, = —(ve — Ÿe) = (we = Wm) < 0, 
ATV, = 0, AŸ, = 0, 
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therefore Af = 0 and A7, = —2(we — Wm) < 0. 


CD-PT, (Uz, Um) is a CD and (Um, Ur) is a PT. In this case v; = Um = (m 
Vmax > Ur = Vr > Ve, We = We > Ww 2 Ur ŻW > W = Wm È Wm, 
(U) > f(U) > q > max{f(Um), f(U,)), Ur € OF, Un € OF and U, € €x. 
CRS%[U,, U,] has at most three waves (U;, U;), (U,,U,.) and (U,,U,) that are a 


S, a NS and a possibly null CD, respectively. As a consequence 


ATV, = 2(v, — ĉe) > 0, AT, = —(Um — 9j) — (we — wą) < 0, 
ATV, = —2(u, — Wm) < 0, AŤ„=0, 


therefore Aj € {0,1} and A7, = —2(v« — vr) — 2(w, — Wm) — 2(w; — w) < 0. 

NS-S; (Ue, Um) is a NS and (Um, U,) is a S. In this case vum > v, = V, > w, 
wW Me > Wm = Wr > w7, JU) = f(Um) = q > fU. and Us Um, Ur E€ Qe: 
CRS%|U,,U,] has two waves (U;,U,) and (U,,U,) that are a NS and a CD, 
respectively. As a consequence ATV, = —2(um — v,) < 0 = ATV, = AT, = 
AŸ,, therefore Af = 0 and AT, = —2(Um — v,) < 0. 

NS-RS, (U;, Um) is a NS and (Um, U,) is a RS. In this case ve < Um = V, < Up € 
Vus WE > Wm = Ùm = Wr > MT) =] (Ue) = f(Um) and Ug, Um, Ur E€ Qe- 
CSS": U;, U,] has two waves (U;, Ù,) and (U,, U,„) that are a NS and either a PT 


or a CD, respectively. As a consequence 


ATV, = 2(ÿ, — v,) 20, AT, = 0, 
ATN = 2(w, — ù) > 0, AT, = —(V, — Um) — (w, — $9) < 0, 


therefore Aj = 0 and A7, = —2(v, — um) < 0. Notice that v, > v, if and only 
if wm = Wr = w^ and ©, = Vmax > Ur > Um = Vd. 

NS-P'T, (U:, Um) is a NS and (Um, Ur) is a PT. In this case Um = Vmax > Ur = 
Vr > Ve, We > Ur > Wr = W > Wm, f(Us) > f(Ue) = f(Um) = q > f(Ur), Ue, 
U, € Q. and Um € QF. CRSR|U;, U,] has two waves (U,,U,) and (U,,U,) that 


are a NS and a CD, respectively. As a consequence 


ATV, = —2(Vmax — vr) < 0, ATU, 
ATV, = —2(w^ — Wm) < 0, AŤ, = 
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therefore Aj = 0 and A7, = —2(Vmax — Ur) — 2(w^ — Wm) < 0. 


This concludes the proof. 


A.6 Proof of Proposition 5.8 


We underline that for any interaction Af < 2” — 1, see Table A.1. 
From what we already show in the proof of Proposition 5.7, see Table A.1, we 
deduce that 
te a + (t) 


strictly decreases after any interaction, except the following cases. 


A CD (Ur, Um) interacts with (Um, U,) and one of the following con- 


ditions is satisfied: 


e (Um, Ur) is a S and wọ —w — 1; 

l (A.1) 
e (Um, U,) is a S and w™ — 1 < we S WT = w; 
e (Un, Ur) is a RS; 

e (U„,U,)isaPTandw,>w.. 


For this reason it remains to bound the number of only the above type of interac- 
tions. We observe that the number of waves of U, do not change after interactions 
as in (A.1). This implies that the number of waves is uniformly bounded. We also 
observe that any interaction as in (A.1) has exactly one incoming CD and exactly 
one outgoing CD. Since no wave can reach any CD from the left (and then possibly 
have with it an interaction as in (A.1)), we have that as long as a CD remains a 
CD (possible further interactions involving it have to be taken into account), it 
can interact only once with another wave W (or with waves generated by further 
interactions involving W), moreover in this case W is slower then such CD and is 
not another CD. Since furthermore we already know that the number of waves is 
uniformly bounded, there can be only finitely many interactions involving CDs. It 
is therefore now clear that also the number of the interactions described in (A.1) 


is bounded. 


Nikodem Dymski 


Technical proofs of Subsection 4.3.1 


In this section we collect the proofs concerning the non-conservative constrained 


ARZ model. 


B.1 Proof of Lemma 4.1 


It is enough to show that the function h,(v) = v + p(q/v) is Lipschitz and convex 


in the interval I = [6(wy), v’(wn)] C (0, +20). For any v € I we have 
q (d q QY , d. n(4 
POS A e m) v3 (av() | tp (2)). 
By (2.12) we have p'(p) > 0 and 2p'(p) + pp”(p) > 0. As a result h7(v) > 0, so 


NEC 


and therefore h, is Lipschitz. Furthermore the inverse of hg restricted to intervals 
[6(wn), 6(w?)] and [v’(w®),v’(wy)] are Lipschitz. Notice that 6(w?) = 5(q) = 


v? (a4). 
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B.2 Proof of Lemma 4.2 


We denote by CD a contact discontinuity, R a rarefaction shock, S a shock and NS 


a non-classical shock. For completeness, we compute below 


A, = TV (v(t?)) — TV (v(t; )), A, = TV (w(t?) — TV (w(t;)), 


A; = TV (6(tż);R_) — TV (@(t7);R_), 
AS TV (v (tT R-) = TV (v (6  R-); 
AU = TV (w(i; ); R_) CIV (w(t; ); R ) , 
Ay= (tH) - 77), ATS rE +r) — AËZ HI) — NG), 
and 
A= >. BU ue SU ue) nios SU t) 
rEJ(tj) z€J(t; ) 


A Assume now that at time t; > 0 exactly one interaction between two waves 
(Ur, Um) and (Um,U,) occurs at x = 0. For notational simplicity let f; = f(U), 
de = 5(wę), v) = v (wi), U =U(UU,) = (we, v.) and so on. We distinguish the 
following cases: 

A1 Assume that (U, Un) is a S and (U,,,U,) an NS. Clearly we = Wm 2 w, 
Um < min {ve,v,} and fr < fm = fr = q. As a consequence (Up, Up) € Qi, see 
Figure B.1. 

Al.a If v, > vw, then At > 0 because RSI, U,] has à fan of Rs (Us, U) and 
a CD (U,U,), see Figure B.1. Obviously 


Ac Rh A, = 0. 


b 


m; we have 


Since Un € Ù < Ve < 6 = Ur, à X dj = We = Wm € w, and v, > Ve = v 
Ay = —2(v, — Um) < 0, Ay = —(v? — Um) < 0, Aw = 0, A =0. 
In conclusion we have AT = —2(v; — v?,) < —2e. 


A1.b If v, € ve, then At € 0 because RS (Ue, U,] has a possibly null S (U;, U) 
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Figure B.1: Case A1. 


and a possibly null CD (U,U,) (but not both null), see Figure B.1. Obviously 


0=A,=A;=A,=A 


UE 


b 


m; we have 


Since Um < 0 = v, < vj, t) X Ÿ = We = Wm, Ù < w, and v) = v 
Ay = —2(vr — Um) < 0, Aw = 0, Ay = — (min [vv] — Um) <0. 
In conclusion we have 


AT = —2 (vr — min fu, u, |) <0. 


A2 Assume that (U;, Um) is a CD and (Um, U,) an NS. Clearly vg = Um = Om < Up, 
Ùw < min{ Wm, wand fm = fr = q. 

A2.a If (U„U,) € Qı, then Af > —1 because RSẸ2,[U;, U,] has a fan of Rs 
(U;, U) and a possibly null CD (U,U,), see Figure B.2. Since 6 = v, > v; = Um, 


We = Wm and u(v,) = w, we have 
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Since U = we < Wm, We € w, and W € min{wm, w, |) we have that 
Aw = (Wr — Wp) — |w, — Wml] < 0, A, = —(wm — We) < 0. 
Since wę < Wm and fm — q, we have that U, > v > Ù 2 ty 2 Um = Um we have 


Az = — (ôe — Um) < 0, AY = —(v 


sf) 
m 
Ay = — (min tv, oe] — Um) <0. 
In conclusion we have that if v, € (Um, 0] then 0; > vp, we € w, < Wm and 
AT < —w,, — 2w, + 3we — 2 (6; — min (v2, v, }) — 2 (v, — v) < —3e, 
while if v, > © then wy AW, We < Wm, Ôp = Ù = v and 


AT < —w,,—2w43w,—2 (u, — min [v5 v. ]) X —e. 


f f 


U, U, 
28 Un 7 : 
| i | 
p p 
R CD R CD 
U 0 
U, U, 
Um 
CD NS 


Figure B.2: Case A2.a. 


A2.b If (Up, U,) € Q and à < we < Wm, then Af > 0 because RS%pz[U;, U,] has 
a fan of Rs (Uz, 0) and a NS (U, U,), see Figure B.3, left. Since v, > tg > v; = Um 
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and w(t) = De = we < we = Wm, we have 


Since we = We < Wm we have that 
A, = |w, — Wel — ((Wm — we) + [wr — wml) <0, A, = —(wm — we) < 0. 


Since uU < w < Wm and fm = q, we have that v, > vb > Ôp > Ve = Um = Vm we 


have 


A = — (ĉe — ve) < 0, A = c (in — vt) <0, 


In conclusion we have that Af > 0 and 


AT < 343 + 2(4, + 45 FA, — 47) 
< [us — we) e 2(0). =) + 2 (min TA v.) — min DA v, |) 
X — (Wm — We) € —E. 


U, 
p 
R NS S R NS S 
Ü, U, 
U, U, U U, U, U, 
Um Um Un 
CD NS CD NS S S 


Figure B.3: Cases A2.b, A2.c and A3. 


A2.c If (U„U,) € Q2 and w < Wm < we, then At = 0 because RS%h,[U;, U,] has 
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a S (Ur, Up) and a NS (Op, U,), see Figure B.3, right. Since v, > Um = v, > te and 


no R is involved, we have 
Ay = 2(v, — be) > 0, A^ = 0. 
Since Wm < We = We we have that 
Aw = |w, — Wel — ((we — Wm) + fur — wml) <0, A; = —(we-— wm) < 0. 


Since w < Wm < w, and fm = q, we have that UŁ > u, > Un = tm = Ve > 0, we 
have 


A; ——(vw—9) «0,  Aj=-fluż— vh) <0, 


v 


Ay = (ve — de) + (min {up vr} — min CA) > 0. 


In conclusion we have that Ag = 0 and AT < A, + 2A; = 0. 
A3 Assume that (U, Um) and (Um, U,) are two shocks. Clearly v, < Um < w, 
We = Wm = w, and max{ fe, fm, fr} € q. Hence, the solution consists of a shock 


(U,,U,), see Figure B.3. In this case we have 


0 =A, = ^ = 45 = 45 = 4y = A7 = A, = AT. 


v 


A4 Assume that (U;, Um) is a CD and (Um, U,) is a S. Clearly Wm = Wr, v; = 
Um > U, = Ù and max{ fm, fr} € q- 

A4.a If (U,,U,) € Qı, ie. min{ fe, f} < q, then Af = 0 because RS Ur Ur] 
has a S (U2, Ü) and a CD (© U,.), see Figure B.4. Since v, = Un > v, = 0, We = Ù, 


Wm = w, and neither R or NS are involved, we have 


(SA, 54, = = OF, A, = —|we — w,| < 0, 


w 


A; = —[6, — 6,| < 0, Ap = lu; — | < 0. 


In conclusion we have that Af = 0 and AT < 0. 
A4.b If (U,,U,) € Qa, namely min{ fe, f} > q, then we = W > Wm = Wr, Ur > 6, 
and Af > 0 because RS (Us, U,] has a S (Up, Ue), a NS (DU) and a possibly 


Nikodem Dymski 


Section B.2 Page 147 


Figure B.4: Cases A4.a and A4.b. 


null CD (U,,U,), see Figure B.4. Since wm = wp < dy < we = d, 0; < 0 — v, = 


Dp < Up = Um < v and no R is involved, we have ty < 0, < v? « v and 


— (0, — de) <0, A, = — (04 — v?) <0, 


A 2 A; = 
Aw = 0, A, = —(we — Wm) < 0, Ay = (v, — dy) > 0, 
0 


In conclusion we have that Aq > 0 and AT < A, + 2(A; — Ay) = —2(6, — ĉe) < 
—2€. 

A5 Assume that (Ur, Um) is a NS and (U,U,) is a S. Clearly wm = W, Um > 
max (ve, v.) and f, < fm =q = fe. 

A5.a If (U„U,) € Qi, then v, < w, max{f,,f} < fe = q = fm. Moreover 
At < 0 because RSG, (Us, U,] has a S (Uz, U) and a possibly null CD (U, U,.), see 
Figure B.5. Clearly 


0 =A; = A3 = A3 = An 


Since v, = 0 < vg < Um, Ve X v, we = W and wy = w,, we have 
A, = —2(um — ve) < 0, Ay = — [min{u?, Um} — vel « 0, Au =Q. 
In conclusion we have that Af < 0 and AT = —2 [Um — minfv;, Um } | <0. 
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Figure B.5: Cases A5.a, A5.b and A6. 


A55 If (Up, D.) E€ Os, then ve < Ur, fr < fe So = In < f. Moreover Af = 0 
because RS [Ur U,] has a NS (U, U.) and a CD (Ur, U,.), see Figure B.5. Clearly 


0 = 47 = A3 = A3 = Ar 


Since uj < Up = Up < Um < v and Wm = w, < w, < wy, we have 
A, = —2(um — Ur) < 0, Nem T Ay = —(Um — Ur) < 0. 


In conclusion we have that Aj = 0 and AT — 0. 

A6 Assume that (U, U„) is a NS and (U,,,U,) a R. Clearly w, = w, < w, 
Um < v, — € and q = fe = fm < fr. Moreover Af = 0 because RS og [U,, U,] has a 
NS (U,,U,) and a CD (U,,U,), see Figure B.5. Clearly 


0 = A5 =A, = Au 
Since v > ©, = Ur > Um > Ve and we > Wp = Wm = W, > Wy, we have 
A,=0, Ay = 2(w, — dr) > 0, Ay = (up — Um) > 0, A, = — (w, — à) < 0. 


In conclusion we have that Af = 0 and AT = —2(v, — um) < —2e. 
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A7 Assume that (Us, Um) is a R and (Um, Ur) a S. Clearly we = wp = w, and 
ve > Um — €. It is easy to check that (Up, U,) € Q4. Thus RSS, [Ur, U,] is a S 
(Us, U,) and Af = —1, see Figure B.6. Clearly 

0=A, =Az =Az =Az=Ay, A, 


vU 


A 
© 


Since Un > Ve > v, we have 
A, = —2(Um — ve) <0. 


In conclusion we have that Aj = —1 and AT < 0. 


Figure B.6: Cases A7, A8.a and A8.b. 


A8 Assume that (Ur, Um) is a CD and (Um, Up) is a R. Clearly v; = um < w, 
Wm = Wr and mini fm, fr} <S q. 

A8.a If (U„U,) € Q1, namely max{ fe, f) < q, then R55 [U;, U,] has a single R 
(Us, U) and a CD (U,U,), see Figure B.6. Clearly 


A= NA 81] fu>0l 4 A- -u — vh] ify >0 m" 
0 if v, = 0 ° 0 if v; = 0 


À —|we—Wm| if ve > 0 <0, Ay=0. 
0 if ve = 0 
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Since Ve = Um < Up = 0, W = we and Wm = Wr, we have 
Ng = 0; AWZ A^ =: 


In conclusion we have that Ag = 0 and AT < 0. 

A8.b If (Us, U,) € Qa, i.e. max( fe, f} > q, then Af > —1 because S55; Us, U,] 
has a possibly null S (Uz, Us), a NS (Ue, U,) and a possibly null CD (U,, U,), see 
Figure B.6. Since W, > w, — €, Wm = Wr < We = We, 0; X; Ve = Ug < Vp = Up = 


0 € v), 0; < bm = 0, and vè >Y? =v", we have 
£3 £ m T? 


2 p — Dy. if m= |< jr 

A= 2(v — be) 20, A= he ZĘ 
0 otherwise 

A;— -—(0,—9, «0, Aj=—(lw — w) < 0, 


ro Y pus Ti : “ < < . 
A, = —(u; = v) « 0, A= (Wm wr) if 6(Wmax) < Um < Up SV Z 0, 
0 otherwise 


Ay= (Ur = de) > 0. 


In conclusion we have that Af > —1 and 


AT < A, — 2A7+ Av +2A, 


2(w, — Wm) =0 iffn=q<f, 
Lex gua (Wr — Wm) nee ile a 
0 otherwise 


A9 Assume that (U2, Um) is a S and (Um, Up) is a R. Clearly Um < v, < w and 
We = Wm = w,. By observing that f(RSarz[U2, U,])(0*) = min{ fe, f.) < q, we 
have that (U,,U,) € Q1. Hence RS lUe, U,| = RSaRz(U, U,] has a S (U,,U,) 


and Aj = —1, see Figure B.7. As a consequence A, < 0, Ay = 0 because no 
NS is involved, A,, = A; = A; = A’, = 0 because no CD is involved. Since 
Um < Um + € € Up < v, we have A, = —2(v, — Um) < 0. In conclusion Af = —1 
and AT « 0. 


B Assume now that at time t; > 0 exactly one wave (U,, U,) reaches z = 0. We 
distinguish the following cases: 

B1 Assume (U;, U,) is a CD. Clearly v, = v, and f, <q. 

Bl.a If (U4,U,) € (4, namely f; € q, then At = 0 because RSS lU U,] has a 
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CD (U4,U,), see Figure B.7, left. Clearly 0 = A, = A, = Ay = A, and 


Ap = -lw -w| S-E, A7 = -lê ô| < e, Ag = loe- vrl S E, 


therefore AT = 3A; +2(A; + A) < —7e « 0. 


f F f 
q z : 
1 U, 
U, q à 
p p 
CD S S 
U; / U; U, 
cD S S 


Figure B.7: Cases A9, Bl.a, B1.b and B2. 


B1.b If (U„U,) € Q, namely f; > q, then Aż > 0 because RS (Us, U,] has 
a shock (Us, Us), a NS (U„U,) and a possibly null CD (U,,U,), see Figure B.7, 


center. Since v, = v, = 6, > 0, and no R is involved, we have 

A, = 2(v, — de) > 0, A, D. 
Since w, < tb, < we = We we have that 

scm 0: A, = —(we — w,) < 0. 
Since w, « wy, we have that v > v? SUD > Ün Ü= s = dy vU and 


A; = —(t, — 6%) < 0, A7 = —(v — v?) < 0, Ay = v — 6, > 0. 
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In conclusion we have that Af > 0 and 
AT = A,+2(A; ŁA, — Ay) +3A, = —2(v) — v?) — 2(0, — de) —3(we—w,) < —Te. 


B2 Assume that (U;,U,) is a S. Clearly we = w, and max{ fe, fr} < q. Hence, 
At = 0 because RS [U,, U,] has a S (U,,U,), see Figure B.7, right. In this case 


we have 


v 


At=0=A,=A, =A; = 45 = Ay= A} = A23 = AT. 


B3 Assume that (U;,U,) is a R. Clearly we = w, and v; Ś v, — €. 


f f 


f 
U, 
Z: AIDS u 
q U, 
p 
R 


Figure B.8: Cases B3.a, B3.b and B3.c. 


B3.a If (U,,U,) € Q1, namely max { fe, fr} < q, then RE zip, U,] has a single 


R (Us, U,), see Figure B.8, center. In this case we have 


Af=0=A,=A,=A;=A,=Ay=A,=A,=AT. 


vU 


B3.b If (U,,U,) € Q9 and f, =q < fe, then Aj = 1 because RS Te; U,] has 
a S (U,, U.) and a NS (U De), see Figure B.8, right. Since w = W, = w, and 
bp < Ù < v, < v, = V}, we have 


A, (uy — 0) > 0, 0m A, mA; =A, =A, =A, Aye (v0): 0. 


vU 
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In conclusion we have Af = 1 and AT = —2(v, — ue) < —2e. 
B3.c If (U,,U,) € Qz and f, = q < fr, then Aj = 1 because RS [Us U, has 
a NS (Uy, U.) and a CD (U); see Figure B.6. Since w, < w = we = we and 


Ve < Up = Ur < v, we have 


0=A,=A;, =A, = Ay, A, = —(we — i) < 0, 
Ay = (v, — u) > 0, Ap = 20h — Ùr) AS. 


In conclusion we have Aj = 1 and AT = —2(v, — ug) € —2e. 


C At last, assume that at time t; > 0 exactly one interaction between two waves 
occurs at x # 0. In this case A, < 0, the number of the waves does not increase, 
At < 0, the size of the jumps in the w-coordinate does not change, hence 0 = 
Ay = A A7 Az; no R is created, hence A, < 0, and clearly no NSs are 


w v 


involved, hence Ay. As a consequence Af < 0 and AT < 0. 
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